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FOREWORD 


This document Is submitted by the McDonnell Douglas Astronautics Company to the 
National Aeronautics and Space Administration, George C. Marshall Space Flight 
Center, In accordance with the requirements of contract NAS8-27571. Results 
are presented of a program involving the analysis, design, fabrication, and 
test of a multipurpose full-scale liquid hydrogen acquisition and thermal 
control system, termed the Interface Demonstration Unit (IDU). The IDU was 
developed for inclusion in a 105- in. -diameter liquid hydrogen tank as part of 
the NASA/MSFC auxiliary propulsion system breadboard. 

This program was sponsored by the Fluid Mechanics and Dynamics Branch, 

Propulsion and Thermodynamic Division, Astronautics Laboratories, NASA/MSFC 
under the direction of Mr. George M. Young, III. The initial phase of this 
program was performed under the technical direction of Mr. Leon J. Hastings, 
NASA/MSFC. 

This program was supported by the following personnel. J. N. Castle and 
P. Maruschak assisted in the system design and fabrication. D. W. Kendle and 
Professor Ivan Catton (Consultant) assisted in the analysis of the thermody- 
namic vent system. The computer codes were developed by D. W. Kendle. 

W. N. Geiger was instrumental in the fabrication of the IDU. Assembly and test- 
ing were performed under the direction of Roger Yeaman. The contributions and 
suggestions made by G. W. Burge, E. C. Cady, C. R. Easton, and C. E. Schroeder 
are appreciated. The support of Mr. Ralph Adams (NASA/MSFC Astrionics Labora- 
tory) in the area of capacitance-level sensor evaluation is gratefully 
acknowledged. 
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Section 1 

INTRODUCTION AND SUMMARY 


Future space missions will require cryogenic fluid storage and expulsion 
subsystems capable of providing efficient long-term subcritical storage, 
reliable multiple low-g restarts, and low-g liquid expulsion for auxiliary 
propulsion, life support systems, and in-orbit propellant transfer. These 
requirements demand significant technology improvements in the areas of 
cryogenic propellant acquisition, venting, and thermal control relative to 
the concepts used in present operational vehicles. For most applications 
requiring low-g engine restarts, it is imperative that a passive propel- 
lant control system be used which supplies 100-percent liquid under a 
wide range of adverse acceleration levels at any point in the long-duration 
low-g coast. Use of present techniques for low-g engine restarts, such as 
auxiliary propulsion system (APS) propellant settling or idle-mode engine 
restart, involves high weight penalties and may not be sufficiently flexible 
to meet advanced mission requirements. Low-g thermal control techniques, 
involving the venting of ullage gas from a stratified cryogen and requiring 
either frequent or continuous propellant settling, also result in weight 
penalties and decreases in the overall mission time. 

Presently, screen retention devices are considered prime contenders for 
cryogenic propellant acquisition, because of advantages of low weight, 
multicycle capability, reusability, long life, and simplicity. Efforts by 
the NASA and aerospace companies have resulted in a variety of screen reten- 
tion device designs and an extensive background in terms of demonstration 
and development tests, and fabrication experience. However, cryogenic 
propellant screen retention devices can be subject to certain fluid dynamic 
and thermodynamic failure modes due to transient flow, pressure decay, or 
heat- transfer-induced "breakdown".* As a result, a concept termed the 


*Loss of screen retention of liquid. 



"Integrated start tank" was conceived at McDonnell Douglas Astronautics 
Company (MDAC) In July 1968 (References 1, 2, and 3) which was designed to 
avoid problems associated with startup and shutdown transients and to Isolate 
the propellant from tank pressure fluctuations and heat transfer. 

The start-tank concept Involves a screen retention device Installed within 
a separately pressurized propellant tank located In the main tank. A line 
with a shutoff valve transmits LHg from the retention screen directly to 
the pump. A second line with a shutoff valve delivers propellant from 
the main tank to the Internal tank. In low g, when liquid In the main tank 
Is unsettled, cryogenic propellant Is supplied to the pump from the acquisi- 
tion device. When the liquid In the main tank Is settled, following some 
period of relatively high positive vehicle acceleration, (fcO.l g e ) the retention 
device can be refilled by opening the Internal tank refill valve and vent- 
ing the ullage gas In this tank, yet maintaining constant pressure in the 
main tank. In this manner, the start tank can be refilled while a constant 
flow of propellant is maintained to the pump and engine. 

Cryogenic screen acquisition devices subjected to heat transfer or tank 
pressure variations (such as pressure collapse below the saturation pressure 
of the retained liquid after engine shutdown) face potential failure modes 
due to bubble generation within the device and capillary "breakdown" due 
to screen drying. The start-tank concept minimizes these problems because 
the retained propellant is pressurized with cold helium, the tank is 
pressure isolated from main-tank pressure fluctuations, and the heat transfer 
into the start tank can be reduced by insulation, or eliminated by a wall- 
mounted thermodynamic vent system, which would intercept all heat flow into 
the start tank. 

Retention failure ("breakdown"), due to startup or shutdown feedline pressure 
transients. Is eliminated whenever the start tank is filled or the screen 
device Itself is submerged within the propellant. 
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The IDU as designed, fabricated, and tested under this contract, incorporated 
the above features in terms of propellant acquisition and thermal control. In 
addition, It serves as a full-scale test article having a wide range of 
capabilities for Investigating cryogenic propellant screen acquisition system 
operational modes. The IDU will also be part of planned NASA/MSFC liquid 
hydrogen APS breadboard tests, which will investigate the operation of a 
combined turbopump, accumulator, tankage, and acquisition system. 

To meet these multiple objectives, the IDU was designed to supply liquid 
hydrogen to a turbopump system at flowrates up to 7.5 lb/sec for a period of 
approximately 5 seconds before refill. Refill is accomplished in less than 
10 seconds, during which a constant flowrate can be maintained to the pumps. 

3 

The IDU shown in Figure 1-1 is comprised of a 10-ft cylindrical tank, 29 inches 
in diameter and 35 inches in length, with three Parker submerged liquid hydro- 
gen 2-inch ball valves, an externally mounted thermodynamic vent system with 
coolant flow control valves and orifices providing flowrates of at least 
0.2 to 5 Ib/hr, and an annular cylindrical screen device 24 inches in diameter 
and 24 inches in length, with 250 x 1370 fine mesh stainless steel screen 
within the tank. A bypass valve is provided to allow expulsion directly from 
the 105-inch diameter tank. Instrumentation provided includes seven carbon- 
point liquid-level sensors, three capacitance-point liquid-level sensors, a 
full-length capacitance liquid-level sensor, twenty platinum resistance 
temperature transducers, and four heat-flux gauges. 

Subsystem acceptance and functional tests, including acquisition subsystem 
expulsion, pressurization, and refill, were performed in the MDAC cryogenic 
test facility. Additional tests were performed to investigate operational 
modes and possible limitations of surface tension screen-acquisition devices. 
These tests included: thermodynamic vent system operation, warm-gas pressuri- 

zation effects, hydraulic pressure surge effects, screen device operational 
limitations, feedline vapor control, and two-phase refill. However, due to 
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repeated malfunctions of the Parker submerged liquid hydrogen valves, test 
results were limited and certain other tests were not accomplished. 

The tests achieved the following results: 

• The screen device was shown to perform flawlessly for the cold-gas 
pressurization design condition with a higher "bubble point" or 
breakdown pressure difference than expected and with no indication 
of degradation throughout the test program. 

• Both the overhead and submerged pressurization diffusers were 
used successfully. 

t The data acquisition and instrumentation functioned properly. 

• Flowrate control of the thermodynamic vent system was achieved 
over the entire range, and agreement between predicted and actual 
flowrates was good. 

t Nc component failures were observed, although the 2-inch submerged 
liquid hydrogen ball valves malfunctioned, as discussed in detail 
in this report. 

• The complete operational mode associated with a simulated low-g 
engine restart involving IDU outflow, followed by refill with 
simultaneous outflow was demonstrated. 

• Breakdown data were obtained, showing a marked decrease in retention 
capability when warm hydrogen pressurant was used. Warm helium 
also decreased the retention capability causing premature breakdown, 
but the effect was not as great as with warm hydrogen. 

• Use of a polyurethane adhesive to seal minor holes in the screen 
was successful; a conducting epoxy used to fill gaps in the TVS 
fillet on the IDU tank wall showed no adverse effect due to 
exposure to LHg and repeated thermal cycles. 
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• Finally, the integrity of the IDU subsystem in terms of operating 
pressures, flowrates, and soundness of welds, seals, and connections 
was demonstrated. 

Analyses were performed of the liquid flow losses, autogenous (warm gaseous 
hydrogen) pressurization, and the thermodynamic vent system; a mathematical 
model of the thermodynamic vent system (TVS) was developed, programmed, 
utilized, and documented; and an existing MDAC tank-pressurization computer 
code was modified and used in assessing autogenous pressurization effects in 
the IDU. 

Data obtained from the tests has been analyzed and the results correlated. 
However, because of feedline valve malfunctions, results in some tests are 
of marginal utility as discussed for each specific test in the following report. 
A corrective modification to the feedline valves has been recommended, based 
on a preliminary investigation of their operational factors of safety.* 

An improvement to the basic start-tank concept was conceived but has not yet 
been reduced to practice. The principle of the "jet pump" can be used 
during start-tank refill to induce gas flow out of the start tank, with 
concurrent liquid inflow, while the main tank is pressurized and while a 
net flow to the turbopumps from the main tank through the start tank is 
sustained. This concept completely eliminates the need for overboard venting 
of start-tank ullage gas during refill, thereby decreasing the total system 
weight. The concept deserves additional consideration not only because of 
the improvements to the start-tank principle, but also because of possible 
additional applications to screen-device refill for both cryogenic and storable 
propellant systems. The external plumbing of the IDU pressurization system 
could be modified to incorporate a jet pump for refill witnout overboard 
venting. 

An improved version of the cone screen housing, used for feedline vapor control, 
is described in Section 4, Test T; this version will ensure a liquid-filled 
feedline, thus diminishing the probability of vapor ingestion in the turbopump. 


*Operational factor of safety is the ratio of applied actuation pressure to 
minimum actuation pressure. 
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This program was performed concurrently with another NASA/MSFC sponsored 
program. Contract NAS8-27685 (Reference 4) with MDAC; these two programs were 
compl ementary because: (1) fabrication experience gained on the IDU was 

beneficial in the fabrication of full-scale screen channel devices; (2) a screen 
device, developed and fabricated to simulate propellant retention for a cryo- 
genic advanced space propulsion module (ASPM), is mounted inside the IDU for 
full-scale tests in the NASA/MSFC 105-inch tank; and (3) common, integrated 
feedline plumbing will be utilized between the IDU (with either the annular 
screen device or the ASPM device) and the full-scale screen channel for the 
NASA/MSFC APS breadboard tests. Finally, the IDU tests, verifying the 
operational aspects of the "start tank" concept, demonstrate the validity 
of the localized pressure isolated channel (LPIC) acquisition concept as well 
as the ASPM propellant acquisition concept selected in the system design study 
of Contract NAS8-27685. 

Section 2 of this report encompasses a description of the IDU and TVS assembly 
and installation. Section 3 presents the supporting analyses. Section 4 is 
comprised of the test description and results; supporting annotated data from 
the tests is presented in a Supplemental Data Document. Conclusions are 
presented in Section 5. Appendices support the technical discussion. 
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Section 2 

INTERFACE DEMONSTRATION UNIT DESIGN AND ASSEMBLY 


2.1 DESIGN 

This task included the design of the screen device, tank, and integrated system 
consisting of valves, components, instrumentation, line routing, and the thermo- 
dynamic vent system. Each of these elements is discussed in this subsection. 

2.1.1 Screen Device Assembly 

The basic requirements of the screen device assembly design are that it have 
an annular configuration and that it be capable of sustaining an outflow rate 
of 7.5 lb/sec of liquid hydrogen for approximately 5 seconds against 1 g before 
initiating refill. 

Details of the screen assembly design were established using an MDAC computer 
code (Subsection 3.3) which determines, for specific configurations, screen 
breakdown level as a function of flowrate and includes all pertinent flow 
losses. Results showed that a 250 x 1370 stainless steel screen with a backup 
plate and coarse mesh standoff screen would meet the necessary requirements, 
and that the basic configuration should be a 24 inch outside diameter 
annular screen, 24 inches in height. 

Detailed design specifications and drawings were submitted to NASA/MSFC by 
MDAC and, following approval, competitive bids were sought i d a source 
selection made. The screen device, termed the Liquid Hydrogen Acquisition 
Device (LHAD), was subsequently fabricated by Western Filter Company; an 
assembly drawing is shown in Figure 2-1. 

Additional detailed drawings have been submitted to NASA/MSFC. 

2.1.2 Tank Design 

The tank is comprised of a 6061 T6 aluminum cylindrical pressure shell 
34.5 inches in length and 2$ inches in diameter, a 2-inch-thick, 29-inch 
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diameter plate, and additional hardware as shown In Figure 2-2. The LHAD Is 
bolted on the top of the plate. All tank penetrations for valves and ancillary 
plumbing are Incorporated In the tank bottom plate. The cylindrical cover Is 
equipped with a helical wall-mounted dip-brazed 0.25 inch outside diameter 
aluminum tube 70 feet In length. A summary of the design is presented below. 
Additional details are given In MDAC fabrication drawings submitted to NASA/MSFC. 

2.1.3 Cover 

The cylindrical pressure shell of the IDU is designed for 50 psid at ambient 
conditions. Crushing pressure is greater than 200 psid. Cryogenic allowables 
Increase the cryogenic operating pressures by approximately 40 percent over the 
ambient conditions. Indium tin seals are used on the lip of the cover:, these 
seals provide leakproof joints and are inexpensive and easily installed. 

There are no other seals, joints, or penetrations through the cover. The 
flange Is internal to the cover to maximize tank volume within the envelope of 
30 inch diameter. 

2.1.4 Plate 

The supporting plate contains ill tank penetrations. The main feedline valve 
and cone screen housing were mounted on the plate along the tank centerline. 

The refill valve is mounted flush with the bottom. The total length of the 
hardware outside of the cover and plate is minimized. Bulkhead fittings 
are used for the diffusers and pressure sensing line and Borders/Physical 
Sciences cryogenic feed-throughs are used for wiring and coaxial cables. 

2.1.5 Diffusers 

The diffuser connections are composed of aluminum bulkhead flared tube fittings 
with K-seals. The overhead diffuser is a lateral multijet formed by two rows 
of ten 1/4-inch holes. The tube diameter is 1 1/4 inch outside diameter to 
allow sufficient vent outflow rate for tank refill within 10 seconds. The 
submerged helium diffuser is a horizontal circular tube, termed a ring-tube 
diffuser, with one row of forty 1/8- inch holes in the bottom. Locating the 
holes in the bottom of the ring- tube minimizes the possibility of having 
liquid fill the inside of the tube and affect gas outflow. 


10 


j 


i 




11 





d«SMI 




















* > 


*■**¥,«* r . 




JL 










././ \ D/A 6HOLPS 


PEEP THRU 




O - n o 
Q \ 


GENERAL NOTES : 

LM^lOSS OT. ■*#*/,££ SPECIE teo: 

/ a4/jw iVfc? pep w/i-v#*//. Pt*er*AvT /*spfcr 
P£P ***+•!• <,Hh>T*P€1L PAU&6PAPH molds t66% Pep 
***** * J*5P’#JTS **JP££T #A9**&*A*tv To **«. -/* -v*,* Jl‘7»C 
2 PEST 90c'*e* TO 7TD PSH* tc **+ W*7» ' * 

SvMt'Hi/ o* CfApf A HPw+nw Am-Ptiopf cw «**•,*>« 1’Vt 

^ - «•*■ *x*b*i$ W9 sz£ A ***STi/*6 c*v tb*/T pP 

2+ t P^V J- SLAV'S /Wj) A **AAiVSU \0V‘VtH4TUe 
^'o*a^**d** c*«ns*T a* )**•**? *Siu#» r&t J pP o*u. 
f PPCPPPf r*PF f«3S Pfc S7AOf99~t 

■* CLf**/ Rex ?TP O1S9-01 **» ****** Pie H STPottp-ao* 
r <»DvtT**v BOAte p**. *u -9-Tees e*ce* r *au 06#ap*«. 
.Arret* re* ?tps jsv-* 

b ^S*< r&r M*e*FD JoetTS T* rxsor* ScCS. 

7 c SEALS P*Q*P JUS P/A A/DhjfA-Ymi vine 
w *€PiAP CMOS 

9 aacpshoSop AfpL* mr* "**pt < ue *«t* «ts *m 


*>* / 


WJ/AK - 9 0OLT 
**J bODEK */ 4 Si iRfc 
4fi A Rpfif 0 


I 


*>J 

‘ JU* ' ) 



£*} 4* Ones I 



r-L- 



f/t 



\ \ 

At->oK\ “ 

\ »y*> 

/ r ^ \ 
/ § 9 C \ 

t / 

” '*- 47^ I 

G y 

rrRn) 

r / 


c-z/j 




MTA.L -3 


.*<* 8*o*es 

f3T W<w ) 


g ^T ’^r T“ Lfy ACQUISITION 

^;53S3 ASSEMBLY- 

eomwmitmou*. IntMawnlwmn in 








2.1.6 Tee 

The bypass valve and main feedline valve are connected to the main outflow 
line through a tee connection. Several options were considered for the tee 
design since the 3-1nch tank-outlet line must be expanded to the 4 Inch outside 
diameter pump-inlet line. Use of a standard expansion joint (approximately 
7 degrees total Included angle) would Increase the total connection length. 

Use of a standard tee and expansion flange Increases the total pressure drop. 

As a compromise, a tee joint design which Incorporates the flow expansion Is 
used, as shown In MDAC Drawing 1T42474 (submitted to NASA/MSFC). 

2.1.7 Integrated System 

In the design of the tee, refill, overhead, and submerged helium diffusers, 
numerous flowrate and pressure drop calculations were made to assure 
practical operation. For hydrogen pressurization, a line diameter of at 
least 1/2 Inch Is required to meet the required volume flowrate of 

3 

1.7 ft /sec at 40°R with choked flow. In addition, the 3/8-1nch Valcor 
solenoid valve must be bypassed to achieve the maximum vent rate. Therefore, 
a 1-1/4- Inch-diameter bypass line Is provided. 

The refill diffuser has a pressure drop of 0.5 psld through 96 holes, each 

1/2 Inch in diameter. The corresponding Impact pressure of the liquid 

hydrogen flowing Into the tank and against the screen Is difficult to 

0 

determine analytically, but the inertial pressure ( P s P^-) Is of the 
order of 1 psld even If the flow from the diffuser Is assumed to expand 
to four times Its outflow area; this pressure exceeds the screen bubble- 
point pressure. Therefore, a shield Is used to stop possible bubble entrain- 
ment Into the screen device. 

The total pressure drop through the screen device Is 0.06 psl. The loss 
through the valve Is 0.17 psl at 7.5 Ib/sec. The flow loss through the 
tee Is approximately 0.5 psl. 

Estimates of the overhead pressurization effects on the LH 2 Interface 
were made, based on the results of visual pressurization tests reported In 


1 


13 




Reference 5. No significant cavity will be formed due to incoming flow 
except for the liquid surface very close (' 1 inch) to the diffuser. 

2.1.8 Instrumentation 

The complete Instrumentation list is shown in Table 2-1. 

2.1.9 Thermodynamic Vent System 

The thermodynamic vent system (TVS) was added to the IDU after the design of 
the IDU had been completed, and the fabrication and assembly operations 
selected. Therefore, the TVS design was compromised to meet the existing 
design constraints. An external wall-mounted, helical-coil tube design 
was selected to be continuously dip brazed on the cylindrical portion of the 
IDU wall. Tubing was not placed on the top of the IDU since heat flux through 
the top would be relatively low in ell cases and a separate brazing process 
would have added substantial costs. Furthermore, the wall-mounted TVS can 
be used with a slightly higher outflow rate and lower coolant temperature 
to remove the heat that enters the top plate of the IDU. 

The mounting arrangement of the TVS offers several alternative flow paths. 

Two, 1/4 inch diameter tubes, 35 ft long, are helically wound on the outside of 
the IDU, and connected by a tee manifold at the top of the cover. Coolant 
flow on the IDU wall can be routed as follows: 

A. A 70 ft continuous path is provided. Liquid enters at the bottom 
of the IDU, flows up through 35 ft of tubing to the top of the IDU, 
and then flows down through 35 ft of tubing to the bottom of the 
IDU. Coil separation is 3.5 inches on the IDU wall, with alternating 
uphill and downhill flow. 

B. Flow can start at the top of the IDU and flow down through one 
35 ft length of tubing (coil separation distance of 7 inches). 

C. Flow can start at the bottom of the IDU, and flow up through one 
35 ft length of tubing, exiting at the tee (coil separation 
distance of 7 inches). 

D. Flow can enter at the top of the IDU through the tee and flow 
down both 35 ft tubes in parallel paths (coil separation distance 
of 3.5 inches). 
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Level sensor/ Continuous level sensing 1 NASA-KSFC-GFP CL4 Attached to centra* 100 pipe (3 in. dia/(coax) 

capacitance (oscillograph, digital volt 

meter) 
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Temperature sensor/ IDU outer wall temperature 1 Thermal T19 Bonded to cylindrical IDU wall 70° to -424' 
platinum resistance Systems 4 wires - main tank 500 ohms at 
patch 5001-19 273.2° K 
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HF3-A surface temperature 1 International T?3 Bonded to HF3-A 

{needed for calibration) Thermal 2 wires - I DU 

Instrument 
Model -A 
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These options allow pressure drop through the tube to be varied by varying 
tubing length, and tube separation distances of 7 Inches or 3.5 Inches can be 
tested. The latter flow configuration (0) was selected for these tests since 
cooling capacity was maximum and the pressure drop was minimum. 

The flo.rate control required for the TVS was 0.3 to 3 Ib/hr of hydrogen. 

This was accomplished using five vlscojet.s orifices, selected to provide 
flowrates of approximately 0.2, 0.4, 0.8, 1.4, 1.8 lb/hr of hydrogen with 
approximately 15-psl pressure difference across the viscojets. Increasing 
or decreasing the pressure difference, and selecting the appropriate com- 
bination of viscojets allowed the flowrates to be controlled within a band 
exceeding that required (i.e, 0.2 to 5 lb/hr). 

2.2 ASSEMBLY 

The sequence of operations leading to the final assembly of the IDU is 
presented in this subsection. Figure 2-3 gives a side view of the Liquid 
Hydrogen Acquisition Device (LHAD) after its completion by the subcontractor. 
Western Filter Company, Inc., Chatswortn, California. 

A series of bubble-point tests were performed following fabrication of the 
LHAD. The film bubble-point test procedure was employed. The LHAD ;.as 
partially submerged in ethyl alcohol and then a spray of alcohol was used to 
wet the inner and outer screens. Pressure was applied to the LHAD as it was 
rotated while partially submerged. Bubbles released from imperfection points 
on the screen were noted and marked for later sealing. This type of film 
bubble-point test procedure is necessary for screen devices of large size 
relative to the head of test liquid which can be supported, and the procedure 
used for the LHAD can be easily extended for testing of vehicle screen devices 
in-situ. In the first test, several small holes were located and sealed with 
Allstate 430 Cadmium-free solder. Although this particular solder has been 
used successfully with filter elements at liquid helium temperatures, the 
possibility of crystallization exists. Therefore, a special polyurethane 
adhesive (Spec-lP20075 DPS-32330) which has been used successfully in a 
number of cryogenic applications was applied over the solder spots as a 
precaution. The device was then shipped to Garwood Laboratories for final 
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cleaning to NASA/MSFC Spec 164 and an acceptance bubble test, at which time 
approximately five small areas were sealed with the polyurethane adhesive to 
halt premature bubble breakthrough. The final bubble-point test was performed 
with isopropyl alcohol (68°F, o= 21.8 dynes/cm» p= 0.785 g/cc). The bubble 
point pressure was approximately 14 inches of water column, whereas an immediate 
and total loss of retention did not occur until 16 inches water column. The 
cone '-oreen device, used for feedline vapor control (see Test T), met a bubble- 
point pressure of 3.2 inches of water column, with isopron>l alcohol (68°F, 
o = 21.8 dynes/cm, p * 0.785 g/cc). 

Figure 2-4 shows the interior of the LHAD in an inverted position relative to 
its installation in the IDU. The protective layer of bolting cloth used on 
the outer layer of the LHAD is shown, as well as ports provided on both the 
top and bottom circular channels for inst'umentation. A 3/4-inch AN plug 
(screen fitting) drilled through and covered with 250 x 1370 SS screen, (not 
shown), is located to facilitate refill of the LHAD following breakdown. This 
screen fitting was bubble-point tested.* A teflon plug was used in the top 
of the control pipe to eliminate ullage space and minimize heat transfer 
from the ullage space to temperature probe T4. 

Figure 2-5 shows the tube routing during the early stages of assembly. The 
lines on the left are for pneumatic operation of the Parker Submerged LH^ Ball 
Valves. Figure 2-6 shows the initial installation of the IDU internal hardware 
prior to the placement of the LHAD. The overhead diffuser is removed before 
positioning of the LHAD, but is shown in place for completeness. The sub- 
merged helium pressurization diffuser has 40 equally-spaced 1/8-inch holes 
underneath. Placing the holes underneath assures that no liquid remains 
entrapped in the diffuser. The NASA GFP continuous capacitance-level sensor 
is shown. The cylindrical baffle guards against direct impingement of two- 
phase, high-velocity flow against the screen. Installation of the LHAD and 
much of the instrumentation is shown in Figure 2-7. The NASA GFP capacitance 
point-level sensors are shown; two are mounted on the central pipe; the top 
capacitance -level sensor is mounted on the pressure sens-ing line. Backup 
carbon point-level sensors are shown. Punch marks were placed on the LHAD, 
refill baffle, and plate to assist in proper alignment for subsequent reassembly. 


* Bubble point achieved was 22.4 inches of water column, using denatured 
ethyl alcohol (DPH 514) at 63.5°F, with a = 22.6 dynes/cm. 
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Figure 2-4 LIQUID HYDROGEN ACQUISITION DEVICE 
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Figure 2-8 shows the IDU before installation of the tank cover. Figure 2-9 
shows various operatio , involves in tank closure. The top left-hand photo- 
graph shows the IDU ready for installation of the cover. The top right-hand 
photograph shows the cover suspended by a crane. The test engineer is shown 
installing guide rods which position the cover as it is lowered onto the IDU 
plate. The lower left-hand photograph shows the lowering of the cover onto 
the plate. The lower right-hand photograph shows the installation of the Indium 
tin seal. 

Figure 2-10 shows the instrumentation wiring and heat flux gauges (4-inch by 
4-inch area) mounted on the inside of the cover. Immediately before the 
lowering of the cover, an electrical technician connects the instrumentation 
wires from the heat-flux gauges to the wires passing out through the electri- 
cal feedthrough mounted on the IDU plate. All other instrumentation wires 
were previously connected to the feedthroughs. The cover is then lowered 
and bolted. 

The IDU refill diffuser (1T42470) was tested on 2/1/73 to evaluate its 
effectiveness, prior to assembly. Water was used, with the setup shown in 
Figure 2-11. The test was run so that the degree of splashing could be 
observed over the range of fill depths. The maximum flowrate obtainable 

•5 

was 0.272 ft /sec. Figure 2-11 shows the diffuser installed upside down. 

Since the diffuser did not fill and allowed the water to flow out at 
approximately one half the height of the diffuser, the diffuser was only 
half as effective as it would have been in the proper orientation. This 
orientation appears to have maximized the degree of splashing, although the 
water flowrate was limited by the supply. Dynamic similarity conditions 
were only approximately achieved. 

The degree of splashing and disruption of the interface should be governed 
predominantly by the Froude number, Fr = V^/2g. With the IDU filled at the 
maximum LH^ inflow rate (0.5 ft /sec), the Froude number is 0.38, assuming 
the entire refill diffuser hole area (96 hole, 1/2-inch diameter) is 
effectively diffusing the flow. The Froude number for the water test was 

0.266. consequently, the LH2 inflow should be similar to that of the water 
test, even though the Froude numbers for the LHg and water are not equal. 
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A principal observation of the test was that when the liquid level reaches 
the diffuser, the degree of splashing and surface disruption Is greatly 
attenuated. When the diffuser was completely covered, there were only 
Insignificant ripples on the surface. Since the IDU normally operates with 
the refill diffuser submerged, it was concluded that the diffuser design Is 
satisfactory. 

2.2.2 Thermodynamic Vent System (TVS) Assembly 

The TVS coll Is a 70 ft length of 1/4-1 nch outside diameter tubing wound In 
parallel helical paths along the cylindrical tank wall, with a spacing between 
tubes of 3.5 inches. The tubing was attached to the wall by supporting clips, 
spot welded to the wall, and then the tube was dip brazed. The supporting 
clips were spaced approximately 10 inches apart along the length of the tube 
to secure the tube during the dip braze operation. However, thermal expansion 
caused the tube to separate from the wall in certain areas, and portions of 
the tube were not brazed to the wall. These areas, ranging in length from 
1 to 10 inches for a total length of approximately 10 feet out of 70 ft were 
filled with Wakefield Delta Bond which is a thermally conducting epoxy, and 
tnerefore, good thermal contact is maintained along the length of the tube. 

The TVS was originally configured so that liquid would enter at the bottom, 
flow upwards through 35 ft of tubing, and then return via the remaining 
35 ft to exit at the bottom. 

However, calculations of the two-phase pressure drops for the maximum expected 
flowrate (5 Ib/hr) indicated that for the 70 ft length path, the pressure drop 
could be as high as 12 to 15 psid. Moreover, there was some question regarding 
flow separation effects for the fluid flowing uphill through the 0.25 inch 
outside diameter TVS tube for the first 35 ft. By a slight modification of the 
routing, the inlet flow was brought directly up to the top of the IDU, where 
the flow splits and runs down in parallel paths. The running length is, there- 
fore, 35 ft, rather than 70, and the maximum flowrate in each tube is 
2. 1 b/ hr; the maximum pressure drop is, therefore, reduced by a factor of 

eight, since pressure drop is directly proportional to length and proportional 
to the square of the mass flowrate. 
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Portions of the TVS line leading up to the TVS heat exchanger coll are not 
brazed to the IDU wall. These lengths were Insulated with teflon tape to 
reduce extraneous heat transfer to the coolant fluid. However, this Insula- 
tion Is not a substitute for the foam Insulation necessary to maximize the 
operating efficiency of the IDU/TVS, but Is used to allow reasonable data 
to be obtained of the IDU/TVS operated under conditions of maximum heat flux. 

2.2.3 Instrumentation 

All Instrumentation wiring In the IDU was checked to verify continuity. The 
heat-flux gauges were tested with an Imposed hea\, flux to verify polarity and 
proper operation. The thermocouples were tested to verify operation. No 
problems were encountered. 

2.3 INSTALLATION 

The IDU is shown in Figure 2-12 immediately before installation in the 
260-gallon tank. Figure 2-13 shows the tank suspended from the BEMCO vacuum 
chamber lid. Figure 2-14 shows the tank, wrapped with high performance 
insulation, immediately before installation in the BEMCO. 

2.4 TEST HARDWARE LAYOUT, PLUMBING, AND INSTRUMENTATION 

The schematic drawing of the IDU/TVS layout In the MDAC test facility Is 
shown In Figure 2-15. The layout includes the IDU and associated plumbing 
and Instrumentation. A heat exchanger composed of 130 ft of 1/2-inch outside 
diameter tubing Is Included In the IDU pressurization line for use when cold 
(40°R) pressurant Is required. The TVS viscojet arrangement Is shown. Both 
GH ? and GHe pressurization can be supplied directly from storage tanks, and 
where necessary, a 1.5 ft pressure bottle can be used so that the approximate 
gas Inflow rate can be determined by measurement of the Initial and final 
pressure and temperature of the bottle. The mass flowrate measurement system, 
consisting of two hot-wire anemometer devices. Is Included. IDU/TVS and 
facility fill, drain, vent, evacuation, and purge components are shown. 

A schematic of the cold traps added to the system to prevent any contaminants 
in the helium from reaching the 2 inch ball valves- is shown in Figure 2-16. 
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Section 3 
ANALYSIS 


The principle analyses performed for this study are presented in this section. 

An analysis of the IDU wall-mounted thermodynamic vent system (TVS) was per- 
formed and a computer program written and documented. An analysis was per- 
formed of the transient pressure, temperature, and heat transfer response of 
the IDU, with internal hardware, to warm GHg pressurization. The computer 
program used (H431) was a modified version of that in use at MDAC for several 
years. The analysis and computer code used in sizing the IDU screen device are 
presented. 

3.1 ANALYSIS OF THE IDU WALL-MOUNTED THERMODYNAMIC VENT SYSTEM 

3.1.1 Objectives 

A TVS with a wall-mounted heat exchanger was incorporated in the IDU design 
to demonstrate the operation and control of this type of TVS. The specific 
objectives were to demonstrate (1) both steady-state and transient operation 
of the TVS, (2) control and predictability of the vented flowrate, (3) that the 
TVS could provide coolant flow for additional hardware such as feedlines and 
turbopumps, and (4) the fabrlcability of the wall-mounted heat exchanger 
subsystem. 

To meet these objectives, the TVS was configured with a number of flow-control 
orifices (Viscojets®) in parallel, which gave a vent flowrate variable from 
0.2 to 5 lb/hr. Bypass flow was also provided for turbopump or feedline 
cooling. The TVS heat exchanger coils on the IDU were configured to provide 
a number of alternate flow paths (as shown In Figure 2-15). Successful heat 
interception operation of the TVS required that the wall-mounted heat exchanger 
be insulated from the main tank fluid (to prevent condensation and loss of 
cooling capacity). However, funds were not available to provide the insulation, 
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and therefore, the performance tests of the TVS were made without inflation 
to determine flowrate as a function of applied pressure, and to verify flow- 
rate control. Heat transfer/cooling tests showed that insulation would be 
required. 

3.1.2 Operation 

The theory of operation of the TVS has been discussed in the literature 
(References 6 to 10) and is straightforward. The process is shown in 
Figure 3-1. The vented liquid is expanded isenthalpically to a lower pressure 
and temperature (A to B) and becomes a two-phase fluid at a temperature below 
that of the tanked fluid. This two-phase mixture flows through a heat 
exchanger at an approximately constant pressure where it boils by absorbing 
heat from the warmer fluid in the IDU, and/or by intercepting heat from out- 
side the IDU (B to C). Under steady-state conditions, a "design flowrate" will 
be just adequate to absorb all the heat entering the system, so that the vented 
fluid is completely boiled and exits the system as vapor. Additional heat 
transfer can increase the vapor temperature to equal or exceed the tanked 
liquid temperature, (C to D), which is equivalent to oriented vapor phase 
venting. This vented cold vapor can also be used to provide cooling of 
additional hardware, such as feedlines and turbopumps. 

In addition to performance characterization during steady-state operation, the 
transient performance of the TVS, while achieving steady state from system 
start, was evaluated. The thermal performance of the TVS depends on convective/ 
conductive heat transfer from the IDU fluid to the IDU tank wall, conduction 
along the tank wall to the heat exchanger tubing, and convection/boiling heat 
transfer inside the tube. Analysis of these heat transfer processes was 
necessary to evaluate the TVS performance. 

3.1.3 TVS Analysis 

The analysis of the wall-mounted TVS is presented in Appendix A for the case 
of a uniform heat-transfer coefficient between the heat-exchanger fluid and 
the inner TVS tube wall. Appendix B is an extension of this analysis which 
allows for the two heat-transfer coefficients encountered along the tube. In 
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the initial region, the heat-transfer mechanism is governed by annular flow 
with forced convection heat transfer and boiling in the liquid film and the 
Chen correlation is assumed for the heat transfer coefficients. At a critical 
quality, X^, the heat transfer mechanism is transformed to forced convection 
of a two-phase "mist" flow. This second region is termed the liquid deficient 
regime. 

The Dittus-Boelter heat transfer coefficient is used in this regime. Details 
of the heat transfer coefficients are given in Appendix C. 

The principal unknown in the TVS analysis is the fluid quality at which 
transition from annular flow to liquid deficient flow occurs. Various 
reference data indicated that x c r -jt< ca -| the order 0 f 0.85 for hydrogen 

in 0.25-inch tubes and, therefore, this value is assumed. The IDU/TVS program 
was used to generate toe required coolant flowrate for an uninsulated tank as 
a function of x c r ^ ca ]» all other parameters being identical, and the results 
shown in Figure 3-2. It is seen that the total flowrate varies from 5.4 to 
6.2 Ib/hr for x critical var Y in 9 from 0.7 to 0.9, and therefore, the variation of 



CRITICAL QUALITY (TRANSITION TO LIQUID - DEFICIENT FLOW) 
FIGURE 3-2. COOLANT FLOW RATE DEPENDENCE ON CRITICAL QUALITY, X 
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Critical not ^ ave a 9 reat e ^ ect on the coolant flowrate. However, 
there Is a significant variation In the lengths of the annular and liquid 
deficient region. For example, at Xg 8 0.65, the annular region Is 7.72 ft 
(out of a total tube length of 35 ft) whereas, for Xg 8 0.85, the annular 
region Is 16.16 ft out of 35 ft. As a result, when the IDU Is partially 
filled, the decreased heat transfer In the gaseous ullage region, coupled with 
a low critical quality, will lower the effective overall cooling capability of 
the TVS. This can occur In Its present configuration, since the annular flow 
(high Internal heat transfer coefficient) will occur In the upper part of the 
tank where the wall heat-transfer coefficient Is decreased due to vapor. In 
addition, heat entering through the bottom portion of the tank, which is 
convected to the liquid Interface, would not be removed efficiently if the 
TVS In that region has liquid deficient flow. 

Two test conditions are analyzed using the IDU/TVS program. The first case 
corresponds to the IDU surrounded by GH£ at a pressure greater than the 
saturation pressure of the IDU LHg. The IDU is not insulated, and therefore, 
condensation occurs. The required coolant flowrate for two IDU LHg saturation 
pressures is plotted in Figure 3-3. Details of the results are given in 
Table 3-1, The second case is for free convection heat transfer from warm 
GHg, surrounding the outer IDU wall. The coolant flowrate required to maintain 
the IDU liquid at a constant temperature is shown in Figure 3-4, and tabulated 
results are given In Table 3-2. Insulation conditions are shown for comparison. 

Although the flowrates required of the TVS are of the order of 4 lb/hr, the 
system is capable of flowrates of the order of 10 lb/hr with a pressure drop 
of approximately 2 psi in the TVS line. However, the corresponding pressure 
drop through the viscojets would be approximately 70 psi. Details of the 
pressure drop versus mass flowrate are given in Appendix D. 

Appendix E describes the computer code, and the listing is given in Appendix F. 

Complete results generated by the IDU/TVS computer program are given in 
Appendix G for the IDU insulated with 0. 25-inch and 0.5-inch of Saturn S-IVB 
type foam. These results are included so that comparisons can be made between 
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I DU/TVS COOLANT FLOWRATE (LBS/HR) 



Figure 3-3 


I DU/TVS COOLANT REQUIREMENTS - CONDENSATION 
ON 1DU EXTERIOR - NO INSULATION 



Table ?-l 

PREDICTED OPERATIONAL REQUIREMENTS FOR IDU/TVS - CONDENSATION ON IDU EXTERIOR, NO INSULATION 
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the operating characteristics of the uninsulated and insulated IDU/7V3. The 
results of Appendix G clearly show the need for Insulation, condensation heat 
transfer Is eliminated, and the net thermal resistance Is increased. Use of 
1/4-in. foam thickness reouces the coolant flowrate to approximately one-fifth 
the values shown In figure 3-4, and use of 1 /2-in . foam thickness reduces the 
coolant flowrate to one-tenth the values of Figure 3-4. The IDU/TVS is capable 
of providing flowrates of the order of 5 lb/hr, whereas only approximately 
1 Ib/hr coolant flowrate is required for the IDU insulated with 1/2-in. of 
foam, and surrounded by 100°R gaseous hydrogen. 

3.2 IDU AUTOGENC-5 PRESSURIZATION 

A versatile tank pressurization analysis computer code has been in regular 
use at MDAC for several years. This modular code was designed to facilitate 
modification and adaptation to specialized studies of heat and mass transfer 
which determine the course of propellant tank thermodynamics. Computer 
output includes the time-variable histories of the temperature distributions, 
heat and mass transfer rates, tank pressure, and pressurant flowrate. 

3.2,1 Basic Tank Pressurization Computer Code 

The analysis of tank pressurization performance is based on a one-dimensional 
model. Spatial variations in the system variables can occur or.ly dlong the 
vertical tank axis; there are ro radial or circumferential variations. Tank 
pressurization computer programs based on this type of model have been compared 
extensively with experimental data and found to be valid. Buoyancy tends to 
produce a stable thermal stratification in the gas and liquid, giving a 
temperature distribution which is essentially one-dimensional. 

The thermal system for this analysis consists of the tank wall, internal 
hardware, propellant liquid, and ullage gas. Any size and configuration may 
be specified for the tankage. The ullage gas may be pure propellant vapor or 
a mixture with helium. Tabulated variable properties are used to describe 
the thermodynamic behavior of all materials. 

The nomenclature and mathematical model for the code are delineated in Table 3-3. 
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Table 3-3. TANK PRESSURIZATION COMPUTER CODE NOMENCLATURE 
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The computations are based on a finite-difference representation of the 
physical system. The tank wall, internal hardware, propellant, and ullage are 
each divided by horizontal planes into a number of nodes, the properties 
within each node being uniform. The gas and liquid are divided into nodes whose 
thickness and location can vary with time. The tank wall and hardware nodes 
are of equal axial thickness and are stationary. The size and number of these 
nodes is sufficient to give an adequate step function approximation to the 
continuous axial variation of the system variables. Gas and liquid noc'es may 
be subdivided or combined as required to meet solution accuracy criteria. 

The volume of each liquid and gas node is bounded by the top and bottom 
boundary planes and by the solid surface of the tank wall and the internal 
hardware. Heat transfer takes place between each gas node and the solid 
surfaces with which it is in contact. The physically simultaneous processes 
of heat transfer and pressure change are assumed to take place sequentially 
as isobaric heat transfer and isentropic pressure change. The numerical 
solution is obtained by calculating the change in the state of each node in 
the system during each successive time step throughout the total solution 
time. The state of each node is determined from equilibrium, conservation 
relationships. 

Duty-cycle control data is input to the program in time-variable tables. These 
include propellant outflow rate, vehicle acceleration, pressurant inlet 
temperature and composition, tank pressure, and pressurant inflow rate. Heat- 
transfer coefficients may be input or calculated internally from free con- 
vection or other . elationships. The solution may be computed in two modes: 

(1) either the tank pressure schedule is specified and the required pressurant 
flowrate is calculated or (2) the pressurant flowrate is specified a;id the 
pressure calculated. Initial conditions specified at the start of the duty 
cycle are the ullage fraction, tank pressure, and all temperature 
distributions. 
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3.2.2 Computer Code with Conduction 

This basic computer code Is available In various specialized versions with 
added capabilities. One such version analyzes the transient heat-transfer 
processes throughout the mission. Thermal conduction equations were added 
for one-dimensional heat transfer axially down the tank wall, the liquid 
and the ullage gas. The propellant remains settled. Heat is transferred 
between the gas and wall by free convection. Heat conducted down the wall to 
the liquid surface leve* is transferred directly to the surface liquid node. 

The excess heat conducted from the gas to the liquid surface over that con- 
ducted from the liquid surface to the bulk liquid results in liquid vaporiza- 
tion. Conversely, a deficiency in this exchange will result in vapor 
condensation. Heat added to the liquid surface node by the wall conduction 
raises the liquid temperature. If the surface liquid node reaches saturation 
temperature, further heat addition causes vaporization. 

With these assumptions, the tankage system is virtually always in a transient 
nonequilibrium state when prepressurization and propellant outflow events 
occur. The duty cycle can include continuous pressure control, giving 
pressurant requirements to maintain tank pressure, or pressure decay, giving 
the tank pressure history following the shutdown of the pressurization system. 

3.2.3 Computer Code Modifications for IDU Analysis 

The specialized version of the pressurization computer code used for the IDU 
analysis was further modified for this application. First, the code was 
simplified by removing components that were of no use to this study; this 
included the diffusion equations, applicable to a two-component ullage, and 
the associated interface heat and mass transfer model for use with helium. 

Then the available methods of numerical solution were reviewed to select the 
most efficient for this study. Both Gauss-Jordan and Gauss-Seidel methods had 
been adopted to this code. The simple Gauss-Seidel method was found to give 
identical results to the more reliable, but more expensive (in core storage 
and computing time) Gauss-Jordan method for typical test cases, so the former 
was adopted for this application. 
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The solid conduction routines were called to operate on the internal hardware 
as well as the tank wall. The various parts of the internal hardware (screen 
device) were lumped into a single component of equivalent mass and area. 

Different items of internal hardware could be treated individually with 
relatively simple modifications to the code, but this appeared to offer little 
advantage and would have greatly increased core storage and computing time. 

The area and mass of the flat tank-top is accurately represented in the heat 
transfer calculations. However, the top must be treated as having a uniform 
temperature. Conduction between this mass and the cylindrical wall occurs 
according to the one-dimensional model with no provision for radial conduction 
analysis in the tank top. 

Since the conditions of the I DU operation would permit portions of the wall 
and hardware to be at temperatures less than the saturation temperature of 
the vapor, it was necessary to add the computations for condensation heat 
transfer. Since condensation heat transfer rates are much greater than free 
convection rates, this led to an added complication. With the low specific 
heat of the metal at liquid hydrogen temperature, condensation heat transfer 
would give a rapid temperature rise during a single time step. The final 
temperature could greatly exceed the saturation temperature making the use 
of the condensation heat-transfer coefficient invalid for part of the computation. 

To correct this inaccuracy, the final temperature of the wall or hardware node 
is checked at each time step when condensation ' r cjrs. If it exceeds the 
saturation temperature, the computation is repe in two steps: first, the 

time required to reach the saturation temperature in the metal node is 
determined at the condensation heat transfer rate; then, the free convection 
heat-transfer coefficient is calculated and the metal node temperature increase, 
starting from the saturation value, is calculated for the remaining time of 
the time-step interval. This two-part computation removes the possibility of 
a significant error in the use of condensation heat-transfer rates. 
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When condensation occurs at any point on the wall or hardware, the condensed 
mass Is removed from the gas at the point of condensation and Is added to the 
surface liquid node. The manner In which the condensed liquid travels this 
distance Is not treated. 

When the screen device does not contain liquid, It Is treated as Internal 
hardware In the conventional manner; that Is, Its temperature distribution 
will change as heat and Is transferred to It from the warmer gas. To simulate 
the liquid-filled screen device, It Is treated as a constant-temperature heat 
sink. This "constant" temperature may be either the initial temperature 
distribution or a uniform saturation temperature. In the former case, 
condensation will occur at any point that is below the saturation temperature; 
in the latter case, condensation will not occur, but the temperature itself 
will change to follow any tank pressure changes. 

3.2.4 IDU Autogenous Pressurization Results 

The MDAC H431 tank pressurization computer code, modified to model the IDU, 
has been used to generate results for two warm GHg tank prepressurization 
cases. In Figure 3-5, the GH., pressurant inflow rate and ullage temperature 
at the top of the tank are plotted versus time. A 5-second prepressurization 
rise from 15 to 40 psia is assumed. The temperature profile results of 
Figure 3-6 show that the initial ullage is compressed by the incoming 
pressurant and undergoes a temperature increase of 10°P, After approximately 
400 seconds, the tank temperature gradient varies slowly, and the incoming 
pressurant mass flowrate is less than 10 lb/sec. 

Results shown in Figures 3-5 and 3-6 are for an adiabatic tank corresponding 
to the IDU with insulation. Figures 3-7 and 3-8 correspond to the conditions 
for a noninsulated tank with free-convection heat transfer from the outer tank 
wall to 37°R Gl^, (Heat transfer to the LH^ through the tank wall from the 
surrounding GHg is neglected, since it has little effect on ullage conditions.) 
The characteristic time for the ullage temperature profile and incoming gas to 
reach a steady-state condition is somewhat longer than for the adiabatic tank 
wall case. 
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Figure 3-5 TRANSIENT PRESSURIZATION CONDITIONS IN IDU 
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Flqure 3-6 PREDICTED TEMPERATURE DISTRIBUTION FOR IDU PREPRESSURIZATION 
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Flour#* 3-A PREDICTED TEMPERATURE DISTRIBUTION FOR IDU PREPRES3URIZATI0N 
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3.3 SCREEN SIZING COMPUTER PROGRAM 

The expulsion characteristics of an annular screen acquisition device are 
predicted by the computer program in Appendix F. The total AP across the 
screen in steady-state operation is the sum of several contributing factors 
evaluated by the program: (1) flow loss accompanying the movement of liquid 

through the screen into the annular outflow passage, (2) hydrostatic head loss 
as the liquid moves upward away from the bulk liquid pool, (3) viscous flow 
losses within the annular passage, and (4) dynamic pressure loss (1/2 pv ) due 
to fluid motion at the most critical point in the passage. Liquid expulsion 
can continue only as long as these flow losses do not depress the pressure 
within the outflow passage to the point where pressurant gas can break into 
the passage. 

One of two experimental correlations can be used to determine the pressure 
loss for flow through the screen. If the screen mesh forming the walls of 
the passage is 200 x 1400 mesh or coarser, flow correlations established for 
cryogenic fluids in Reference 11 are used: 

A,L kJL ~ 

AP =-L- pV + ~4r p\r (3-1) 

D* 


where 


V = average flow velocity through screen 
D = average screen pore diameter 
L * screen thickness 
A^A^ = empirical constants 

If a finer mesh screen is present (or if a comparison is made between fine and 
coarse mesh screen performance), the correlation developed in Reference 12 is 
used. Ambient temperature fluids were used to establish this correlations. 

iP . LIm! (c-^+c,) (3- 

e^O 1 v c 
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where 

e s screen volume void fraction 
a » surface area to unit volume ratio 
Lfa s empirical constants 

Backup material to this screen (such as perforated plate) can be Introduced 

! Into the mathematical model described by the computer program. 

! 

5 

; The computer program requires as Input the acquisition device dimensions, 

outflow rate, and screen characteristics. It then computes the most 
critical pressure differential across the screen liquid-gas Interface as 
a function of the liquid level within the tank containing the acquisition 
device. At each step this critical AP Is compared with the bubble-point 
pressure of the screen to determine whether breakdown has occurred. 
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Section 4 
TESTING 

The objectives of this phase of the program were to: (1) perform a component 
and overall system checkout of the IDU/TVS and (2) explore certain fluid 
dynamic and thermodynamic phenomena associated with screen device propellant 
acquisition. The basic checkout tests Included component functional tests, 
tank-leakage and proof-pressure tests, determinations of the LHAD screen 
retention breakdown ("bubble point") condition as a function of flowrate, 

TVS flowrate control, and a thorough demonstration of the operational 
characteristics of the "start tank" principles. Including outflow and 
refill with simultaneous outflow. 

Additional tests were perfo»med to evaluate the IDU/TVS thermal control, 
transient thermal response, effects of warm hydrogen and helium pressuri- 
zation on breakdown, effectiveness of the cone screen for feedllne vapor 
control, two-phase refill, and hydraulic pressure transients due to valve 
actuatlon/deactuatlon. 

As detailed In Test G and Appendix J, problems were encountered with the 
three Parker LH^ ball valves which Impaired the capability to meet the test 
goals. Extensive effort was made to correct this situation, and a solution 
Is proposed which Involves relatively minor modification to the valve 
pressurization system. The valve problem caused severe test schedule 
problems, directly affected test operation, and eliminated two tests 
(Tests R and W). However, useful data were obtained In the major areas 
investigated, and the successful operation of the IDU/TVS was demonstrated 
in all other respects. 
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Detailed, annotated data are presented in a separate document, the Supplemental 
Data Document (SDD), which has been submitted to NA5A/MSFC. This section 
discusses the test results and outlines the data available. More detailed 
evaluation of Individual tests is provided in the SDD. 

4.1 TEST A - SOLENOID ACTUATION 

The solenoid valves were actuated under ambient conditions, and current and 
voltage recorded. All valves drew 0.8 to 0.85 amperes at actuation, with 
28 VDC applied. 

4.2 TEST B - BALL VALVE ACTUATION - AMBIENT 

The three Parker submerged LH^ ball valves were actuated five times to verify 
their operation, using 500 psig GHe. The output of the position indicators 
was verified using ohm meters. The valves deactuated properly without back 
pressure. On 11 December 1973, the minimum pressure to actuate the valves 
was determined for ambient conditions. No line pressure was applied, and 
therefore, the initial deactuation pressure was approximately the same as 
the actuation pressure supplied. The minimum actuation pressures were 
150 psi for PV1 and 175 psi for PV2 and PV3. 

4.3 TEST C - PROOF PRESSURE TEST - AMBIENT 

The IDU tank was proof tested at 75 psig with GHe for five minutes and 
inspected. There were no indications of deformation. 

4.4 TEST C - LEAK TEST - AMBIENT 

The I.’U was pressurized with GHe to 40 psig and all joints sniffed with a 
mass spectrometer. Leaks at the tee connections used for the platinum 
temperature transducers were corrected by tightening the fittinqs until the 
mass spectrometer revealed no leakage. The sensitivity of the mass spectro- 
meter was shown to be 5.8 x 10”® sees. 

_5 

Two leaks were detected. A very low leak (10 sees) was found at the K-Seal 
on the bottom of the plate, at the overhead diffuser fitting. A second leak 
was detected at the cone screen housing which bolts to the plate arid forms 
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a portion of the main IDU feedllne. Upon further testing, It was found that 
the -7 port of the housing had a small crack In the weldment. This was 
rewelded, the Indium tin seais changed, and the housing reassembled. Sub- 
sequent leak tests revealed no Indication of leakage anywhere In the system, 
except for the K-Seal. This low leakage rate Is small compared to the nominal, 
unavoidable leakage through the Parker valves, and therefore, no attempt was 
made to correct this leak. 

4.5 TEST E - CONDITION LEVEL SENSORS (PRIOR TO ASSEMBLY) 

The carbon-resistor, point-level sensors were conditioned by dunking In LH 2 
and examined for damage. No damage was detects 

4.6 TEST F - SOLENOID ACTUATION - CRYOGENIC 

All solenoid valves were actuated at cryogenic temperatures during the test 
series. A variable vcltage power supply was used to control the current to 
the solenoid valvc-s to assure that the current to each solenoid was approxi- 
mately 1 ampere. At ambient conditions, the solenoid valves (Valcor SV 90) 
use 28 volts at 0.8 ampere. At reduced temperatures, the decrease in 
resistance allows a larger current to be supplied, thus increasing the 
magnetic force, w..le maintaining the same rate of electrical energy dissl- 
patlon, I R. At -420°F, current of 0.8 ampere corresponded to a supply voltage 
of 2 volts. In practice, 1 to 2 amperes at approximately 3 to 5 volts was 
normally supplied to the pairs of solenoid pilot valves to actuate the Parker 
submerged liquid hydrogen valves. At pressures above 500 psl, and at tempera- 
ture of the order of -420°F, the solenoid pilot valves opened sluggishly, and 
at 750 psl would often fail to open. 

4.7 TEST G - BALL VALVE ACTUATION 

On 13 November 1973, tests were Initiated and the 260-gallon main tank was 
filled. Following a wait period of several h^urs, the Parker submerged LH^ 
ball valves (PV1, PV2, and PV3) were actuated, but did not respond properly. 
Numerous attempts to actuate these valves were made, until the vent check 
valve (V43) froze shut due to cold helium being expelled from PV1 , PV2, and 
PV3. With V43 frozen shut, the Parker valves could not be vented. Check 
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valve V43 was freed by blowing air over it, and additional attempts were made 
to actuate the valve. The valve signals ("talk bacK") indicated the valves 
could not be closed completely. During a main-tank fill operation, liquid 
entered the IDU, proving that either PV1 and PV3, and/or PV2 were at least 
partially open. Anticipating that the problem Uy with the pilot valves 
(SV11, 12, 13, 14, 15, and 16) not opening properly, a solenoid valve was 
installed upstream of V43 (see Figure 2-15), this valve was used to apply 
pressure to PV1, PV2, and PV3 simultaneously. Again, the PV valves could not 
be made to operate properly; actuation pressures were increased to 750 and 
then to 1000 psi, but to no avail. On 14 November, the PV valves again did 
not actua e, and therefore the tests involving the PV valves were postponed 
and the thermodynamic vent system (TVS) tests were initiated. On 15 November, 
the system was shutdown and, when the temperature reached -110°F on 

19 November, the valves were actuated properly. 

During consultation with George Young and Haschal Hyde of NASA/MSFC on 

20 November, it was c. ermined that extremely pure helium was required to 
actuate the valves, and it was agreed that samples of helium would be tested. 

Proper operation at -110°F was an indication that ice was probably not present 
in the valves; although thermal expansion could have caused affected parts to 
break loose of ice, this was not considered likely. Subseque t mass s{. ro- 
rraph analysis of the helium disclosed no measurable amounts of water vapor. 

The failure of the PV valves to operate when the pilot valves were not actuated 
but the external solenoid valve upstream of V43 was used, indicated that the 
problem lay with the PV valves. Although it was possible that the pilot valves 
had failed, subsequent testing confirmed that the solenoid pilot valves (SV11, 
12, 13, 14. 15, and 16) were operating properly with 500 psi helium. 

It was agreed that a long-duration vacuum purge coupled with numerous valve 
operations with clean helium would be performed to remove any contaminants 
in the valves and lines. 
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A helium sample was taken near the PV valve pressurant line inlet to the vacuum 
chamber lid on 20 November, and on 21 November the IDU system was pumped down 
for the entire day. The helium sample results were received on 26 November; 
the mass spectrograph analysis indicated the presence of 0.7-percent and 
0.2-percent O^, which indicated that the helium supply system was contaminated 
with air. 

A second sample was taken within the helium supply system, which Indicated 
0.3-percent N 2 and 0.03-percent 0 2< Two additional samples were taken of the 
helium in the high-pressure bottles of the supply trailer. These samples 
indicated (a) 0.?4-percent N 2 , 0-percent 0 2 , and (b) 0.1-percent N 2 , 0.02-percent 
Og. Since the samples from the helium trailer indicated contaminant levels 
within the experimental accuracy of the mass < .'sctrograph analysis and 
sampling technique, it was assumed that the helium in the trailer was 
sufficiently pure. A helium pressurization line was run directly from the 
helium to the Parker valve pressurization system. 

On 29 November 1973, it was noted that a leak somewhere in the 500-psi PV 
valve pressurization system was causing a loss of 3.15 SCFM of helium when 
the lines were pressurized to 500 psi. It was decided to continue testing, 
but to minimize the helium loss wherever practical by decreasing the line 
pressure after actuating the PV valves. It was later determined that a 
silver soldered 1/4-inch outside diameter joint at the top of the 260-gallon 
tank was the leak site. 

Tests were initiated on 5 December, following a full day of vacuum purging 
and valve operations. At -125°F, during the tank chilldown and loading, the 
PV valves were actuated but responded sluggishly, which may have been due to 
differential of thermal expansion problems, caused by n t having reached a 
steady-state temperature in the valves. Both the external and internal pilot 
valves were used, but the PV valve response was the same. After approximately 
one hour at -420°F, the valves were again actuated, but would only open part 
way, followed by erroenous "closed" signals with pressure applied. However, 
following another hour of thermal adjustment, the valves were actuated with 
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both the internal (pilot) solenoid valves and the external valve. Six 
successful PV actuations were recorded. Following a 45-minute period, test- 
ing was initiated on the IDU fill and drain, but the valves began to behave 
erratically. One slow expulsion was obtained through the refill valve (PV2) 
since the main feedline valve, PV1, would not open. Following this, none of 
the valves could be made to operate, even at 750 psi. However, test data were 
obtained for calibration of the level sensors, even though PV1 was frozen 
closed and PV2 a-.d PV3 were frozen open. Testing was halted on 5 December. 

By 7 December, the IDU temperature had reached -250° F. The PV valves were 
actuated several times and the signals indicated proper operation. However, 
when the IDU was vented or pressurized, the main-tank pressure remained equal 
to the IDU pressure, which either indicated a severe leak or the PV valves 
were not completely closed. Subsequent testing proved that there were no 
leak paths between the IDU and main tank, other than the PV valves. On 
10 December, the IDU temperature was above 32°F. All valve signals indicated 
that the valves were closed. Both the IDU and main tanks had 10 psig. The 
IDU was vented, and again the main tank pressure followed along. The valves 
were actuated several times, and both tanks pressurized to 10 psig. The 
IDU was then vented, but the main tank pressure did not decrease. Therefore, 
the microswitches were apparently not properly adjusted to indicate complete 
closure of the PV valves, or are inherently unreliable at liquid hydrogen 
temperatures.* 

A liquid hydrogen cold trap was then installed to freeze out any trace 
impurities in the helium supply, and the separate line connected directly 
with the helium supply trailer used as the source. The iDU was vacuum 
purged for a full day, with intermittent valve actuation (at least 10 valve 
actuations for each valve), for a total of five complete vacuum purges. 


*C.E, Schroeder, MDAC Saturn Propulsion, has investigated microswitch per- 
formance for valve position indications at liquid hydrogen temoeratures, 
and has found them to be unreliable. Microswitches have indicated valves 
open, when they were closed, and vice-versa, and have failed to close when 
the va'Wes were partially open or completely open or closed, giving an open 
circuit indication. Similar observations were made of the Parker valves. 
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On the next test day, 12 December, the cold trap was first filled with LHgi 
helium was then passed through the heat exchanger to freeze out any foreign 
gases, and the Parker valves actuated 10 times with the purified helium. 

This operation served to clean out the Parker valves to an additional level, 
before filling the 260-gallon tank with LH 2 . 

During the main-tank fill, when T6 indicated -387°F, PV3 was closed per the 
test procedure. This was required to ensure that the decreased pressure in 
the valve pressurization line below the check valve, V43, would be minimized, 
and therefore, decrease the probability of air flow through the check valve, 

| V43, in the event it was leaking. PV3 was then commanded open with 200 psi 

! line pressure, but did not actuate. The line pressure was increased to 

1 500 psi and the valve opened upon command. All subsequent operations to 

diminish the line pressure decrease were performed by opening and closing 
Vll , to pressurize the line up to V62 with 500-psi helium. V62 was then 
opened so that the helium trapped upstream would flow into the cold line, but 
would not cause the Parker valves to open. This operation assured that no air 
could leak past V43 into the helium pressurization line to the Parker valves. 

During the IDU leak test of 12 December, valve V63 was closed because the leak 
in the high pressure helium line (inside the main tank) to the Parker 2-inch ball 
valves was causing an excessive lo f ^ of helium. PV3, which had been maintained 
open, did not close completely, but could be completely opened. It was left 
in the partially closed condition with V63 closed to conserve helium. 

The erratic behavior of the PV valves outlined above continued throughout the 
tests. Significant valve failures are discussed in conjunction with each test 
in the following sections. Numerous attempts were made to induce the valves 
to operate properly, including increasing the supply pressure to 750 psi, 
decreasing it to 250 psi, increasing the current to the solenoid pilot valves 
to cause more rapid opening and closing, and frequent valve cycling as well 
as valve cycling with long periods in between. None of these attempts were 
successful. As discussed in the following sections, the valves would sometimes 
fail to open, or would fail to open all the way, and at no point did PV2 and 
PV1 or PV3 close completely, thus pressure isolating the IDU from the main tank. 
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The reasons for the problems encountered with the Parker valves have not been 
fully ascertained. A number of factors have been identified which are 
potential sources or contributory sources of the erratic valve behavior 
observed, as listed below. 

o Parker baked the valves at 165°F for two hours to remove moisture 
before their qualification tests. This was not done with the valves 
before assembly into the I DU. 

• Tube routing was performed in a shop area, rather than in a clean 
room. Pressurant tubing was cleaned before assembly but assembly 
progressed in several stages, with modifications which required the 
removal and replacement of tubing. Contaminants, therefore, could 
have been introduced into the high-pressure helium pressurization lines. 

• The facility helium supply was evidently contaminated and appeared 
to be a major cause of the initial problems encountered. However, 
subsequent purging and use of helium directly from the trailers with 
a LH^ cold trap did not solve the problem. 

• Each Parker valve was pressurized and vented by a pair of two-way 
solenoid valves, since a LH 2 submersible three-way solenoid valve was 
not available. The valve closure required that residua 1 pressure and 
spring force close the valve, with both of these forces diminishing 

as the valve approached the closed position. In principle, the valves 
can be operated to use supply pressure (500 to 750 psi) for both open- 
ing and closing (dual pressurization and venting) but this would have 
required the use of additional solenoid valves and cables. An analysis 
showing the comparative operational factors of safety is given in 
Appendix J, which shows that dual pressurization and venting significantly 
increases the net actuation/deactuation forces. 

• The solenoid pilot valves selected were rated for 1,000 psi but would 
sometimes fail to open or would open slowly at pressures of 750 psi. 

In principle, the use of higher actuating pressures should have helped 
in the 2-in. ball valve operation. However, we found no improvement 
throughout the test program as a result of the frequent use of 750-psi 
actuation pressure, rather than 500-psi pressure. 
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• All three valves were observed not to operate after a 1-hour thermal soak 
period at -420° F, but in the one case appeared to actuate properly for 

a total of six times following a 2-hour soak period. (However, the 
position indicators may have been giving erroneous readings.) Thermal 
transients and differential coefficients of expansion could, therefore, 
have been a contributing factor to the observed failures. During the 
testing, it was necessary to sometimes drop the main-tank liquid level 
beneath the valves, exposing them to the main-tank ullage gas. However, 
the temperature differences were small (usually of the order of 10°F to 
50°F), and the effect would be expected to be minor. 

• The piston seal effectiveness under cryogenic conditions is not known. 
Defective seals would cause more rapid pressure equalization across 
the piston, which decreases the net actuation force. 

4.7.1 Conclusions 

The valve behavior may be due to: (1) peculiarities of the valves (i.e., history 

of usage or de. ign, (2) contamination, or (3) use of contaminated helium. How- 
ever, the analysis of operational factors of safety, given in Appendix J, 
indicates that modifying the valve pressurization system to pressurize the 
actuator during opening and closing should increase the net actuation forces 
by at least a factor of two; this modification appears to present the greatest 
probability of solving the valve malfunction problem, and is therefore, the 
recommended "least cost" approach, 

4.8 TEST H - LEVEL SENSOR EVALUATION 

A number of tests were performed to compare the response of the carbon point- 
level sensors, capacitance point-level sensors, and the capacitance continuous- 
level sensor. A calibrated curve of the liquid-level versus continuous-level 
sensor signal voltage and corresponding level indications from the carbon 
point-level sensors is given in Figure 4-1. The capacitance continuous-level 
output (CL4) is linear above the 6-inch level and the carbon point-level 
sensors, RL1, RL2, RL3, RL4, RL5, and capacitance point- level sensors, CL1 , 

CL2, CL3, accurately indicate liquid lpvel, as shown by comparison with CL4. 

Below the 6-inch level, there is an apparent nonlinearity. In addition, the 
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voltage output and capacitance are nonlinear, as shown In Figure 4-2, but 
since the level range Is In the linear portion, this effect should not have 
caused the nonlinear voltage versus height portion of the CL4 curve. Since 
the voltage reading with no liquid In the tank Is 1.38 volts, and the voltage 
reading when static breakdown occurs Is 1.4 volts, at a level 2. 25-Inches 
above the tank bottom, there Is an unexplained nonlinearity In this region. 
Unfortunately, there are no level sensors below 5.5 Inches, and therefore, no 
Independent measurements can be made. 

However, the calibration curve In this region can be inferred. Since liquid 
dropout occurs at breakdown, the liquid held in the annular screen region of 
the LHAD (above the free liquid surface), flows out, increasing the height of 
liquid measured by the probe a known amount. In most cases with liquid 
hydrogen saturated at 15 psia, the measured voltage at breakdown increases 
from 1.5 to 1. 9-2.0 volts, which corresponds to the dropout of liquid in the 
LHAD. The volume of liquid depends on the degree to which the uppermost 
region of the LHAD is full, and this is indeterminent. However, the approximate 
amount of liquid corresponds to that in the annulus, above a height of 2.25 to 

2.5 inches from the bottom of the I DU, where breakdown occurred in most cases. 

3 

The annulus volume is 1,26 ft and the height increase after dropout is 

3.6 inches. Therefore, the cases for which breakdown as measured by CL4 
occurs at 1.5 volts, followed by an increase to 1.9 to 2.0 volts due to liquid 
dropout, correspond to breakdown at 2.25 inches (the bottom of the screen) 
followed by the dropout height increase to 5.85 to 6.1 inches. The calibration 
curve thus exhibits the nonlinearity shown in Figure 4-1. 

The carbon point-level sensors responded rapidly to liquid level during fill, 
indicating a liquid level which corresponded closely to both the capacitance 
point-level sensors and the capacitance continuous-level sensor. However, 
during expulsion, the carbon point-level sensors sometimes gave late indica- 
tions of exposure to ullage gas, probably due to a film of droplets of liquid 
hydrogen being retained on the resistor element. Examples of the response of 
RL1 , RL2, RL3, RL4, and RL5 are shown in Figure 4-1. 
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Test experience Indicates that CL1, CL2, CL3, and CL4 provide fast response, 
accurate measurements of liquid level during fill or expulsion. Submerged 
helium pressurization causes a very slight (±0.5 Inch) variation In Indicated 
liquid level due to the bubble and surface disturbances of the liquid hydrogen. 

4.9 TEST I - TANK LEAKAGE - CRYOGENIC 

The leak test of the IDU was first performed on 13 November by pressurizing 
the IDU to 52 pslg, with the main tank vented to atmosphere. The main tank 
was partially filled with liquid hydrogen. Figure 4-3 shows the pressure decay 
starting approximately five minutes after the IDU was pressurized. The tem- 
perature of the IDU and main tank near the top of the IDU was approximately 
90°R, and constant during the test; the temperature beneath the IDU was 
approximately 40°R. The pressure decay was, therefore, due primarily to 
leakage. 

The loss corresponded to 0.3 lb of helium in 2,500 sec, for a loss rate of 
-4 

1.2 x 10 lb/sec at approximately 40°R, or 0.6 SCFM. This leakage is much 
higher than would be expected with the fittings used in the IDU lines and the 
nominal leakage of the PV valves would be expected to be 10 SCIM according to 
Parker data. Further, the mass spectrometer tests showed no significant leak- 
age of any of the connections. Therefore It was concluded that the leakage 
is through the Parker valves which could not be completely closed. 

On 12 December 1973, the second IDU leak test was initiated. The IDU was 
initially empty of liquid and purged with helium. The PV1 and PV2 valves 
were closed at ambient conditions and remained closed during the test. 

Following the main-tank fill to a level above the IDU (indicated by RL10), 
the main-tank pressure was increased co 50 psig and controlled at this level. 
The IDU was locked up at 3 psig and the pressure monitored for any increase. 
Temperatures inside the IDU were monitored to assure steady-state conditions 
were maintained. The pressure increased to 3,65 psig in 30 minutes, for a 
rate of approximately 3.51 x 10”^ psi/sec. The top temperature (T3) indicated 
an increase from 45.5 to 56°R, whereas the bottom temperatures were maintained 
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at approximately 40°R, The effective pressure Increase due to temperature 
Increase In the top portion of the 1DU would be of the order of 2,8 psi, 
assuming a linear temperature Increase from the bottom of the tank to the 
top. The actual temperature profile for a typical stratified ullage would 
not be linear over the IDU tank length since liquid in the main tank Is above 
the bottom of the IDU; however, the level is indeterminate. It can only be 
concluded that the pressure increase Is probably due primarily to the tempera- 
ture increase in the top region of the IDU, which Is the result of a lowering 
of the main-tank liquid level from boiloff. Neglecting the temperature 
Increase effect, the observed pressure increase of 0.65 psi would correspond 
to approximately 0.035 lb of hydrogen gas in 1,800 sec, for a rate of 
2 x 10 lb/ ec or 0.2 SCFM of hydrogen, which Indicates that the observed 
pressure increase is due to heat flux into the IDU, since valve closure was 
complete unoer ambient conditions. Also, the leakage based on the 10 SCIM 
nominal rate at 40°R, as well as the ambient leak test data, are much lower 
than che 0,2 SCFM required to produce the observed pressure increase. 

The DU was then pressurized to 50.2 psig and the main tank vented to 
atmosphere and the pressure decay rate of the IDU observed. Results are 
shown in Figure 4-3. Temperature readings were not obtained during the first 
30 minutes of the test due to a faulty connection in the vacuum chamber, but 
subsequent temperature readings indicated steady-state conditions were 
obtained. The pressure decay rate during the period from 12:13 to 12:33 
indicated a loss rate approximately the same as that observed during the 
13 November tests. 

4.10 TEST 0 - SCREEN BREAKDOWN, PRESSURIZATION, AND IDU REFILL 
The objective of this series of tests was to determine the breakdown point of 
the IDU liquid hydrogen acquisition device (LHAD) for various flowrates, 
including "static breakdown" conditions, for both cold helium submerged and 
overhead pressurization, as well as cold and ambient overhead hydrogen 
pressurization. A series of tests at a maximum flowrate of 7.5 Ib/sec was 
scheduled. Because of the difficulties encountered with the Parker submerged 
liquid hydrogen valve, the test plan was modified, and high flowrate tests 
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could not be performed. However, except for the problems encountered with 
sporadic valve opening and closing malfunctions, the basic response of the 
system was as expected. I DU outflow and breakdown levels indicated th»t the 
screen device static bubble point was not adversely affected by the tests, 
successful screen cevices ( LHAD) refill was accomplished following breakdown, 
and, with one exception, all instrumentation functioned properly. A carbon 
point-level sensor (RL7) (located inside the LHAD annulus at the top of the 
device, between two of the 1-inch diameter tubes from the central pipe) did 
not properly indicate breakdown in some cases, apparently because of liquid 
"hang up." However, breakdown readings were obtained in all cases because of 
a characteristic hump in the capacitance continuous-level sensor curve as well 
as measurements of temperature increases within the LHAD. 

Details of each outflow test are summarized in Table 4-1 and discussed below. 
Detailed annotated data are presented in the Supplemental Data Document and 
examples of test data are presented in Appendix H. Selected data for breakdown 
level at various outflow rates are compared in Figure 4-4 with the analytical 
predictions obtained from the flow-loss computer program of Appendix F. It 
should be noted that the screen device performs better than predicted in the 
majority of tests with the cold pressurant, as shown, for example by Test J-19. 
The use of warm pressurant, however, resulted in premature breakdown. 

Certain tests were not performed due to the vilve malfunctions, and data on 
flowrate versus pressure drop across the valves is of little use since the 
valves would not open fully. Another difficulty encountered which limited 
the rate of outflow from the IDU was the limitation on the rate of venting the 
main tank. The main-tank vent (Port B of Figure 2-15) was used as a conduit 
for the electrical wires, and the reduction in flow area limited the rate at 
which the main tank could be vented. As a result, the initial outflow rate 
from the IDU to the main tank was higher than the final rate before breakdown, 
for the high flowrate cases. For example, the initial outflow rate for 
Test 0-25 was 1.5 lb/sec, whereas the final outflow rate was 0.73 ib/sec. 
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Different methods of IDU pressurization were used to evaluate the effects on 
screen device performance. The methods used for IDU pressurization werei 

• Cold submerged helium 

• Cold overhead helium 

• Warm overhead helium 

• Cold overhead hydrogen 

• Warm overhead hydrogen 

Tests J-25A, J-25B, J-33A, and J-33B showed that breakdown occurred near the 
bottom of the LHAD (2.25 to 2.85 Inches above the bottom of the tank) with a 
submerged screen depth no greater than 0.5 Inch, for flowrates from 0.5 to 
0.73 lb/sec. The helium was passed through the pressurization gas heat 
exchanger located below the IDU, and additional cooling was achieved by the 
gas bubbling up through the LH 2 in the IDU. Heat transfer effects on break- 
down were of negligible importance in these cases. Since the flowrates 
obtained were relatively low, viscous and dynamic pressure losses were low, 
and hence breakdown occurred very close to the bottom of the LHAD screen. 

Test J-25B shows a spurious Indication of breakdown at 6 inches, as determined 
by the "hump" in the capacitance probe liquid height versus time curve, but 
this is probably not an actual breakdown since expulsion continues to the 
2.5-1nch level. 

Tests J-3, J-ll, J-19, J-27A, J-27B, J-35B, and J-35C also showed that break- 
down occurred near the bottom of the IDU, for flowrates ranging from 0.25 to 
0.78 lb/sec. Tests J-35B and J-35C showed that helium temperatures of the 
order of 60°R, as measured by temperature probe T4 inside the LHAD (at the top 
of central column), did not cause an observable degradation in breakdown 
level of the LHAD. Since T4 is inside the LHAD, its measured gas temperature 
wcuid be less than the gas temperature surrounding the LHAD. The outflow 
periods were of the order of one minute; it is not clear what the effects on 
the screen device would have been with longer periods of exposure of the LHAD 
to the 50 to 60° R warm gas. 
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The data obtained for the overhead and submerged cold helium pressurization 
Indicates that either of these is satisfactory for Ll^ outflow. However, 
at higher outflow rates, the higher mass flow of helium entering the IDU 
through the submerged diffuser would increase the voltage fluctuation or 
“noise" obtained from the continuous capacitance level sensor, as well as 
distorting the readings from the capacitance point- level sensors. Therefore, 
use of the overhead diffuser is recommended, with cold helium. 

Tests with both warm gaseous hydrogen and warm helium pressurant showed a 
marked decrease in effective screen bubble point. As expected, the heat 
flux caused screen drying and premature breakdown. The warm gaseous hydrogen 
(autogenous) pressurization caused a greater reduction in effective screen 
bubble point than warm helium. This may be due to the higher heat-transfer 
coefficient for hydrogen condensing on the upper solid portions of the LHAD 
compared with the lower free convection heat-transfer coefficient for the 
helium near the top of the LHAD. 

Test J-35A, with helium temperatures ranging from approximately 70° to over 
100°R, showed breakdown at a liquid level of 6.4 inches above the bottom of 
the IDU. By contrast. Test J-5C, with 56°R gaseous hydrogen, showed break- 
down at 6 inches; Tests J-7, J-31, and J-39, with approximately 100°R gaseous 
hydrogen, showed breakdown at approximately 15 inches. 

These tests do not prove that warm pressurant cannot be used with screen 
devices since the LHAD has a relatively large area of exposed stainless steel 
annular channel to which the 250 x 1370 screen is welded. Heat transfer 
through this region could cause boiling and screen dryout that could perhaps 
be avoided by a design which had a minimum of exposed sheet metal surface. 

The LHAD was originally intended to supply liquid hydrogen, using cold helium 
pressurant, hence heat transfer effects were not an important consideration in 
its design. Consequently, the LHAD was designed to minimize fabrication costs 
rather than to optimize its performance with warm pressurant. 
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IDU refill test data are summarized In Table 4-2, and the annotated data from 
the strip chart presented In the Supplemental Data Document. The maximum 
Inflow rate obtained was 8.35 Ib/sec (Test J-26A) which was obtained with V4 
(see Figure 2-15) wide open. 

Two types of refill were demonstrated: (1) refill through PV2, which Is the 

procedure for "start tank" refill, and (2) refill through PV1 and PV3. The 
latter procedure "backfilled" the IDU by flowing liquid up the control IDU 
column, through the four communication tubes to the annulus. Refill directly 
through PV2 brought liquid into the IDU through the refill diffuser. Both 
methods refilled the IDU completely, as determined by RL7 and T4 data. In both 
cases, the screen port (see Figure 2-6) appears to have allowed gas inside the 
LHAD to pass cjt of the annular region until the liquid level, both inside and 
outside the annular region, covered the screen port, wetting and sealing it. 

If the LHAD had not refilled completely, T4 and/or RL7 would have indicated 
the presence of the gas, at least in some cases. Neither indication was 
observed, and therefore, it appears that either refill method can be used 
successfully. The preferred procedure is to refill through PV2. Refill 
through PV1 and PV3 results in flow out of the screen device, which places an 

Table 4-2 

TEST J - IDU REFILL 


Test No. 

IDU Initial 
Liquid Height 
(in.) 

IDU Final 
Liquid Height 
(in.) 

Refill 

Time 

(sec) 

Refill 

Rate 

(lb/sec) 

J-2 

18.5 

31.5 

42 

0.47 

J-6A 

4.2 

31.5 

38 

1.1 

J-6B 

4 

31.5 

24 

1.76 

J-10 

6.2 

31.5 

75 

0.52 

J- 18 

16 

31.5 

12 

3.24 

J-24 

10 

31.5 

28 

1.18 

J-26A 

20.5 

31 

2 

8.35 

J-26B 

5.8 

31 

28 

1.38 

J-28 

14.4 

31 

17 

1.5 

J-32 

2.6 

22.8 

16 

1.94 

J-34A 

5.8 

31 

26 

1.5 

J-34B 

9.6 

30.2 

25 

1.27 

J-38 

14.4 

31 

23 

1.11 
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outward pressure on the Inner cylindrical screen, which Is not supported. 
Therefore, continued refills In this manner at high rates could conceivably 
cause failure of the screen weldment to the supporting annulus, or bubble-point 
degradation. Further, refill through PV1 and PV3 can introduce gas into the 
LHAD, which is minimized by refill through PV2. It was necessary in our tests 
to both verify the refill procedure using PV1 and PV3 and to use this technique 

when PV2 would not open (and vice-versa). Refill through PV2 is a proven and 

preferred technique and should be used in general. 

The response of the continuous-level sensor to IDU refill was satisfactory in 
all cases. Refill through PV1 and PV3 resulted in a smoother output voltage 
versus time curve than for refill through PV2. Refill through PV2, even with 
the refill diffuser, results in surface sloshing, as shown in the water tests 
of Figure 4-2. As the liquid level in the IDU rises, the surface disturbance 

decreases. In the worst case (Test J-18), a voltage fluctuation of ±0.1 volts 

was observed, which corresponds to a measured liquid-level variation of 
±0.75 inch. Most of the data showed a liquid-level variation of less than 
±0.5 inch. 

4.11 TEST K - EXPULSION/REFILL 

The objective of this test was to demonstrate the complete IDU sequence of 
operations involved in the start tank procedure for engine start, refill, and 
continued expulsion. The maximum expulsion outflow rate obtainable through 
the relatively small (1 in.) port of the main tank (with the additional flow 
losses associated with the liquid hydrogen vent system) was 0.75 Ib/'sec. The 
complete outflow and refill procedure was demonstrated, as shown by the data 
of Test K-2 of the Supplemental Data Document. Outflow through PV1 was 
initiated with the IDU full, it continued to a liquid level of 5.8 inches, 
where PV2 was opened and refill, along with outflow, continued until the main 
tank was depleted of liquid hydrogen. The IDU was refilled to a level of 
24 inches. Breakdown did not occur during this test. 
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4.12 TEST L - VALVE SEQUENCING BY OPERATOR 

The objective of this test was to determine whether the test conductor could 
manually actuate the PV valves during IDU expulsion before breakdown for 
refill or for halting expulsion. Because of the difficulties with the PV 
valves, this test could not be performed as planned. Experience at the flow 
rates of approximately 1 to 2.5 lb/sec (e.g., Tests N-l and N-12) revealed no 
difficulty with manual flow control, and It Is expected that manual control 
could be used for flow rates of the order of 4,5 Ib/sec. 

4.13 TEST M - THERMODYNAMIC VENT SYSTEM FLOWRATE VERSUS PRESSURE 

The objective of this test was to determine the TVS mass flowrate as a 
function of pressure for the various combinations of viscojet flow control 
orifices, and to compare the results with analytical predictions. These 
tests were performed with a TVS line pressure of approximately one atmosphere, 
and an IDU pressure of approximately two atmospheres, 

Tests were performed on 14 November 1973, The IDU and main tank were partially 
full of LH 2 ; the IDU was pressurized with cold helium through the overhead 

diffuser. The TVS flow control valves, SV1, SV2 SV6, shown in 

Figure 2-15, were opened in a sequence that allowed the mass flowrate to 
increase with each step, except in three cases. In all cases in which the 
successive opening of a valve led to an expected flowrate increase, the flow- 
rate agreed well with predicted results. However, in those cases where certain 
valves were opened and closed to decrease the flowrate, the measured outflow 
rate would stabilize at a higher value than expected and in the case of 
Test M-27, not reach a steady flowrate for periods in excess of 20 minutes. 

This inconsistency was not noted during Tests M6 and M7 since the flowrates 
appeared to represent steady-state values for periods of 5 to 10 minutes. In 
restrospect, it is seen that the flowrates in these cases were high due to the 
presence of additional liquid in the TVS tube which was evaporating, and con- 
tributing to the total flowrate almost as if additional flow control valves 
were open. The time to reach equilibrium would be shorter if the heat flux to 
the TVS were higher. In these tests, heat flux to the TVS was quite low since 
the IDU was submerged in LH 2 over most of its height. However, heat-flux data 
was not taken during these tests. 
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Data are compiled In Tables 4-3 and 4-4 and test results are shown In 
Figure 4-5 where the analytical results are given for comparison. 

Using data from Reference 10, the correlation of flowrate and pressure drop 
for saturated (x>0) versus subcooled hydrogen (x=0) In the viscojet was found 
to be improved If the effective density of the two-phase fluid was used 
Instead of the liquid density. The average effective density of the two-phase 
fluid in the viscojet is given by 

p TPF = ^ " $ p LH 2 + T p LH 2 (4_1) 

where x is the vapor quantity at the viscojet exit. Figure 4-6 shows data for 
the VJ-1 viscojet (P/N 38VL1-CM) which was useu in earlier MDAC I RAD tests. 
Agreement between the subcooled and saturated conditions (inside the viscojet) 
are seen to compare well. Over the range of conditions tested with the IDU/ 
TVS, data were found to correlate adequately using pressure drop times liquid 
density versus flowrate, since the quality was less than 0.1 in all cases. 

As discussed in Appendix D, the flow resistance of the viscojets is given 
in terms of LOHMS, with 


L (LOHMS) = 1270 - n/aP p , 

• L 


pi corresponds to single 
phase, subcooled fluids 


(4-2) 


where LOHM is a flow resistance unit corresponding to the flow of 100 gallons/ 
minute of 80°F water with a pressure drop of 25 psi. The factor f is required 
for specific fluids. Previous tests (Reference 13) showed this value to be 
f * 0.823 for liquid hydrogen. 


The data confirmed the expectation that total flowrate for multiple parallel 
paths is given by. 


m = 1 ,270 f 



since the flowrates from viscojets in parallel are additive. 
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Table 4-3 

THERMODYNAMIC VENT SYSTEM FLOWRATE VS PRESSURE 
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Table 4-4 

THERMODYNAMIC VENT SYSTEM FLOWRATE VS IDU AND TVS PRESSURE 
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M-32 SV1,SV2,SV3,SV4, No data since all liquid 

SV5,SV6 drained from IDU 
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Figure 4-6 
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Where o TpF Is a more general correlation parameter than p^. Equation 4-3 Is 
recoranended for predicting flowrate as a function of uallty, x, and pressure 
drop across the viscojet. In the analysis of Appendix I, the pressure drop 
across the viscojet Is shown to be two orders of magnitude greater than the 
pressure drop along the TVS tube (see Figure 0-1), and therefore, the pressure 
difference between the I0U ullage and the TVS exit should adequately represent 
the pressure drop across the viscojets. However, no measurements were made 
of the pressure drop In the TVS line, and therefore the measured total pressure 
drop may exceed the actual pressure drop across the viscojets alone. This may 
explain the fact that some of the data show a smaller than expected maximum 
flow rate for the measured pressure drops. 

4.14 TEST N - HYDRAULIC PRESSURE SURGE 

Exploratory tests were performed to determine the effects of pressure 
transients Induced by valve opening/closing on screen device retention capa- 
bility, Five tests were performed In which valve opening/closing was partially 
successful. Numerous tests were attempted, without success, due to valve mal- 
function. In particular, PV1 failed to open in the majority of cases. 

Because of the relatively low flow rates (1 Ib/sec), the pressure surges 
encountered at the screen were relatively low. The expected pressure ,e 
level was analyzed as a part of the parallel effort of NAS8-27685, ;nce 1. 

The MDAC Liquid Propulsion Feed System Dynamic Analysis Program, H 172, was 
used to determine pressure and flow conditions caused by valve opening and 
closing. The program uses a nonlinear description of components at discrete 
junctions with solutions to the flow between junctions obtained by the method 
of characteristics. Using a mathematical model of the IDU, results were 
obtained for pressure transients of both a rigid and compliant annular screen 
device. For a steady-state flow rate of 6.5 Ib/sec, the pressure excursion 
at the top of the screen device was found to vary from 1 psi above the ullage 
pressure to 2 psi below the ullage pressure assuming a valve opening time of 
100 Ms. However, the maximum flow rates obtained during the test was 1 Ib/sec, 
and since tne pressure surge Is proportional to the steady-state velocity of 
the flowing fluid, the maximum pressure surges at the screen would be expected 
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to be of the order of 0.15 to 0.3 psl. Attenuation of the pressure transient 
due to screen compliance would be expected to further reduce the pressure 
surge. The critical condition for serein device retention capability is the 
valve-opening induced rarefaction wave effect, rather than the valve closing 
Induced compression wave effect since the resulting pressure drop inside the 
screen device can cause gas Ingestion through the screen (Reference 14). The 
total flowrate out of the IDU was not necessarily through the screen device 
and PV1. Since the refill valve often gave indications of incomplete closure, 
some of the outflow was undoubt'. Ily through PV2, Therefore, the pressure 
decrease due to PV1 opening would be less than if PV2 were not leakinq. In 
spite of these difficulties and low flowrates, there were three cases in 
which opening of PV1 was followed by an apparent premature breakdown. 
Oscillograph data for the tests are given in the Supplemental Data Document. 
Results are summarized in Table 4-5 Tests N-l, N-2, and N-3A show breakdown 
occurring 0.5 second after PV1 is opened, but the liquid level cannot be 
determined accurately near the bottom of the tank, because of the nonlinearity 
in voltage versus liquid level in this region. Therefore the results are 
questionable. 

Opening and closing PV1 with the IDU liquid level 8 to 20 inches above the 
plate never resulted in breakdown at the flowrates obtainable, and in test. 
N-11 breakdown did not occur on opening PV1 with the liquid level at 
13.5 inches and an apparent outflow rate of 1.04 Ib/sec. This flowrate is 
sufficiently high to cause an unattenuated pressure decrease of the order of 
0.3 lb/in. , which could cause noticeable premature breakdown. The absence 
of pressure transient induced breakdown significantly preceding the nominal 
breakdown level could be explained by (1) valve closure is sluggish, thus 
decreasing the pressure transient (2) the vapor region in the top of the I.HAD 
attenuates the wave, or (3) the actual feedline outflow rate was less than 
the apparent outflow rate due to leakage from the ,'DU to the main tank throuqh 
PV2. Structural compliance could not contribute significantly to the wave 
attenuation since the wave travels through the relatively ridig feedline, 
central IDU pipe, across the four arms of the ’DU, and to the top of the 
screen device where breakdown occurs. 
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Table 4-5 

HYDRAULIC PRESSURE SURGE TESTS 
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4.15 TEST 0 - TVS FLOWRATE WITH VARIABLE TVS LINE PRESSURE 
Based on the analysis of TVS flowrate versus pressure difference, it would not 
be expected that TVS line pressure would have a first-order effect on flow- 
rate for a given pressure difference. Tests were performed to confirm this 
expectation; the TVS line pressure was controlled at values ranging from 
5 to 15 psia, and flowrates measured and compared with theory. Originally, 
these tests were to be performed at 5 psia only, but two considerations led 
to varying the TVS line pressure over the stated range. First, for steady- 
state thermal isolation ( i . e. , heat interception), the TVS is required to 
maintain both a given flowrate and a given temperature inside the TVS line. 
Maintaining a given temperature of the coolant implies that a given, constant 
pressure be maintained. Based on the analytical results shown in Appendix G, 
in most applications of the IDU this pressure will be of the order of 10 psia, 
and a valve downstream of the TVS coil will be required to throttle the flow 
to the required TVS line pressure. 

A second consideration was that the vacuum pumps available for the tests could 
not maintain 5 psia in the TVS at the higher flowrates, and therefore, some 
variation in TVS line pressure was unavoidable. 

The TVS schematic drawing of Figure 2-15 shows that line pressure is controlled 
by needle valves V23 and V24. These were adjusted, as necessary, to vary the 
line pressure, and hence partially control flowrate for selected viscojets and 
a constant IDU pressure. 

Data obtained from the tests is shown in Table 4-6 and plotted in Figure 4-5. 

The dependence of flowrates on the square root of pressure difference times 
the effective two-phase fluid density is partially confirmed within the range of 
conditions tested, but the flowrate obtained for a given total pressure drop 
(including the TVS line pressure drop), is lower than predicted. This may be 
due to the two-phase pressure drop through the TVS being higher than predicted 
by the analysis of Appendix D. It would be expected that the pressure drop 
through the TVS would be greater with the lower line pressure and hence higher 
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Table 4-6 

THERMODYNAMIC VENT SYSTEM FLOWRATE VS IDU AND TVS PRESSURE 
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quality obtained In Test 0. Since there were no provisions for measuring 
TVS line pressure at the viscojet exit, this could not be confirmed. The 
TVS line was bypassed In Tests 0-4, 0-5, and 0-10 In an attempt to evaluate 
pressure drops directly through the viscojets. Flow proceeded directly to 
the overboard vent, as shown in Figure 2-15, through a l/2-1nch outside 
diameter line. The flowrate In these cases was also low compared to the 
analytical value. Further, no comparison test was performed, using the same 
viscojet (e.g., VJ4, or VJ5) and flowing through the TVS line, and therefore, 
an approximation to the TVS line pressure drop cannot be made. 

The test data obtained are sufficient to empirically determine flowrate vs. 
pressure, and the analytical prediction for viscojet flowrate/pressure depen- 
dence still appears to be valid, although the actual pressure drop across the 
viscojets is not known. It is therefore recommended that at least one addi- 
tional pressure pick-off point be included immediately down stream of the 
viscojets in any future tests. 

The data obtained is particularly useful since the TVS line pressures corres- 
ponds to values required for steady state thermal control of the IDU. 

4.16 TEST P - STEADY-STATE THERMAL CONTROL OF IDU 

Thermal control of the IDU, or any type of similar tank with a wall mounted 
TVS heat exchanger, can be accomplished by one of two basic methods: 

(1) steady-state operation in which the total heat entering the propellant is 
zero at all times, or (2) transient operation in which the tank pressure is 
maintained within certain limits by intermittent flow through the TVS. 
Steady-state operation offers the advantage of minimum transient heat transfer 
and thermal control problems. (In principal, heat-flux meters can be used to 
monitor heat flux into the tank. Maintaining TVS outflow so that the overall 
net heat flux is zero allows rapid response to changing heat flux.) The 
thermal stratification in the fluid is minimized and the total mass vented is 
always equal to the amount required to intercept the heat into the tank. The 
pertinent boundary condition is that the total surface integral of the heat 
transfer, taken at the inner walls of the tank, be zero. This condition 
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Implies that the maximum temperature between the TVS colls not be equal to 
desired steady-state temperature of the liquid In the tank. Rather, the 
maximum temperature Is somewhat higher than the tank liquid. Thus, some 
heat enters the liquid In the vicinity of the midpoint between colls, whereas 
an equal amount of heat Is absorbed from the tank liquid near the colls so 
that no net heat enters the propellant. This procedure Implies that the maximum 
temperature of the tank liquid approaches the maximum temperature between the 
colls, The further implication Is that boiling will occur unless the total 
pressure in the tank exceeds the vapor pressure associated with the maximum 
temperature between the colls. Therefore, the coil separation distance will 
be governed by the maximum acceptable wall temperature (for a given heat flux). 
This condition on T mgx between the coils allows the total TVS tube length and 
spacing to be determined. Otherwise, a wide variety of combinations of tube 
length, spacing, and coolant temperatures are possible. 

The IDU has been designed so that the steady-state operation can be demonstrated 
with two coil separation distances. This test flexibility is obtained by the 
addition of a "Tee" connection in the IDU line, resulting in two spiral loops 
covering the tank length of 35 ft and 70 ft. In the current test plan, the 
two 35-ft loops are used in parallel. Additional capability can be obtained 
by the addition of an SV-90 solenoid valve and an extra TVS outlet line; it is 
recommended as a desirable addition for any subsequent tests. 

The main objective of this test was to demonstrate steady-state thermal control 
by flowing sufficient coolant, at the required TVS line pressure (and hence 
coolant temperature) to maintain complete heat interception. Heat flux from 
ihp. main tank was to be sustained by maintaining a gaseous hydrogen ullage 
pressure slightly higher than the saturation pressure of the liquid in the 
IDU, with the result that condensation heat transfer would occur on the IDU 
tank wall and coils. It must be emphasized that the IDU/TVS is designed so 
that an external foam should be applied to the tank wall to alleviate conden- 
sation heat transfer, but that the foam could not be applied until after these 


94 


i 


i 



Initial tests. The presence of a relatively cold TVS coil with a surface 
temperature several degrees below the saturation temperature of the main tank 
GHg results in a relatively high, localized condensation heat-transfer coefficient; 
this results in high overall heat flux to the coolant, which tends to diminish the 
effectiveness. Insulation, however, would alleviate condensation and the TVS 
would then be able to maintain a steady-state thermal condition without the 
additional complexity of local, high-heat flux to the TVS tube itself. 

In these tests, it was originally intended that the IOU be filled with liquid, 
the PV valves closed, and the proper TVS coolant flowrate and temperature 
adjusted to maintain a steady-state temperature and pressure inside the IDU. 
However, the PV valves could not be closed completely, in spite of repeated 

attempts, and hydrogen continuously leaked from the IDU to the main tank. 

As a partial demonstration of thermal control, it was decided that the main- 
tank liquid level would be set at the bottom of the IDU, and the main-tank 
pressure controlled to approximately 0.2 to 0.3 psi above the IDU pressure. 

The IDU pressure would be vented periodically to maintain the pressure dif- 
ference. The net result, it was hoped, would be that liquid in the main tank 
would tend to flow into the IDU, maintaining a full tank during the test, while 
the higher external 6 H 2 pressure would result in heat transfer to th? IDU. The 
TVS coolant flowrate would be adjusted to maintain steady state. 


The test was performed as described for over 1 hour, and an overall thermal 
steady state was demonstrated. The measured heat fluxes showed a higher heat 
flux into the tank near the top of the IDU, and an approximately zero net 
heat flux into the IDU near the bottom of the IDU. The results are shown in 
the data plots, taken directly from the DYMEC Data System, in Figures 4-7 to 
4-9. Due to leakage of liquid into and out of the IDU, as well as the IDU 
venting, the total liquid lost exceeded the expected liquid loss based on the 
measured flowrate and total test duration. Table 4-7 compares the measured 
and expected liquid levels. 
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Figure 4-9 -TEST P-1 f PLOT NO. 3-1 
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It Is recognized that these test results are of limited value from a heat 
transfer standpoint. However, a steady-state coolant outflow rate was 
demonstrated, a relatively steady-state heat flow was observed, and all 
IDU/TVS Instrumentation was observed to function properly. 


4.17 TEST Q - EFFECTS OF DOWNSTREAM FLOW CONTROL 

A basic question concerning the operation of the TVS Is the method by which 
coolant flowrate and temperature can be varied to meet the thermal control 
requirements and the TVS response. A variable downstream orifice can be used 
to control the TVS line pressure and hence coolant temperature, while flowrate 
can be controlled by varying the selected orifice and the I DU pressure with 
respect to the TVS line pressure. When steady-state conditions are not main- 
tained, transient flow effects occur. The objective of this test was to 
demonstrate some of the transient phenomena resulting from (1) closing the 
downstream flow control valve, V24, thus halting the outflow, (2} achieving 
a steady state, and then (3) reopening the valve. Immediately after Test P, 

V24 was closed, but the solenoid valves used to provide a controlled flowrate 
of 2.5 lb/hr were left open. The IDU and TVS conditions were then observed 
until a steady state was reached. The steady state was observed for approxi- 
mately 1 hour, and then the TVS flow was reinitiated and conditions monitored. 
Throughout this period, the main tank and IDU pressures were maintained at 
25 psia. The IDU liquid level was approximately 16 inches. 

Referring to the selected data of Figures 4-10 to 4-16, the phenomena observed 
are described below in chronological order. At 12:45:00, valve V24 (see 
Figure 4-10) was closed, halting all TVS outflow, and entrapping some subcooled 
liquid in the TVS tube. For a period of approximately 10 to 15 minutes, the 
coolant in the upper part of the TVS coll was vaporized by heat flux from the 
main tank ullage, and the vapor and liquid forced back into the IDU. 
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Figure 4-12 - TEST Q-l, PLOT .40. 3-1 
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Figure 4-15 - TEST Q-2, PLOT No. 3 
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Figure 4-16 MAIN TANK TEMPERATURE, T5, DURING TESTS Q-l AND Q-2 


During the first 4 minutes, the TVS inlet temperature, T7, increased from 
40°R to 41.0°R, and remained approximately constant (+0 5°R) during the entire 
steady-state period. The TVS outlet temperature also remained constant, at 
approximately 41.0 ± 0.5°R for the 1.2 hour steady-state period. The tempera- 
ture, T14, immediately downstream of viscojet orifice FJ5 increased in the 
first six minutes from 29° to 41.5°R, and regained approximately constant the 
the 1.2 hour period. Thus, the TVS equilibrated in approximately 5 to 
6 minutes. The IDU ullage temperature (T15) measured at HF2-A required approxi- 
mately 20 minutes to reach a stpady-state condition, increasing from 40°R to an 
approximately constant 44°R, The ID'! ullage temperature (T16) measured near 
the top (HF3-A location) was warmer, increasing in approximately 15 minutes 
From 45° to 52°R (Figure 4-11). The IDU outer-wall temperature near the top 
also warmed slightly, from 36° to 41°R in approximately 10 minutes. During 
this initial period, the heat flux into the IDU near the bottom (HF1-B) was 
approximately zero, since both sides of the IDU wall were immersed in liquid. 

The heat flux near tne middle of the IDU was measured bv the interior heat- 

2 

flux gauge, HF2-A, was approximately 1.5 Btu/ft hr, and was relatively 
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constant. The heat flux near the top of the I DU, as measured by the interior 

2 

heat-flux meter, HF3-A, showed a rapid decrease from 10 to 15 Btu/ft hr to 

2 

approximately 1.0 ± 0,5 Btu/ft in approximately 15 minutes. Recall that this 
heat-flux meter is positioned to straddle the position of a TVS co^i I'ccted 
on the exterior of the I DU cover. Positive heat flux is directed from the 
IDU to the TVS coil, and therefore, the cold TVS fluid, located at the top of 
the IDU, was cooling the interior of the IDU as it vaporized and flowed ^ k 
into the IDU through the VJ5 orifice. The heat flux measured by HF4-B, 
located between the TVS coils on the outside of the IDU cover, showed a 
relatively low heat flux during both the initial and steady-state periods. 

After the steady-state period had been observed for over 1 hour, the TVS 
downstream valve (V24) was reopened, and within 1 minute a steady outflow 
rate of 2.5 lb/hr was incurred. Within approximately 6 minutes, all IDU 
temperature measurements showed a marked decrease to steady-state values, 
which were maintained for approximately 30 minutes, as shown by Figures 4-13 
and 4-14, for Test Q-2. 

Figure 4-14 shows that the heat flux near the top of the IDU, measured by 

2 2 

HF3-A, climbs rapidly from approximately 1.0 ± 0.5 Btu/ft hr to 9 Btu/ft hr. 
This is, of course, due to the TVS cooling the IDU interior. Similarly, in 
Figure 4-15, HF4-B shows that the heat flux from the main tank to the top of 
the IDU cover (as measured by the exterior heat-flux gauge located between the 

n 

TVS coils) increases from approximately 1 to 2 Btu/ft 6 hr to approximately 
2 

8 to 9 Btu/ft hr. Again, the TVS coolant, having lowered the IDU wall tempera- 
ture, has caused a significant increase in heat flux to the coil from both 
sides of the IDU. Figure 4-16 shows that the main tank ullage temperature 
momentarily decreased immediately after TVS flow is reinitiated, but the 
incoming warm GH^ soon brought the ullage temperataure back up to 70°R. As 
predicted by the TVS analysis, the uninsulated IDU/TVS causes an increase in 
the overall heat flux by increasing the temperature difference between the main 
tank and the IDU, but more importantly, by causing condensation heat transfer 
to occur, rather than free-convection heat transfer. Tests Q-l and Q-2, there- 
fore, demonstrate the need for insulation on the IDU exterior if the TVS heat 
interception capability is to be fully exploited. 
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4.18 TEST R - TVS TRANSIENT OPERATION 

The principle of transient TVS operation Is that heat entering the propellant 
tank is allowed to induce stratification until the total tank pressure exceeds 
a predetermined value, at which time coolant flow is initiated. An immediate 
concern is the response of the ullage pressure to the initiation of chilldown 
flow. This test was to demonstrate the response rate for an external heating 
condition similar to that demonstrated with Test P with the IDU approximately 
filled with liquid. IDU pressure would be allowed to increase until a facility 
pressure switch exceeds a set value (e.g., 28 psia). At this time, the TVS 
valve for maximum flowrate would be opened, and flow maintained until IDU 
pressure decreases to a previously determined lower value (e.g., 20 psia). 
However, tank leakage due to partial closing of the PV valves could not be 
alleviated, and therefore, this test was scrubbed following repeated unsuccessful 
attempts to induce complete valve closure. 

4.19 TEST S - AUTOGENOUS PRESSURIZATION 

The original objective of Test S was to evaluate the response of the IDU to 
warm gaseous hydrogen prepressurization, without outflow, and to determine the 
pressurization collapse factors for two cases: (1) the screen device (LHAD) 

full of liquid, and (2) the '..HAD empty. Heat transfer out of the IDU was to 
be minimized by evacuating the main tank. However, because of the Parker 
valve malfunctions, this test was modified to further determine the effect of 
warm gaseous hydrogen pressurization on screen retention. 

The tests were run in the same manner as Test J, with warm hydrogen pressurant 
exclusively. Outflow from the IDU to the main tank occurred in every case, 
due to incomplete closure of PV1 and/or PV2. The IDU was pressurized to 
40 psia in all tests. 

Table 4-8 summarizes the results. The data of Tests S-1A and S-1B show that 
vapor is formed inside the LHAD before breakdown, due to exposure of the upper 
part of the LHAD to the warm pressurant. Note that the initial liquid level 
in Test S-1A was 14 inches, thus initially exposing the upper portion of the 
LHAD to warm pressurant. In Test S-1B, the LHAD is submerged in liquid 
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Initially, and vapor formation occurs 20 sec before breakdown, as compared to 
40 sec for Test S-1A, In Test S-1A, breakdown does not occur until the liquid 
level is 8,6 inches. This is significantly above the breakdown level expected 
with cold pressurant, but shows significantly higher retention capability than 
Test S-1B, with breakdown occurring at a liquid level of 20.5 inches, which is 
immediately below the top part of the screen device, Test S-1D gives results 
similar to S-1B; breakdown occurs at a level of approximately 21 inches, but 
in this case the IDU was initially filled and therefore the time the LHAO was 
exposed to warm ullage before breakdown was only 5 seconds. This extremely 
rapid, and unusual, breakdown response is not seen in other tests; it can be 
surmised that the LHAD was not completely filled with liquid initially, and 
the presence of vapor in the top region led to immediate breakdown. Tests S-1F 
and S-1G again show premature breakdown, due to warm gaseous hydrogen, but no 
vapor is indicated inside the IDU by either T4 or RL7. 

In Tests S-1A and S-1B, breakdown is explainable as being due to heat flux 
through the solid central pipe, communications ducts, and cylindrical channel 
to which the screen is attached. The heat flux would be expected to vaporize 
the LH£, filling the upper region with vapor, and thus causing breakdown by 
exposing the screen. However, in Test S-1F, there is no indication of vapor 
forming in the LHAD before breakdown; in fact, T4 indicates vapor 5 seconds 
after breakdown. In Test S-1G, T4 indicates vapor only 1 to 2 seconds before 
the CL4 indication ( i . e. , no further outflow and liquid dropout from the LHAD). 
These two tests and especially Test S-1F, appear to represent a degradation in 
screen retention capability due to heat flux to the wetted screen. 

In all cases, however, the use of warm gaseous hydrogen pressurization caused 
a very pronounced premature loss of liquid retention which appears to be due 
to heat flux. A discussion of a possible model explaining the loss of reten- 
tion for a column of liquid hydrogen supported by a screen exposed to heat flux 
is given in Reference 1, and will not be repeated here. The model does present 
a possible reason for screen device retention failure, independent of the vapor 
formation mechanism due to heat flux directly through solid metal walls, such 
as the top of LHAD, above the screen weldment. 
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4.20 TEST T - FEEDLINE VAPOR CONTROL 

The Interface between the main tank or a start tank and the feedline poses a 
potential problem with screen devices due to vapor generated In the feedllne, 
(above the feedllne valve) passing Into the screen device Iteself. A feed- 
line vapor barrier screen has been included In the IDU design (see Figure 4-17) 
which alleviates this problem by allowing vapor to bypass the screen barrier 
Into the IDU, thus maintaining liquid above the screen. 

The operation of the feedline vapor barrier screen as currently incorporated 
in the IDU is described below. This screen is used to keep gas from passing 
up into the screen device with vaporization in the feedllne above the feedline 
valve. In principle, it would keep the vaporized hydrogen from breaking 
through the screen even if the feedline were completely dry (i.e., all feed- 
line liquid boiled off). 

Consider the vapor barrier screen and bypass line shown in Figure 4-18A. When 
gas is being forced out through the elbow, but the vertical portion is full 
of liquid, then the pressure at the elbow is 

P L = P T + P9 [h LIQ + H 0UT + h ELBOW ^ 

which is approximately equal to the feedline pressure, for the low flowrates 
involved. 


A pressure immediately above the screen is 


P SCREEN “ P T + pg h LIQ 

The liquid remains above the screen, and no gas breaks through, if the bubble 
of pressure is, 


AP BP 2 P L " P SCREEN ~ P9 (h 0UT + h ELBOW ^ 


Thus, h^ujQy must be minimized, which has been accomplished in the tube routing 

of the IDU. To test this device, GHe was introduced into the feedline up to 

*> 

the screen by backflow through the feedllne from a 1.5 ft helium supply bottle. 
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The GHe flow was bypassed through the solenoid valve, SV10, between the cone 
screen housing and the bottom of the IDU (see Figure 2-15). Approximate 
flowrates were measured by noting the 1.5 ft bottle temperature and pressure 
as a function of time. 

A flowrate of 0.004 lb/sec was observed to have no effect on screen retention, 
nor was there any Indication of vapor passing through the cone screen Into 
the LHAD. The flowrate used corresponds to extremely high equivalent heat 
transfer rates into a feedline, assuming all the heat input causes boiloff. 

For a 10-ft section of feedline, 4 inches In Hiameter, a flowrate of 

2 

hydrogen gas of 0.004 Ib/sec corresponds to a heat flux of 275 Btu/ft hr, 
which far exceeds the values easily obtained with foam or HPI-insulated 

feedlines. Even helium purged HPI under ground-hold conditions has heat- 

2 

transfer rates of the order of 200 to 300 Btu/ft' hr, and therefore, the 
screen cone device appeals to be a feasible approach to vapor control under 
the most severe conditions. 

A useful modification to the IDU has been conceived which could be incorporated 
if considered necessary. This modification ensures that the IDU feedline would 
be full of liquid before feedline valve opening and flow to the turbopumps. 

If a separate line were attached to the inlet which allowed gas to pass 
directly into the ullage, conditions would be different, as shown in 
Figure 4-18B. 

Considering Figure 4-18B, the pressure above the screen is, 

^SCREEN = P T + pg h LIQ 

The pressure beneath the screen is P^. Thus, 

^SCREEN " p t " pg ^lIQ’ which is positive, 

and liquid flows down through the screen, keeping the feedline full of liquid. 
This alternate design may serve as a useful addition to the IDU, since vapor 
filling the feedline would normally pass through the pumps, possibly causing 
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damage. As a precautionary measure, the feedline vapor bypass valve could 
be opened momentarily, before Initiation of flow. Thus, any vapor In the 
feedllne would be displaced by liquid. This modification to the IDU would 
require either: 

A. An extension line to the Inlet, as snown In Figure 4-18B, which 
wouid then eliminate the capability for maintaining a nonwetted 
feedllne, or 

B. An additional solenoid valve and extension line, so that the feedllne 
could be maintained full of gas or liquid, depending on the desired 
conditions. 

4.21 TEST U - TWO-PHASE REFILL 

Two-phase refill with high gas volume percentages was simulated to assess any 
limitations of the LHAD. In particular, It was necessary to determine 
whether or not the gas Inflowing with the liquid caused the capacitance probe 
to give spurious liquid-level measurement,, and whether or not liquid would 
be lost through the vent due to the bubbles rising in the bulk liquid. T n e 
tests demonstrated that complete settling and displacement of ail gas hubbies 
from the bottom of the main tank Is not necessary before refill of a start 
tank since no adverse effects were observed. Thus, decreased propellant 
settling times can be used, which decreases the required start tank size and 
weight, and allows a constant, relatively short duration settling time to be 
used with variable propellant percentages in a given tank, rather than 
settling times programmed to the specific settling conditions encountered. 

Two refill tests were performed, U-l and U-2, and data are shown in the 
Supplemental Data Document. In both tests, refill was initiated through PV2, 
with simultaneous helium flow through the submerged diffuser to simulate 
two-phase flow Into the IDU. In Test U-l, the refill rate was a relatively 
constant 0.66 lb/sec, and In Test U-2, the refill rate was 0.805 Ib/sec. The 

3 

helium Inflow rate was approximately 1 ft /sec, which had no deleterious 
effect on refill. The capacitance point-level sensor showed maximum vol.age 
fluctuations of approximately 0.1 to 0.2 volts, corresponding to disruptions 
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of the liquid Interface of approximately 0.8 to 1.6 Inches. These maximum 
fluctuations appeared only at the Initiation of refill, and were strongly 
attenuated for liquid levels above approximately 20 Inches from the IDU bottom. 

4.22 TEST V - kETENTION WITH WARM LH2 

The objective of this test was to determine the screen retention capability 
witn a bulk liquid hydrogen temperature higher than that Investigated in Test J. 
Since surface tension decreases with temperature, the tests were expected to 
result In a decrease In retention capability, wMch was indeed observed. 

The original test plan called for a series of tests with various flowrates and 
bulk liquid temperatures, and both gaseous hydrogen and helium pressurization. 
However, due to the 2-inch (PV) ball valves malfunctioning, only five tests 
could be performed, all with cold helium pressurization (both overhead and sub- 
merged); results are summarized in Table 4-9. In Test V-25, the actual flow- 
rate through the LHAO am. PV1 is suspected to be less than the value calculated 
from the continuous level sensor data since PV2 may not have closed completely, 
thus allowing LH 2 to leak directly from the IDU to the main tank. 

The bulk liquid temperature of the hydrogen was controlled by the fill opera- 
tion. A constant IDU/Main Tank vent pressure was maintained during initial 
fill from the liquid r.ydrogen storage tanks in order to achieve the desired 
bu’k liquid temperature. This procedure was far more rapid and practical than 
waiting for the canked liquid to warm up from, say 37°R to 42°R, and it 
assured a well mixed, destrati fied propellant rather than a stratified 
propellant. 

The breakdown tests were performed as in the Test J cases, and the data arc- 
shown in Fiqure 4-4. The improvement, in terms of rettntion capability for 
the Test V cases is again seen, relative to the expected performance. The 
bulk liquid hydrogen was saturated at 35 psia (42.5°R) and, for example, a 
breakdown height of 7.6 inches above the IDU cover was expected at a flowrate 
of 0.8 Ib/sec, whereas in Test V-3, breakdown did not occur until the liquid 
level reached 4.2 inches, which is a significant improvement. Part of this 
improvement may have been due tc the cooling of the screen interface as 
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hydrogen was vaporized from the free surface of the liquid wetting the screen 
as the screen was suddenly exposed to the ullage gas. Part of this improvement 
was apparently due to the screen breakdown point pressure difference being 
higher than the bubble point pressure difference; i.e., a "catastrophic'' 
breakdown as opposed to single bubble breakthrough. Therefore, Figure 4-4 
shows both of the curves corresponding to the expected breakdown points based 
on the isopropyl alcohol bubble point of 14 inches of water column pressure 
head and the isopropyl breakdown point of 16 inches of water column pressure 
head. Part of this improvement may also have been due to the actual flowrate 
through the LHAD being less than the total flowrate as measured by the capaci- 
tance probe. However, even if the viscous and dynamic pressure losses were 
negligible (corresponding to a negligible flowrate through the LHAD) the hydro- 
static head alone should have caused failure at 5.8 to 7.2 inches. Therefore, 
the observed retention improvement had to be caused, at least in part, by 
effects other than leakage flow through PV2 decreasing the viscous and dynamic 
pressure losses. 

With helium as the pressurant, there is apparently a significant mass trans- 
fer effect taking place at the suddenly exposed liquid surface as the liquid 
drains out of the IDU. As this hydrogen is vaporized and transported into the 
ullage region by diffusion and convection, the free surface temperature would 
tend to drop below the initial bulk liquid temperature, assuming heat transfer 
into the IDU and screen device is low (as was the case for these tests). Thi 3 
effect was apparently present in the corresponding Test J conditions (i.e. 
cold helium pressurization) except that the liquid bulk temperature was 1 

The surface tension is determined by the temperature at the free surface, 
rather than the bulk temperature, and for the cold helium pressurization/ 
expulsion cases, equilibrium conditions were, of course, not met (i.e., tne 
partial pressure of the hydrogen was not constant and equal to the vapor 
pressure corresponding to a bulk liquid temperature). Since this vaporization 
would tend to cool the free surface, the surface tension would increase, thus 
improving the screen retention performance. Conversely, use of saturated 
hydrogen pressurant at a pressure exceeding the saturation pressure of the 
bulk liquid would r suit in condensation of a relatively warm liquid film, or 
condensate, which vr.uld have a rj'-’jced surface tension and, therefore, result 
in decreased screen retention performance. 
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Comparison of breakdown le, ; of Test V-3 and Tests J-25A and 0-27B provides 
some insight into the surface temperature decrease due to vaporization. These 
three tests were performed at approximately the same apparent flowrate, but 
the bulk liquid temperatures were different. 

An estimate as to the surface temperature cooling effect can be inferred by 
interpolating between the appropirate curves shown in Figure 4-4. For example, 
the Test V-3 breakdown point lies between the expected breakdown point for 
(1) the curve corresponding to breakdown with a 16 inch water column bubble 
point in isopropyl alcohol, with LHg properties for both the surface and bulk 
corresponding to 42.5°R, and (2) the curve for breakdown assuming a 16 inch 
water column bubble point and the viscosity and density of the liquid hydrogen 
evaluated at 42.5°R, but the surface tension evaluated at a surface tempera- 
ture of 37.5°R. Linear interpolation gives an expected temperature at the 
surface of 39.2°R, or an apparent cooling due to evaporation of 3.3°R. Even 
if the flow through the LHAD is assumed to be zero (i.e., all outflow is 
through PV2) then the apparent cooling is estimated to be 2.6°. (In this 
case, it is assumed that V-3 breakdown was due to hydrostatic head alone, 
and Fig^ e 4-4 shows that the expected breakdown would be 5.8 inches above 
the IOU bottom not the observed 4.2 inches.) 

In the Test J cases, it is difficult to use the same limiting flowrate 
reasoning, since the assumption of all outflow taking place through PV2 
(i.e., breakdown assumed to be due to hydrostatic head alone) places the 
appropriate Test J data so close to the expected breakdown level of approxi- 
mately 2.5 inches above the IDU bottom that the deviation falls within experi- 
mental accuracy. This deviation is further obscured by the inaccuracy of the 
liquid level calibration curve very near the bottom of the IOU. However, it 
is unlikely that PV1 failed to open at all, and PV2 failed to close completely, 
in cases 0-11 , J-19, J-25A, and J-27B, for example. Since even a decrease by 
’ factor of two in the actual flowrate through the LHAD for the above cases 
still demonstrates an improvement in the breakdown height, the occurrence of 
a surface cooling effect significantly affecting the bubble point in the cold 
helium pressurization cases of Test J as well as Test V is plausible. 
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These results indicate, but by no means prove, that with cold helium 
pressurization an improvement in screen retention performance can be expected 
with non-equilibrium conditions (i.e., outflow). This would be fortuitous 
from an acquisition design standpoint since the worst case conditions are 
normally encountered when outflow losses are combined with the hydrostatic 
level. Additional non-equilibrium outflow cases with rapid outflow rates are 
required to further assess this potential free surface temperature decrease 
and resultant retention performance improvement. Further, this type of testing 
should be performed under "bench test" conditions, for which higher accuracies 
and better control can be obrained, as opposed to large scale, hardware 
development testing, such as with the IDU. 

One other comparison of the Test V and related Test 0 data is given in Fig- 
ure 4-19, which shows the breakdown level for the two bulk temperatures. As 
expected, all of the related case data shows a tendency for breakdown to 
occur more readily with warmer bulk liquid. The question raised by the closer 
examination of the results is whether or not vaporization induced cooling can 
significantly improve retention performance. 



36 37 38 3° 40 41 42 

PULK LIQUID TEMPERATURE, ^ 

Figure 4-19 OBSERVED RETENTION PERFORMANCE DECREASE WITH TEMPERATURE 
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4.23 TEST W - VACUUM VENT/REFILL 

The objective of this test was to demonstrate start- tank refill following a 
vacuum vent. This test had been performed successfully In MDAC IRAD tests 
of the 10-gallon start tank. In principle, the start tank Is first emptied 
of liquid by flowing back Into the main tank. The start tank Is then vented 
to vacuum, the vent valve closed, and the refill valve opened to allow Inflow. 
Complete refill was demonstrated In the MDAC IRAD tests. However, due to the 
Parker valve malfunctions, the IDU could not be evacuated with liquid hydrogen 
In the main tank, and therefore, this test was not attempted. 
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Section 5 
CONCLUSIONS 

The design, fabrication, assembly, and test of the Interface Demonstration 
Unit has demonstrated that the integrated start tank Is a viable candidate 
for the solution of certain crucial thermodynamic and hydrodynamic problems 
associated with multiple restart cryogenic vehicles. In addition, the I r "i 
itself serve r as a flexible test bed, capable of extensive modification and/or 
easy replacement of components, subassemblies, and instrumentation. As such, 
it should be a useful tool in the contiruing development of reliable, light- 
weight cryogenic propellent acquisition devices. 

Specific objectives, results, and conclusions associated with the development 
and test of the IDU are summarized in Table 5-1. Recommendations for the 
further improvement of the IDU, and its use in extending the technology of 
cryogenic acquisition systems, are delineated in Table 5-2. 
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tee, cone screen for ieedline vapor 
control, screen baffle, insulated 
pressurant lines, and refill diffuser. 
System is assembled with tiared tube 
fittings to facilitate modification. 
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*Wall -mounted thermo- Numerical program with steady-state Output results obtained for insulated 
dynamic vent system conditions met by determining coolant and noninsulated tank, with free con- 
analysis outflow conditions so that surface vection and/or condensation. Program 

integral of wall heat flux is zero. includes pressure drop equations sc 
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corresponding coolant temperatures 
are close to required conditions for 
IDU. 
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Modify ball valve pressurization system for dual mode pressurization 
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Appendix A 

THERMAL ANALYSIS OF THE INTERFACE DEMONSTRATION UNIT (IOU) 
THERMODYNAMIC VENT SYSTEM (TVS) 

The derivation of the governing equations for an LH 2 tank equipped with 
thermodynamic vent system (TVS) cooling coils is given below. The tank is 
cylindrical with a helical TVS coil and optional external insulation. Heat 
is transferred to the tank by convection, condensation, and/or radiation. 
Steady-state conditions are assumed. Equations are first obtained for the case 
of uniform heat transfer coefficient between the heat exchange fluid and the 
tube wall over the tube length. The analysis is then extended to the case of 
two regions of different tube wall heat transfer coefficients in Appendix B. 

A.l DERIVATION OF GOVERNING EQUATIONS 

For steady-state one-dimensional heat flow with the finned tube model shown in 
Figure A-l, the heat balance is (see Figure A-l): 

U-Ollx.d** <«r + VVU (A -’> 

The conductive heat flux out is: 

Q o , x = K T A 3x < a ' 2 > 

The conductive heat flux in is: 

Q iU + dx " *T A 3x + & (K T A ax J dx (A ‘ 3) 

The external net radiative heat flux out is assumed to be constant, independent 
of the wall temperature, since the wall temperatures are of the order of 40° 
to 50°R. The radiative heat flux is, therefore, decoupled from the heat trans- 
fer analysis for the tank, which greatly simplifies this problem. The radiative 
heat flux is given by 

Q r ld x “ q r L dx (A ‘ 4) 
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The external heat flux reaching the metal wall due to convection or 
condensation* and Including the effect of the thermal resistance of the 
Insulation Is 

QJ * U L dx (T - T), (A-5) 

e| dx e 


where U, the overall heat transfer coefficient is given by 


1. J. _L 

U-h e K, 


(A-6) 


(If condensation occurs, radiation heat flux will be neglected.) 


Allowing for the presence of liquid and vapor in the IDU tank requires that 
the heat balance be obtained for these two separate regions. The TVS tube 
length adjacent to the liquid is LH^/Hj and for the vapor, the length is 
L(1 - H|_/Hy) , where H L is the height of the liquid in the tank. The tank 
height is Hy. 

The internal convective heat transfer from the liquid to the wall is 


Q.I = h. L dx (T„ - T) J" , (A-7) 

f 'dx,L f L H L H T 

where it is assumed that the bulk liquid temperature T^ outside the thermal 
boundary layer is constant in the liquid; that is, there is no thermal strati- 
fication in the liquid. This assumption is necessary for this one-dimensional 
analysis. 

The internal convective or condensation heat transfer from the vapor to the 
wall is 


Q f I * h f L dx (T h - T) (1 - rk) ( A- 7a ) 

f 'dx,V f V H V "t 
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The governing differential equation for the liquid region is obtained by 
substituting Equations (A-2) to (A-5) and (A-7a) Into Equation (A-l), assuming 
the tube length is L H^/Hy, which gives 


L H. h 2 t H, H. L H, 

\ 4 ( \ • T) - u (T « ‘ T) 


(A-8) 


or, upon simplification and rearranging. 


d 2 T 

3? 


(h f + U) 

R^E ^ T a 


T) 


(A-9) 


where 


T a » 
a L 


q r + \ \ + UT e 

h, + U 
f L 


(A— 10) 


Similarly, the differential equation for the vapor region is 


d 2 ! 

5? 


(h. ♦ U) 

V 

Kyt 


(T a -T) 
a v 


(A-ll) 


where 


V * \ \ + UT e 

BTTTJ 

T v 


(A— 12) 


The boundary conditions for Equation (A-9) are 


T = T r at X * 0 
l L 


£ *°»tX-0 o /2 


( A- 1 3a ) 
(A- 13b; 
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and for Equation (A-ll), 


T s T. at X ■ 0 
c v 

(A-14a) 

* 0 at X ■ D 0 /2 

(A- 14b) 

Nondimensionalizing with 


V T 

L v 

( A-l 5a) 

and 


r - 

5 ^ 

( A-l 5b) 

L-quation (A-9) assumes the classic form 


d 2e L 2 

— jr - \ \ 

dr 

( A- 16) 

where 


r ( V ♦ ui d o 2 

\-J ?EfT- 

(A-l 7) 

and similarly, fquation (A-ll) becomes 


> 

CD 

CSJ 

•crHf 

(A-l 8) 


138 



where 


0 


v 



- T 



and 


(h f + U) D 2 

V 


The solution to Equation (A-16) for the IDU liquid region is 

cosh N l (1 - C) 
e L = cosh 


(A-19) 


(A- 2C) 


(A-21) 


or 


cosh N, (1 - C) 

\ " ^a L ’ ^ T a L ~ V cosOjj 

The solution to Equation (A-18) for the IDU vapor region is 

cosh N v (1 - 0 
8 v ■ coio; 


( A-22) 


(A-23) 


or 





cosh N v (1 - 0 
cosh M v 


(A-24) 
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The heat flux at X * 0 Is assumed to be transferred to the fluid In the TVS 
tube. The heat entering the tube from one side Is given by, for the liquid 
region. 



- ^ 


LtH L dT 

“HfaST 


o 


Using Equation (A-22), 



* 0 



+ “i v’(\ - t c l > \ 


(A-25) 


Similarly, the heat flux entering the tube from one side is given by, for the 
vapor region, using Equation (A-24), 


q ,| x . o = - L 0 - - T C V > / <\ + U > V’ t>nh N v < A - 26 > 

The total heat transfer removed by the TVS flow, Q, must be equal to twice the 
heat transfer entering one side of the TVS coil, or 

Qt s 2 Q. I +2QI , (A-27) 

T L 'x = 0 v, x = 0 


since heat is transferred to the tube from both sides. 


For long-term storage, it may be necessary that steady-state conditions be main- 
tained. In this case, the net heat entering the tank must be zero. If the heat 
exchanger covers the tank wall, and no extraneous heat flows into the tank, steady- 
state conditions can be maintained by requiring that the surface integral of heat 
flux between the wall and the fluid contained in the IDU, Qwall fluid’ is zero: 

H. 


f V 2 


Q | wall 
fluid 


+ 2 


L iU (t h, 

1 i 


- T ) dx + 2 
L Jr 


,v 2 


(A-28) 


L o - i^ )(T H 2 - V dX 
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Substituting Into Equation (A-28) for T L and Ty from Equations (A-22) and (A-24), 
respectively, gives upon Integrating, 


fluid 


f L L ^[ 


0 - rr) 


DJ2 


D (T a, * T c, } Slnh N. 

T < t H 2 ’ \> + cosh N l N ( . JT) 

L 0 


| + 2h f ^ L (1 - ^ -f (T„ 2 - 


V 


(T a.. " T c> Slnh H, 


+ "v y 

cosh N 


v N v (^) 


1 


(A-29) 


Simplifying Equation (A-29) 

H, 0, 

= + 2 h f L 

fluid 


'Ll! = t2h 


. H L °0 ( 
f L L ^T| 


(T o ■ T c 1 

< T H 2 - \> + -IT— tanh N 1 


H. K 

+ 2 h f L (1 - jh T 


{ 


< T a • T c ) ) 

(T H 2 * T a v > * “V M * 


(A- 30) 


However, If extraneous sources of heat flux into the tank are present, then 
the heat exchanger must flow an additional amount of coolant in order to main- 
tain steady state conditions. In this case, the extraneous heat transfer, Q ext 
into the IDU and the wall heat transfer out of the IDU must be equal, giving. 


Q l 


wall ’ ^ext = ® 
fluid 


(A- 31) 


The above Equation (A-31) gives a relationship for T Cl and T Cv which must be 
satisfied in order to maintain steady-state conditions. 
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A.: IDU FILLED WITH LIQUID 

For the IDU/TVS test, the IDU will be full of LHg. Therefore, this case Is 
treated separately, and Is the condition programmed. The appropriate steady- 
state condition for the IDU filled with liquid Is solved below: 

Since 


^wall " ^ext “ 


Equation (A-29) gives 


L D 


+ 2 h. 


° ( t h - T a > + 
H 2 a L 


(T - T ) tanh N, 

d. C. L 


‘L L 

I 


" Q ext * 0 


(A-32) 


T a L - \ 


W n l 


( T H, ~ T a, ^ tanh N," + TTT 

2 L L 2 h f tanh N L 


(A-33) 


The total heat transfer removed by the TVS is given by Equation (A-27) is 


Or 



(A-34) 


Substituting Equation (A-33) into Equation (A-25) and then into Equation (A-34) 
gives 


«T 


«T 


- - 2 L /(h fi ♦ U) K,t - (T^ - y N l 

■ - 2 "L L v/<\ * U > Vf* \ ■ \> + 


+ 



1 


*ext 

"~n r 


2 h. 


* 


«T ' + 2 L < h f L * U > T {' T H 2 ’ \> 
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Since Qj Is the heat removed by the TVS, the minimum mass flow rate of liquid 
In the TVS Is determined by 

|Q t | » m L A h LV 

where m^ Is the Initial fraction of liquid mass flow rate resulting from the 
expansion through the Inlet orifice. 

It Is assumed that m^, the mass flowrate of the liquid In the TVS, Is known 
and that vaporization transforms all liquid in the TVS to vapor but no super- 
heating of the vapor occurs. Knowing m L for the steady state gives Qy for the 
experiment, which allows an overall heat transfer coefficient on the Inner sur- 
face of the IDU to be calculated, if extraneous sources of heat can be deter- 
mined, by using Equation (A- 35). 

The heat flux transferred to the fluid In the TVS assuming the IDU is full of 
liquid Is given In terms of the saturation temperature, T H ^ X of the two-phase 
fluid In the TVS and the tube Internal heat transfer coefficient as: 


- "D h. L (T - T hex ) , (A- 37) 

| liquid tul)e 1 C L HEX 

where it Is assumed that T^ x is a constant, no superheat of the TVS vapor is 

present, there Is no significant pressure drop along the TVS line, and h^ is 

a constant. 


For given heat transfer coefficients, the IDU and external temperatures deter- 
mine Qj and T c ^; the TVS tube length can be evaluated knowing the heat exchanger 
temperature, T^, or vice versa. 

Since L is known for the IDU, the necessary T^ x is determined as 

Q T 

T ucv = T - -= 4 — _ , where T. 1 c given by Equation (A-33) (A- 38) 

Htx c L TrD tu 5 e n 1 l c L 

(The saturation temperature of the TVS fluid must be lowered as Qj Increases.) 
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The above expersslon Is used to establish the TVS pressure required to obtain 
the coolant T^ which corresponds to steady-state r.r? Mtions In the tube. Or, 
knowing Qy, T c ^, and T H £ X , the associated average h^ can be determined. 

From the standpoint of controlling the TVS so as to maintain steady-state con- 
ditions, It Is seen that the mass flowrate must be set such that the total heat 
flux, Qy, can be removed and such that the TVS fluid temperature, T^, deter- 
mined by the tank and TVS pressures must be low enough to allow Qy heat to be 
transferred to the fluid, given h^ and L. However, It should be noted that In 
practice a greater mass flowrate may be used in the TVS than is required for 
the IDU, with the remainder of the liquid used for cooling downstream equipment. 


As a check, the heat transfer to the tank wall from the surrounding region plus 
the extraneous heat is determined, since this must be identically equal In 
magnitude and opposite in size to Qy, the heat transfer removed from the IDU, 
as obtained in Equation (A- 35). Heat transferred to the tank *' given by 


, f v 2 

2 J UL (T - T) dx t Q 


■ext 


ext 
(T a - T. ) 


C D n D u a, “ c, ) 

2 UL I T (T e ' T a l ) * 1 \ — ~ tanh \ j + 


(A- 39) 


(A-40) 


Using Equation (A- 33), 


«T ’ 2 UL {T<Vy T<W *rKf T H e «‘ 


(A-41) 


‘ «■ D o < T e - t h I + + HT-) «ext 

2 f L 

From Equation (A- 35), using the expression for T 


Q_ - + 


D C ^f ^H + 0 

2 l( V u> + K - — sfH] 


(A-42) 
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which gives, for the heat transfer from the IDU, 

«T * - U L D o < T e ’ T H > ' 0 * T,T-) «ext < A ‘ 43 > 

* t l 

Since Qy determined by the heat transfer from the liquid in the IDU Is identical 
In magnitude and opposite in size with that based on the heat transfer from the 
surrounding region, the analysis is consistent. 

Heat transfer coefficients are given in Table A-l based on the compilation of 
Reference A-l. .'t this point, there is no applicable information on the flow 
pattern of liquid or gas with a vertical wail subjected to cyclical temperature 
variations, nor is it known what type of condensate flow occurs. Standard 
equations for the heat transfer equation are used in the analysis which can 
be modified where appropriate following additional TVS tests. The assumptions 
used in selecting the expressions for the heat transfer coefficients are 
summarized below. 

A. 2.1 Condensation - Outer 

Film condensation is assumed and the Nusselt film theory expression is used. 
Condensate is assumed to flow down a characteristic distance, D Q , equal tc the 
tube spacing. The outer condensate temperature is equal to the saturation 
temperature of the surrounding vapor (GH 2 ) and the wall temperauure is approxi- 
mated as an average temperature, (T. + T )/2. If (T + T m 1/2 is greater 
than the saturation temperature, condensation is assumed not to occur. 

A. 2. 2 Free Convection - Gas, Outer 

The free convection he^t. transfer outside the IDU is assumed to be due to con- 
vection cells such as are shown in Table A-l. The characteristic length, x, 
is D q /2 In the heat transfer coefficient expression. The temperature difference 
between the area of width D./2 centered on the tube is approximated as T - T 
which is high, since this area has a temperature gradient and the average wall 
temperature is less than T . 

Between the coils the wall temperature is assumed to be high enough such that 
the temperature difference is T ma „ - T_, corresponding to free convection with 
a warm wall and cold gas. The direction of rotation of the associated cell 
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Table A-l 

IDU HEAT TRANSFER COEFFICIENTS 



i 





pattern is opposite to that for the cold wall /warm gas combination. The 
equations for free convection are valid for Rayleigh numbers between 103 and 10^. 

The calculated heat transfer coefficients are averaged to obtain an overall 
heat transfer coefficient. Thus, 


h « h e>wg..' f .. h e t .c g , 
e 2 


The temperature differences assumed are not expected to have a major effect on 
the values of the heat transfer coefficients, since they are taken to the 1/4 
power in the equation. Moreover, the most significant approximation is the use 
of free convection boundary layer expressions for a complicated ablation flow 
pattern, for which there is no known analysis. 


In some cases, T ma „ will be less than T Q , and therefore, the entire IDU wall is 

niaX 6 

cold with respect to the surrounding fluid. In this case, the characteristic 


length is x = Hj, and the temperature difference assumed in determining h 0 wg 

72. 

L 


< t c l + T max )/ I 2 - 


A. 2. 3 Free Convection - Liquid Inner 

The alternating cellular flow pattern is assumed to exist inside the IDU with 
a characteristic length of X = D Q /2. The temperature difference between the 
cold wall and warm liquid is assumed to T^ - T c ^, which is high. 


The corresponding temperature difference between the warm wall and cold liquid 
is T„„ , - 7 U . Again, these equations are assumed to be valid for R?yl ‘ ! gh 
numbers between 10^ and 1(P* and the overall heat transfer coefficient on the 
inside of the IDU cylindrical wall is the average of the two values calculated: 


h f • (h f + h- )/2. 
T l a.wl T L,cl 
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Table A-2 (Page 1 of 2) 
EXTRANEOUS HEAT TRANSFER - TOP OF I DU TANK 





INSULATION 


IOU TANK TOP 


LIQUID HYDROGEN 


Extraneous Heat Flux Relation Based on Total Thermal Resistance 
Q EXT = U T0 P A c ^ T e ‘ V 

_j_ . i + h. + ^ + 1 

^TOP ^.external *4 K w *4, internal 


Heat Flux Relations for Individual Thermal Resistances 

Free convection heat transfer from external gas to insulated top of tank 
at temperature Tjy j: 

Q EXT = h T, external A c (T e ' T TT,I^ 

Free convection heat transfer fiom external gas to uninsulated top of tank 
at temperature Tjj 

Q EXT = h T, external A c (T e " T TT,M^ 

Conductive heat transfer through insulation 

K I A r 

Q EXT = Tj - (T TT,I ‘ T TT,M* 




) 


} 


Table A-2 (Page 2 of 2) 
EXTRANEOUS HEAT TRANSFER - TOP OF IDU TANK 


Conductive heat transfer through metal tank top 


<EXT 


, Vc (T _ T , 

' TT,M 'TB,M j 


Free convection heat transfer from bottom wall of metal tank top to IDU 
hydrogen vapor 


*EXT 


h T, internal A c * T TB,M ‘ 


h T, external ” K g 


p g p g c p,g 9 (T e ~ t tt,m ) 

Vg 


n 1/3 


h T . . , = 0.08 K 

T, internal g 


p g p g c p,g 9 ^ t tb,m ' Th 2* 


g g 


1/3 
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A.2.4 Extraneous Heat Transfer 

The extraneous heat transfer is not known nor is it practical with the present 
IDU instrumentation to determine it. However, an estimate can be made based 
on the assumption that the major portion of this heat flux enters the IDU 
through the top of the cover. Free convection heat transfer coefficients are 
assumed and insulation may be included, as shown in Table A-2. It is assumed 
that the IDU is not completely full and, therefore, an ullage gas layer exists 
between the liquid surface and the top of the cover. 

A.2.5 Forced Convection — TVS Tube (Mist Flow Regine) 

At relatively high heat fluxes, a mist flow regime is possible since liquid 
droplets no longer "wet" the tube wall, due to the momentum ot the rapidly 
evaporating vapor between the liquid droplet and heat surface supporting drop- 
let. This is referred io as stable film boiling, and occurs at or above the 
Leidenfrost point. Measurements of heat transfer coefficient for this regime 
are close to the corresponding values for the vapor alone. 


As a first approximation, the Dittus-Boelter equation for forced convection 
heat transfer coefficient is used, which gives 


.023 k„ 0.8 0.4 

h 1 ■ -TT-^v' <%> 


where 



4 IT) 




2. 6(. 00202) 
— .0074 


0.71 at T * 31°R 
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However, In normal operations with the IDU/TVS, the heat flux is relatively 
low. Assuming a maximum value of 600 Dtu/hr, distributed over the surface of 
the 70 feet of 0. 25-Inch O.D. tube, the heat flux is 130 Btu/ft 2 hr. The peak 
nucleate boiling heat flux for liquid hydrogen is reported to be 32,000 Btu/ft 2 hr 
In Reference A-2. Since the heat flux for the TVS is far below that required 
for peak nucleate boiling, it Is probable that the heat transfer coefficient 
will correspond to that of the two-phase flow in the initial region. The heat 
transfer calculations in this regime are more involved and require that two heat 
transfer regions be treated. The associated analysis is treated in Appendix B. 

\ 

IDU/TVS HEAT TRANSFER ETFECTS ON LIQUID FLOWING FROM IDU TO ORIFICE MANIFOLD 
When the IDU is completely surrounded by warm ullage gas in the 105" APS tank, 
heat transfer to the liquid hydrogen coolant flowing from the IDU/TVS feedline 
outlet to the orifice manifold through the initial inlet section can cause 
boiling before the liquid reaches the orifices. The severity of this problem 
and the importance of TVS line insulation is determined below. In most appli- 
cations, however, the lower portion of the IDU will be submerged in liquid, 
and therefore, boiling in the initial inlet section will not occur. 

The tube leading from the start tank feedline to the orifice manifold is taken 
to be about 15 inches long (including end fittings). The tube is 1/4 inch 
stainless steel with a wall thickness of 0.020 inch. 

Assume hydrogen properties at 10C°R and 50 psi . 

Density, p= 0.095 lb/ft 3 (data range from 0.088 to 0.102) 

Specific heat, Cp = 2.61 Btu/Lb°R 
Viscosity, p = 1.49 x 10"^ lb/in-sec 
= 1.79 x 10"® Ib/ft-sec 

Thermal Conductivity, k = 0.55 x 10"® Btu/in-sec°R 

= 6.6 x 10"® Btu/ft-sec°R 
Thermal Expansion Coefficient, 0 = 0.022 (° R) - "* 

For a horizontal cylinder in laminar flow, the Nusselt number for free 
convection is 

% = T = 0,53 (R a )1/4 ’ 1q4 < R a <1q8 
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and in turbulent flow 

N = 0.13 (R.) 1/3 , 10 4 < R <10 12 < A ' 45) 

U j a a 

Free convection heat transfer will actually be somewhat lower because of 
the proximity of the cold IDU bottom; it shall be assumed to be about 

\ = °- 4 (R a )V4; \ = °- 1 (R a )V3 (A-46) 

where 

R a = 9 ~ ^ F~ = I * 16 x ] ° 3 at (A-47) 

Since temperature differences of the order of 100 R are assumed, the flow 
is laminar and the Nussult number is 

N u = 2.34 (AT ) 1/4 (A-48) 

The free convection heat transfer rate is 
Q = hAAT = | Nu (ttDLAT) 

Q = 0.218 AT 5/4 Btu/hr (A-49) 

Let the initial state of the liquid hydrogen be saturated at 14.7 psia, 
but pressurized at 30 psia. The liquid bulk temperature is T = 36 . 5°f 

The initial enthalpy, H Q , is -108 Btu/lb, the saturated liquid enthalpy at 
30 psia, H 30 £, is -95 Btu/lb, and the saturated vapor enthalpy at 30 psia, 
H 30 y » is +90 Btu/lb. At steady state, the liquid enthalpy reaching the 
orifices is 

H = H q + Q/W 
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where W is the mass flow rate, Ibs/hr. The mass flow rate must be W 2 Q/13 
in order for pure liquid to reach the orifice for the above case. 

The enthalpy of saturated liquid at 14.7 psi is -108 Btu/lb, and the enthalpy 
of saturated vapor is +82 Btu/lb. The heat of vaporization is, therefore, 

190 Btu/lb, and the quality of liquid reaching the orifices is 


X 



1.09 x 10” 3 



(A- 50) 


if the liquid is stored in the tank at the saturation pressure of 14.7 psi. 


In general, the quality of the liquid reaching the orifice will be, for the 
tank pressurized to 30 psi a. 


X 


0 - W < H 30.* ~ H o> 

Wh. 


( A- 51 ) 


where 


h v = H 30 ,v “ H 30,t 

In order for pure liquid to reach the orifices, X must be zero and, therefore, 
high flowrates and/or high tank pressures must be applied if the heat flux 
into the line is significant. For example, Q = 29 Btu/hr with A T - 50°R, 
which requires a flowrate in excess of 2.2 Ib/hr with a tank pressure of 
30 psia, if pure liquid is to be supplied to the orifice. 

Since it is necessary to be able to supply lower coolant flowrates, the need 
for insulation on the TVS inlet line up to and including the viscojet flow 
control orifices is apparent. In addition, condensation heat transfer rates 
are higher than the free convection rates used in this analysis, further 
indicating the necessity for insulation in the event tne I DU tank bottom is 
exposed to hydrogen ullage. 
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Appendix B 

I DU/TVS ANALYSIS - TWO REGION MODEL 

The analysis of the IDU/TVS presented in Appendix A assumes that the heat transfer 
coefficient between the heat exchanger fluid and the tube wall is uniform over 
the tube length, which corresponds to the case of oversupplying LH 2 to the TVS, 
in order to use the liquid remaining at the TVS exit for cooling of components, 
turbopumps, etc. If the flow rate is such that the fluid quantity exceeds a 
critical value, X c , then there will be two regions in this two- phase tube flow 
with different mechanisms of heat transfer and sigmificantly different transfer 
coefficients. An analysis of this two-region heat transfer model performed by 
D. W. Kendle is developed below. 

The model is shown in Figure B-l. For the purposes of analysis, the tank sur- 
face and the heat exchanger tube length attached to this surface are divided 
into two regions as shown. Coolant fluid flows first through region 1 and then 
region 2; mass flow m and coolant temperature T hex are the same for both regions. 
Two-phase fluid enters region 1 with quality X-j and leaves at Xc (the value of 
X c determines the boundary between the two regions). Entering with Xq, the fluid 
exits region 2 with a quality of Xp. Associated heat transfer coefficients are 
given in Appendix C. 

Since the coolant temperature is the same throughout, the region with the higher 
heat transfer coefficient, assumed to be region 1, will have the lower tank wall 
temperature. In order for the net heat input to the tank contents to be zero 
when there are significant wall temperature differences, heat will pass through 
the warmer wall area (region 2) into the fluid and be removed from the fluid in 
the cooler wall area (region 1). The heat transfer rate between the two regions 
is Q(j. The extraneous heat input, Q ext > enters the tank contents by means other 
than through the cooled wall; Q ext is removed from the fluid by the cooled wall 
in region 1 . 
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FIGURE B-l. - IDU/TVS COMFIGURATION 







With this model, the thermal behavior of each region can be described by the 
same equations derived In Appendix A for a single region. The additional con- 
ditions described above link the two regions and enable a solution for the 
Individual tube lengths, li and L 2 , In each region. 

The fin root (tube wall) temperatures are given by Equation (A-33) In Appendix A 
for th" two regions of this model, they are written 


\ * < t h - \> 


u. 


, . lq .*t * V »1 

tanh N. ** L.h* D tanh H, 

X JL i ^ O X 


N, 


Si N 2 


?C 2 ' % + (Tu " tanh N 2 + " L 2 h f 2 D o tanh 


In the above equations, 


% + h fi t h * u i t e 

* U 1 
X 


(B-l) 


(B-2) 


(B-3) 


from Equation (A-10) of Appendix A and 


N 


i 



(B-4) 


from Equation ( A- 1 7 ) . c Appendix A, where subscript 1 can be either 1 or 2. 
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The Q ext term of Appendix A Is replaced In Equations (B-l) and (B-2) by 
(Q ex t + Q(j) for region 1 and (-Q d ) for region 2. Making this same substitution 
In Equation (A-43) of Appendix A gives 




(B-5) 


\ ■ Wo < T E - V - (* * £ ) «d 


(B-6) 


The heat wxchanger fluid temperature, given by Equation (A-38) In Appendix A, 
Is the same for both regions: 




^r. 


T hex “ % ' 1 D tube L x ^ 


n ^tube ^2 


(B-7) 


Substituting Equations (B-l), (B-2), (B-5), and (B-6) into Equation (B-7) gives 


* 


•.(Hr*)- v,(S) 


(B-8) 


where 


U* D_ ( Tp - T h ) 


B. 


rn + /t - T \ 1 l o E a 

A. ( H A. ) tanh N, 7T D h_ 

* \ 1/ i tube T. 

C • « 


N i 

D q tanh 


if D 


tube 


V 


(B-9) 


(B-10) 
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where subscript 1 Is again 1 or 2. Solving Equation (B-8) for the rate of heat 
transfer Internally from zone 2 to zone 1 gives 


il, - A 0 - Q . 

l 2 L. ext 

(B-ll) 

' (H) 

Adapting Equation (A- 36) from Appendix A. the heat capacity of the fluid in 
each section Is related to the heat load by 

m Ah y (X c - Y z ) . % 

(B-12) 

m Ah y (Xy - X c ) = {Jj 

(B-13) 


where m Is the mass flowrate of heat exchanger fluid and Ah y Is the heat of 
vaporization of thf .’uid. Eliminating mAh y between these two equations and 
substituting Equations (B-5) and (B-6) gives 



(B-14 


where 

Cl ’ <*, - X c ) U 1 D o (t e - V 

C 2 - (X c - V U 2 D, (T e - T„> 
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D 1 ■ B 1 ( *c - V ( J ♦ s; ) 
d 2 • b 2 <s*v(i + ^) ««t 

The total tube length L Is 

L| + Lj 5 L (B-15) 

which is .ubstituted into Equation (B-14) to eliminate 1^ giving a quadratic 
in L| : 

AL 2 + BI^ + c = 0 (B-16) 

A = C 3 - (C (B 2 - B 1 ) 

B = (B 2 - B 1 ) - B^ + C 2 ) - (D x 4 D 2 ) - C 3 L 

C = C^L 2 + D X L 

Since the coefficients of this quadratic are dependent on the temperature of 
the system which are in turn a function of L^, the solution must be obtained 
iteratively. Values for the system temperatures are assumed initially and the 
coefficients are calculated, giving an approximate value of L^. The system 
temperatures are calculated based on this and the iteration is repeated 
until things settle down. 
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Appendix C 

INTERNAL HEAT TRANSFER EQUATIONS FOR 
INTERFACE DEMONSTRATION UNIT THERMODYNAMIC VENT SYSTEM 

C.l ANNULAR FLOW 

Annular flow in the TVS is assumed from the viscojets to a distance along the 
tube at which transition to liquid deficient flow occurs. The model is based 
on the work of Chen (Reference D-l) who used Forster and Zuber's mi croconvecti ve 
heat transfer relation for boiling and the Dittus-Boeltor macroconvective heat 
transfer relation for forced convection, and obtains: 

h = h MIC + h MAC ( C_1 ) 


or 


0.79 0.45 0.49 g 0.25 iT 0.24 ip 0.75 

h = S (0.00122) 0.5"U.29 0.24 ff.24 

0 \ h fg °v 

♦ F (0.023) (Re)?‘ 8 (Pr)°* 4 i (C-2) 

L L U 

where S and F are empirically determined dimensionless functions which ..llow 
for variations in the boiling and forced convection components, respectively. 
The value of S and F are given in Figures C-l and C-2. 

C.2 VAPOR FLOW (LIQUID DEFICIENT REGIME) 

At an initial quality, X c , transition from annular flow to vapor flow occurs, 
and the Dittus-Boelter equation for forced convection in a tube is used: 



0.023 ( 



U 0,8 
-) 



(C-3) 
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Appendix D 

I DU/TVS TWO PHASE FLOW PRESSURE DROP 

An approximate solution to the pressure drop of the TVS coolant flow Is 
obtained using the homogeneous (Fog Flow) flow model (Reference D-l). The 
model assumes equal linear velocities of vapor and liquid, thermodynamic 
equilibrium between the phases, and use of a single-phase friction factor which 
is applicable to the two-phase flow. The pressure drop through the Viscojets 
is calculated separately. 

D-l. TVS LINE PRESSURE DROP 

The total pressure gradient, composed of frictional and dynamic losses, is 
given by (Reference D-l): 



(Z = length of tube, D^ is inner diameter, f is faction factor, 
m is total mass flowrate, X is quality, P^is liquid density, 

P v is vapcr density, P is pressure) 
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PRESSURE DROP PSIA 








The total pressure drop is obtained by integrating Equations (D-2) and (D-3) 
with the quality, X, varying with the tube length, Z. Assuming X increases 
linearly with Z, which implies a constant heat flow, obtain 


x 


Y 4. (Y y ) z is final quality, X. is 

i 1 f i' L initial quality) 


Substituting Equation (0-4) into Equation (D-2) gives, for the frictional 
pressure gradient 



Integrating Equation (D-5) gives, for the frictional pressure drop. 


AP, 


FRICTION 


BL 


(A + B X £ ) L 
c (x f - X.) ' 


In 


1 + 


C (X f - X £ ) 
1 + C X. 

l 


{ C (X. - X.) - (1+ CX.) In 

C 2 (X f - Xi ) ( f 1 1 


1 + 


C (x f - X.) 
1 + C X. 


Sirilarly, the dynamic pressure gradient is 


dP I 
dZ I 

DYNAMIC 




(D-4) 


(D-5) 


(D-6) 


(D-7) 
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Integrating Equation (D-7) gives 

c (x, - xj 


AP s — In 

DYNAMIC C 


r c (x. - x. j 

i 1 . 

1+ CXj 


(D-8) 


By definition. 


E f„ 


2 p i D i 


M 0 . ) 


B •= A 


C = - 


ft-) 


E 

P P 
v 


E — 

\* D j / 


Following Owens, (References D-2), the friction factor of the single-phase liquid 
flow is used for the two-phase friction factor. 


'D. 

l 


f = 


factor is 

given by: 


<3000 , 

( Re Di * Di u ) 

(D-9) 

6k 

Re D. 

1 


(D-10) 

>3000 




f » 


0.316 


K) 


ia 


(D-ll) 
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Since the compressibility factor of vapor Is 0 96 and does not vary appreciably 
with pressure for conditions In the TVS tube, the perfect gas law can l*e used 
in the expression C, giving 





assuming m = 5 lb/hr, 5 0.25 in., the expression for C is found to be of the 
order of 0.01. It should be noted that since C«l, the pressure drop equation 
containing the log term can be simplified using 

lim I In (1 + C) = C 
C * 0 


This simplification illustrates the appropriate linear dependence of pressure 

m o 

with m , for the dynamic pressure drop and turbulent frictional pressure drop, 
and the linear dependence of pressure with m for laminar flow. 

The curves of Figure D-l for conditions corresponding to those of the IDU/TVS 
tests show: (1) that the dynamic pressure drop is of the order of 1 percent of 

the frictional pressure drop; and (2) the total pressure drop is of the order of 
a few psi. 

The heat transfer part of the program gives the required TVS tube pressure, 
^TUBE' must e 9 ua ^ the avera 9 e pressure in the tube. 


- ?ima TVS * 

TUBE * 2 

CALCULATED 


(D-l 2) 


The pressure loss is 

ap loss * p ihlet tvs “ p exit tvs 


(D-l 3) 
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Thus, the Inlet pressure, immediately downstream of the viscojet, must br* 
adjusted such that 


AP 

P = p . LOSS 

r INL2T r TUBE 2 

CALCULATED 


(D-14) 


or 


■ EXIT, 


TVS 


AP 

LOSS 

TUBE “ 2 

CALCULATED 


(D-15) 


D.2 PRESSURE DROP THROUGH THE LEE VISCOdETS 

Lee viscojets are used as flow control orifices since their design provides a 
flow area 10 to 25 times that of an equivalent orifice having the same resistance 
(Reference D-3). Therefore, there is less likelihood of clogging. It is 
necessary to determine the effective resistance of the various viscojets used 
in the IDU/TVS, and these resistances vary over a wide range, depending on 
whether liquid, gas, or a two- phase mixture is flowing. 

A previous MDAC study (Reference D-4) obtained data of flowra e vs pressure 
drop through the viscojets. A reexamination of this data has shown that the 
effective resistance for two phase flow is essentially the same as for liquid 
flow alone, for conditions of < 0.1. The resistance of the gas flow, even 
assuming choked flow, is relatively small. Therefore, the analysis given below 
assumes that the two-phase flow pressure drop/flowrate dependence is given by 
that for all liquid .low. An improvement in the pressure drop/ flowrate cor- 
relation is achieved if the two-phase flow density, p jpp» is used, rather than 
the liquid density, as discussed in Section 4.13. 

Information received on the Lee viscojets (Reference D-5) indicates that flow 
resistance is treated in terms of LOHMS, with 


1270 f 

EU 


V 7 


P Pr 


L (LOHMS) 
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where LOHM Is a flow resistance unit corresponding to the flow of 100 gallons/ 
minute of 80°F water with a pressure drop of 25 psl. The factor f was deter- 
mined In Reference D-4 from flow tests with subcooled LHg (such that vaporization 
did not occur, and X * 0). The value obtained was f = 0.823, which will be used 
In our analysis. 

The values of L for the visojets used In the IDU/TVS are given In Figure D-5 
which shows the expected flowrates versus pressure drop. Note that the cor- 
relation of Equation D-16 Is Improved by using P TpF instead of ° L , which allows 
the equation to be extended into the range of low quality two-phase flow. 

The actual flowrates versus pressure drop compared well with predicted values. 
Knowing m vs P for the viscojets is necessary for the proper setting of the 
IDU tank pressure. The pressure at the TVS line inlet is given by Equa- 
tion (D-14) and the corresponding mass flowrate for given heat transfer con- 
ditions is given by a curve such as that of Figure 3-3 of this report. Knowing 
m, the IDU tank pressure must be set such that Pjqu" p inleT supplies the required 

mass flowrate. 

D.3 REFERENCES 

D-l. Tong, L. S. Boiling Heat Transfer and Two Phase Flow, 1965, John Wiley & 
Sons. 

D-2. Owens, W. L. Two Phase Pressure Gradient, pp 363-368, International 
Developments in Heat Transfer, Pt II, ASME, (1961). 

D-3. Lee Viscojet Catalog, 1972. 

D-4. Castle, J. N. and Cady, E. C. Performance Testing of An Integrated Liquid 
Hydrogen Storage, Acquisition, and Vent System MDAC Report No. MDC G3092, 
June, 1972. 

D-5. Personal Communication. Walt Reynolds (Lee Company) to E. C. Cady (MDAC), 
December 31, 1970. 


169 


i 



? 


Appendix E 

IDU/TVS PERFORMANCE COMPUTER CODE DESCRIPTION 

E. 1 PROGRAM STRUCTURE 
E . 1 . 1 Analysis 

The equations included in this program are derived in Appendices A, B, C, and 

D. Some additional features of this analysis as it is interpreted and incor- 
porated into thv code will be briefly described. 

E. l .1 .1 Multiple Tubes 

Appendices A and B apply to a single tube for the vent flow. The program has 
the capability of parallel flow in multiple, eque.l-lenqth tubes. This is 
accomplished by assuming that Equation ( B- 1 6) for is written in terms of 
a single tube length and then multiplying L-j, L 2 and L by the number of tubes 
in the heat transfer equations to provide the correct overall heat balance. 

E.l. 1.2 Condensation Heat Transfer 

The condensation heat transfer equation is used to calculate the external wall 
coefficient when the average wall temperature is less than the saturation tempei - 
ature at the external pressure. When this occurs, the external temperature T g 
in Equations (A-10) and (A-41 ) is replaced by T , the external saturation 
temperature. 

E.l. 1.3 Insulation 

Throughout the program, options are provided for all computations to be per- 
formed with or without external insulation. Insulation may be specified inde- 
pendently for the wall and the tank top. 

E.l. 1.4 Conductivities 

Constant values of tank wall and insulation conductivity are input to the 
program. Separate variables are used within the code f~ tank top, region 1 
wall and region 2 wall. These conductivity values could be read from a tem- 
perature-dependent table or equations, but such temperature-dependent data 
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were not available with sufficient reliability. This change can be made at 
a later date if desired. 

E.l.1.5 Radiation 

The radiation heat transfer analysis was not developed. However, this variable 
is carried in the code and the subroutine is provided for its evaluation to 
simplify the task of adding this feature at a later date if desired. 

E.l . 1 .6 Tank Top 

The equations given in Table A ? for the extraneous heat flux through the tank 
top are so'vpd by direct substitution. Two solution options are included in 
the code, witn and without tank top insulation. 

E.l. 1.7 Region 1 Tube Length 

Equation (B- 16) is a quadratic in L-j . However, it is possible for the coefficient 
of the squared term to be zero so that the solution would not be operable. To 
circumvent this possibility, this coefficient is checked against a minimum 
lim if its magnitude is less than this limit, the equation is solved approxi- 
mately as a linear equation and corrected iteratively for the small squared term. 
If the coefficient is zero, the linear solution will be exact. 

E.l. 1.8 Single Flow Regimp Solution 

The solution for a single flow regime throughout the tube can be obtained by 
setting the transition quality X c equal to the final quality of the flow in 
the tube to eliminate the mist flow regime. The reverse cannot be done to com- 
pletely eliminate the annular flow regime; however, the annular flow regime can 
be made insignificantly small by setting the value at X £ ven close to that of 
the initial flow quality. 

E.l .1 .9 Curve Fits 

The heat transfer coefficient equation in the annular flow regime of the tube 
flow is discussed in Appendix C. The suppression factor S and the Reynolds 
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number factor F used in equation (C-2) are shown as curves in Figure C-l and 
Figure C-2. These curves are represented by equations in the code as follows: 


1 

1 + 0.305 x Ilf 5 


Y 1 ■ Re L F 


1.25 


and 


F « (!.-»* 2.25 Y,) 


0.8 


* 0.9 p 0 0.5 p,. 0.1 

Y 2 <i4t> $ 


(E-l) 


(E-2) 


An average value of the latter quantity is evaluated by integrating the X-term 

from X T to X_ . 

I c 

E. 1.1. 10 Properties 

Tabular data for the fluid properties are grouped in the BLOCK DATA subroutine. 
Other properties are approximated by linear relationships; these are found in 
the subroutines which calculate the heat transfer coefficients at the inner 
wall, outer wall and inside the coolant tube. 

E.1.2 Solution Technique 

The computer code effects an iterative solution to the quadratic equation 
derived in the preceding analysis for the region 1 flow rsgime (annular two-phase 
flow) tube length. Iteration is necessary because the equation contains many 
temperature dependent terms and these temperatures cannot be determined until 
the region 1 tube length is known. 

The iteration is initiated by assigning estimates to the unknown temperatures. 
Based on these values, the coefficients of Equation (B-16) are evaluated and 
the equation solved for . With this approximate solution, the system tem- 
peratures are then calculated and the process repeated. However, a straight- 
forward application of this simple iteration process frequently does not 
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converge due to extreme fluctuations in some of the variables. To overcome 
this difficulty, the iteration proceeds in the manner illustrated by the 
flow chart in Figure E-l. 

After the first evaluation of L^, that value is held constant while the tem- 
peratures and heat transfer coefficients are iterated to achieve a moderate 
degree of convergence. With the temperat ires thus stabilized, a new value of 
is calculated and the temperature iteration is repeated. In this manner, a 
well-behaved convergence of is attained. This two-level iteration is 
achieved by the two nested D0-loops shown in Figure E-l . Since tube lengths 
are evaluated only on the first pass through the inner loop, their calculation 
is shown logically to be a part of the outer loop. 

As approaches it solution value, ti.e iterative variation of the system tem- 
peratures decreases and fewer temperature/Droperties iterations are required. 
Unnecessary computation is avoided by decreasing the count limit on the inner 
loon as the solution proceeds. In the first four passes through the outer 
DJv-loop, the inner D0-loop is executed 6 , 3, 2 and 1 times; thereafter, the 
entire iteration is performed as a single loop. 

Checks for convergence are initiated on the ti.ird pass through the outer loop. 
The difference between the maximum and minimum values of the three most recent 
flowrates must be less than a specified limit to exit the iteration loop. If 
the error check is not satisfied in ten iterations, an error message is printed 
and the variable values at that point are printed in the normal manner. 

A complete program 'isting, a gloss *7 of principal variable names and a 
diagram of subroutine relationships (calling sequence) are given in Appendix F. 
Comment cards are provided in trie listing to describe the function of each 
subroutine and to identify the principal operations within the code. 

E.2 PROGRAM UTILIZATION 
E . 2 . 1 Capabilitie s a nd L i m it ations 

This code calculates the •, ot rate necessary to maintain a r-ro net heat irn t 
into the tank. The capabilities and lir ‘tatiors ’’nherent in the model and 
analysis are adequately described in the derivation of the equations in 
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Appendices A and B, and in the preceding section on analysis. The tank 
configuration is a flat-end cylinder. The product of the number of tubes, 
single- tube length, and tube spacing must equal the cylindrical wall area. . 

The solution is restricted to the case in which the outer wall of the tank sees 
gas only (but it may have a film of condensation), the inner wall sees liquid 
only, and the tank top sees gas on both sides. 

Since this code treats the transfer of heat from the external environment to 
the tank, it will not operate unless the external temperature is greater than 
the temperature of the tank contents. Also, as mentioned in the preceding 
section on analysis, the code does not permit the complete elimination of 
the annular flow regime by setting the transition quality X £ equal to the ini- 
tial quality Xj. If the input data violate either of these t estrictions, the 
code will skip to the next input case. 

E.2.2 Input 

The input variables are listed in Table E-l . Each variable is identified by 
a number which precedes it on the input card. All cards are written in the 
same format (13, E12.7). The cards may be arranged in any order within an 
input case. The end of an input case is designated by a blank card. The end 
of a run is designated by two blank cards. 

Data must be input for all variables in the first case. Exceptions for which 
default values are supplied are given in Table E-2. In the case of the values 
read from the vapor pressure table, either the pressure or temperature must be 
input for the other value to be calculated. 

For subsequent cases, only those variables are input for which new data are 
specified. The remaining variables retain their values from the previous case. 

A listing of input cards for a series of test cases is given in Appendix G. 

£.2.3 Output 

The output format is illustrated in the test case output in Appendix G. The 
listi g is separated into three sections by blank lines. The first section lists 
the input data, the second lists solution data, and the third lists some of the 
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Table E-l 
INPUT VARIABLES 


ID No. 

Variable Name 

Definition and Units 

1 

DO 

Tube spacing (in.) 

2 

EL 

Single tube length (ft) 

3 

Y 

Number of tubes 

4 

DTUBE 

Outside diameter of tube (in.) 

5 

TWTH 

Tube wall thickness (in.) 

6 

T 

Tank wall thickness (in.) 

7 

TI 

Wall insulation thickness (in.) 

8 

ATT 

O 

Tank top area (ft ) 

9 

TTM 

Tank top thickness (in.) 

10 

TTI 

Tank top insulation thickness (in.) 

11 

HTK 

Tank height (ft) 

12 

AKTx 

Tank conductivity (Btu/hr-ft-°R) 

13 

AKIx 

Insulation conductivity (Btu/hr-ft-°R) 

14 

TE 

External temperature (°R) 

15 

TSE 

External saturation temperature (°R) 

16 

P 

External pressure (Psia) 

17 

TH2 

Internal temperature (°R) 

18 

PTT 

Internal pressure (Psia) 

19 

XI 

Initial quality in tube (fraction) 

20 

XC 

Transition quality in tube (fraction) 

21 

XF 

Exit quality in tube (fraction) 

22 

ACC 

Acceleration (g) 


Notes: Each input card is occupied by a single variable. Format for 

all inputs is I3.E12.7. Input No. 12 is assigned to AKT1, AKT2, 
and AKTT; Input No. 13 is assigned to AKI1, AKI2, and AKIT. 
End-of-case is indicated by a blank card; end-of-run by two 
blank cards. 
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Table E-2 

INPUT DATA DEFAULT VALUES 


Variable Name 

Default Value 

Y 

1.0 

AKIx 

0.1 

TE 

TSE 

TSE 

TSE<P> 

P 

P<T$E> 

TH2 

TH2<PTT> 

PTT 

PTT<TH2> 

XI 

o.os 

XF 

1.00 

ACC 

1.00 

Note: The notation A<B> indicates the value of A is 
read from the vapor pressure table as a 
function of B by linear interpolation. 
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quantities calculated in the course of the solution. All data are labeled with 
units (B/FDH and B/F2DH are Btu/ft^-oR-hr respectively). 

The input data are identified by their variable names which are defined in 
Table E-l. They are printed in the numerical order of their ID numbers. The 
initial data, describing the tank configuration, are recessed; the remaining 
data specify operating conditions. New input data in each case are identified 
by a leading asterisk. 

In the remaining sections, the single-value printouts are identified by their 
variable names used in the code. The double-value printouts are for regions 1 
and 2 respectively and are identified by the common heading of their variable 
names, with the terminal number removed (i.e., El denotes the variables ELI and 
EL2). The last line is an exception, where DPFD labels the frictional and 
dynamic pressure lrop in the tube. These outputs are defined in Table E-3 in 
their order of appearance. 
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Table E-3 
OUTPUT VARIABLES 


Label Definition 

EMD Total mass flow rate In all tubes 

THEX Temperature of two-phase fluid In heat exchanger tube 

PHEX Pressure in heat exchanger tube 

QEXT Extraneous heat input to tank through top 

QR Radiation heat flux to wall 

EL Flow regime tube length 

QT Heat transfer rate to tube 

TS Insulation surface temperature 

TC Minimum tank wall temperature 

TM Maximum tank wall temperature 

HE Heat transfer coefficient at external wall/insulation surface 

HF Heat transfer coefficient at internal wall surface 

HT Heat transfer coefficient at internal tube wall 

U Combined external heat transfer coefficient 

TA Characteristics temperature 

EN Dimensionless parameter 

DP F D Frictional and dynamic pressure drop in the tube 

Note: Labels EL through EN apply to two output data for the two flow 
regime regions. The variable names in the code for these two 
data are formed by adding 1 and 2 to the label. 
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Appendix F 

I DU/TVS PERFORMANCE COMPUTER CODE LISTING 

A listing of the main program, followed by the subroutines in alphabetical 
order, is given on the following pages. Preceding the listing is the sub- 
routine calling sequence in Table F-l and a glossary of major variables (those 
appearing in C0MM0N) in Table F-2. A printout of main program is shown in 
Figure F-2. 
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Table F-l 

SUBROUTINE RELATIONSHIPS IN PROGRAM STRUCTURE 


MAIN 

INPUT 

U 

PTABLE 

EXTHT 

FRECTT 

GENTAB 


RADHT 

H0UTER 

GENTAB 

FREC0N 

HINNER 

GENTAB 

FRECON 

HTUBE 

GENTAB 

OUTPUT 

U 

W12 

PDROP 
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Table F-2 (Page 1 of 3) 
GLOSSARY OF COMMON BLOCK VARIABLES 


Block 01 

General Variables 

TE 

External temperature 

TH2 

Internal temperature 

T 

Tank wall thickness 

11 

Tank wall Insulation thickness 

DO 

Tube spacing on tank wall 

EL 

Single tube length 

DTUBE 

Tube outside diameter 

P 

External pressure 

TSE 

Saturation temperature at external pressure 

QR 

Radiation heat flux at tank wall 

QEXT 

Extraneous heat transfer rate through tank top 

EMD 

Total mass flow rate in all tubes 

THEX 

Fluid temperature inside tube 

Y 

Number of tubes 

TTM 

Tank top thickness 

TTI 

Tank top insulation thickness 

ATT 

Area of tank top 

AKIT 

Insulation conductivity at tank top 

AKTT 

Tank metal conductivity at tank top 

PTT 

Internal pressure 

HTK 

Tank height 

ACC 

Acceleration 

XI 

Initial quality in tube flow 

XC 

Transition quality in tube flo.-; 

XF 

Exit quality in tube flow 

DI 

Tube inside diameter 

PHEX 

Pressure inside tube 

CHI FAC 

Average value of chi-function used to evaluate F in 
subroutine HTUBE 

RETA 

Reynold's number in annular flow region of tube 
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Table F-2 (Page 2 of 3) 
GLOSSARY OF COMMON BLOCK VARIABLES 


Block 02 

Region 1 and 2 Variables 

Variables for Regions 1 and 2 are denoted by their terminal numbers. 
These variables are listed below in the form Vx to indicate that VI and 
V2 are the variable names for Regions 1 and 2. 

AKIx 

Insulation conductivity 

AKTx 

Tank wall conductivity 

HEx 

External heat transfer coefficient 

Ux 

Combined external heat transfer coefficient 

HFx 

Internal heat transfer coefficient 

HTx 

Tube flow internal heat transfer coefficient 

ENx 

Dimensionless variable, N 

TAx 

Characteristic temperature, T/\ 

TCx 

Minimum wall temperature 

TMx 

Maximum wall temperature 

ELx 

Flow regime tube length 

QTx 

Heat transfer rate 

TSx 

Insulation surface temperature 

Block 03 

Properties 

VPR(20) 

Vapor pressure 

TSAT(20) 

Saturation temperatures 

NVPR 

Number of values in vapor pressure table 

THV(20) 

Temperature for heat of vaporization 

HV (20) 

Heat of vaporization 

NHV 

Number of values in heat of vaporization table 

TDN(20) 

Temperatures for liquid density 

DN(20) 

Liquid density 

NDN 

Number of values in liquid density table 

NCC 

Number of values in gas conductivity table 

CS(10) 

Gas conductivity 

CT (10) 

Temperatures for gas conductivity 

NVV 

Number of values in gas viscosity table 


1C3 
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Table F-2 (Page 3 of 3) 
GLOSSARY OF COMMON BLOCK VARiABLES 


VS(10) 
VT(IO) 
CP( 1 20) 
TCP (20) 
PCP(6) 
NCPP 
NCPT 
CPI (20) 
CP2(20) 


Gas viscosity 

Temperatures for gas viscosity 

Gas specific heat, function of pressure .1 temper 

Temperatures for gas specific heat 

Pressures for gas specific heat 

Number of pressures In gas specific heat table 

Number of temperatures In gas specific heat table 

Specific heat table at P, function of temperature 

Specific heat table at PTT, function of temperature 
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FF0G{JAM'Knf:TIWnT70ti^U , r.TArr?»!WPUT l TAFE6*nUTFUT) 

CC-MCO^ / 01 / Tt,TK2,T,Tl #Cn,r.t,DTUHfe,P,lSE,UR | f5fcXT,EMD.TW|X,Y 
? ,TT^,TT!,ATT,AKIT l AKTT,f , TT#KTK, ACC,Xl,XC,XF ,D! »PH6 X,CH|FACiHETA 
COMfeCN / TJT '/ A'Kn,ARTr.hFi,Ul,wFl,KTi # EM|TAl,TCl,TMl l EU,QTl,TSl 
*. ,AKU , ,AKT2.Mt2.U2i''F2,hl^,FN2 ( TA? l T(;2|TM2,ei.2,0T2,TSZ 

UOMPCN / C3 / VP«(20>iTSaT(20)|NVFR,T < V(ZO) # HV(?U),NmV # 

S TtN<2C)iDN<20)#^DI J #‘iCC,CSaf.'),CT(10)»^vV,VS(in),VT(io> 

* ,lf tl2C).TCP(20),rcP(6)iNrPK,r!CFT.CPi(20)iCP2(20) 
lc FC'RPAT { J5P:iTbHATIo.M LO^P CCP r LE TFD • t'-WOH «,b9,2 ) 
lul FCRPAT ( iHi ) 
i)«Mb(\& JCN V ( 3 ' 

^ R I T b ( 6,101 ) 

UC CGNTjM'f 
CALL INPUT 

fet « 3,14159 . DTU&F/12, 

CO lo« 1 « 1,10 

c*»«* F 1 0* PEOlPE Ti.Eb l E f ’u T P 7 ANT f PRCK CP* CP AkE CALCULATED ON EACH 
C PASb Tr'POLGP T ML OUTER LC'T tr.fe R CR CPPCK STARTS UN TplRfl PaS$> 

M * P Ay 0 ( 6 / 1 # 1 ) 

;c ire ¥ ■ i,kl 

n-MPLPATU'-ES, PRUPLRTlfcS, HEAT YPANbFPW ClfcFF ! C I b‘iTS AND RATfcS 

arf c Alcul at^p oh each pa^s Tv-much t.i inner loop 

«* tVALuATfc MEAT TRANSFER C&LFF I C l P N T $ A >C b 3 UATJJN TtRvE 
IF ( T I ,lt. % ) GO Tc He 
CALL POlTFR ( Pbl, TSl # T^i# <F-. ) 

CALL ►•Cl'TF R ( *L2, TS2. T r ’2* ► * '> ) 

v«0 Ft 111 

115 SCNTiMjL 

CALL PCI TU ( Pbl, TCI, T f: l # M 1 > 

CALL POtTFR ( PP 2 , TCP, T'2» *F 2 ) 

ill CONTINU 

Ul « PEl*AKll / ( AK I 1 ♦ Mfci • TI/12, ) 

l? * Pt 2 • AK 1 2 / ( AK J r ♦ HEP * T l / 1 ? , ) 

CALL P^NFR ( PF 1 , TCI, T^l ) 

CALL P I f NFS ( PF 2, TC2, T M 2 ) 

TT1 » Th 

IF ( KFi ,FC. 1 ) TTj, * T5E 

tAl a ( GP' ♦ i-F 1* Th2* l ' 1 • T T 1 > / I HF 1 ♦ LI ) 

TT2 * TP 

IF ( *VZ ,EC. 1 ) TT2 « T*E 

TA2 * ( UP ♦ PF2*TH?* u?*TT?) / < PF2 a L? > 

bN-1 * Sf.RT ( ( HfJ ♦ Hi J / ( A*Tl * T/i2. > ) * 00^24 # 

Ef-2 « SCRT ( ( PF? ♦ 1J2 ) / < AX 7 2 * T/i2, ) ) • DO/24, 

TANPM * T A \ P ( t M > 

TANp\a s TAPp(E fJ Z) 
bl * LJJ/12 , * f TTl- TP2 > 

tP » CO/12. * ( TT2- Tm2 ) 

b3 * 1, ♦ LI / HF 1 
£4 « 1, ♦ U2 / HF2 

IF < £Pr 0, > OTl 8 Y,EL1*'J1*E1 ♦ '= ( j*utXT 

Call PTUPt f pTI* TCi, l ) 

CALL PTlbF. { PT 2, TC2, 2 ) 

A1 « TA1*(TH2,TA1)*EM/TaPHA 1-.11*F1 / (tO*HTl> 

A2 » TA2MTP2-TA* )*EF 2/TA h H K 2-L?*F2 / <EU*PT2) 
bl » £:;1 / (TaM'P1*mf 18DQ/15, ) ♦ £3 / (F 0*hTi) 

be « bN2 / (TaNHN 2* u F 2*Do/le» ) ♦ 64 / IPU*M2J 

IF I X ,gt, X ) G^ TO 107 REPRODUCIBILITY of the 

ORIGINAL PAGE IS POOR 


Figure F-l 
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C BEGIN CALCULATION Of FLO* RFGIHE TUBE LENGTHS 

— et~» xjjrsrn*** ^TWmXCT^ r r^ ('^~^7rr ~ 

Cl * «XC-XI> * U2 • El / P.5 

C2 ■ (XF-XC) * I'i • E2 / E5 

di * cot • ~s r tit rr'v * ei > 

;l? * LE*T • B2 / E5 * (X-XC) • 63 t t v • EL ) 

AA » 1, • (Ci*C2)» (H2*Bl) 

PF ■ Ci*(C2-BiJ - !?!• « Cr*“CT2» - ^i*r2) - t, 

CC * Ci • 61 ♦ Hi 
HADlC * b» f • • 2 * 4 »*Aa*GC 

if x apsuaf .LTr'AnsigTTi^/TA.tTn j *~ti,ooi > co to io:< 

ElX 8 { -BE - SOB T ( RADIO ) > / < 2,MA I 
go u io4 
It 3 ELX • C. 

*/!’ ly5 « 1 # *. 

ICb EUX « ( -CC - ELa«* 2«AA ) / BP 
10« jr f El* .IT, ,99«9 > aft To Iftft ' 

ELI * fcL 
ti . 2 « 

60 TO 1C? 

1 Ct> ELI * EL* • tL 

el 2 * Pi - fcu 

c • *. end calculation nr rinw rfgtnf tope lengths 

1*7 CLNljMlf 

C •* ElALLATF hfc AT TWA -SFrR WaT 6? fc‘.n SYSTl-N IEMPEPaTUHES 

21 * t ( A1-A2 > *V - B1/FL1*REXTJ • FL2 f (Bi*EL?/FLi ♦ BP) 
THE* * A l - f,i • » G I ♦flEXT > / < fc L 1 • V ) 

JT 1 a V8tH*lJ«El ♦ 

2T? s y*er 2*i?*F2 - H 3 ** fTI 

TCI * THLX ♦ cTl/<t3*FLl*"Tl * v ) 

TC2 a TCI 

IF t tl? , 3 t , 0# 1 T C? *“TH?X * CT2/t r 0*EL2»HT2*Y1 
TK1 8 TAl - { TA1 • TCl ) / C 0 S M ( F N 1 > 

Tf>2 * Th 2 * ( TA2 - TC2 ) / CtSHtEN?) 

if ( ti ,le; j'ro’td 1*6 ■ 

TW 1 8 ( TCl ♦ TH 1 > 8 c.5 

TS1 8 TH ♦ T I /12 « *l;l/A*< 1 1 • ( TE - T H ) 

TW2 * ( TC2 * TH2 ) * !J.5 

Tb2 8 TK> 4 T I /12 1 *U2/ AK 1 2 • I Tl • M ) 

1*r CdNTINviL 

C •• EFRCT5 CFECK ' 

L « I - << 1-1 )/**>•» 3 

*<L) 8 f-nr 

IF { l ,LT, 3 ) GO TO XD2 
A VV s ( V(l> ♦ V(^) ♦ V C 3 ) 1 / 3 > 

fcf- * <4F AXl<Y(i)#y<2># *M3>>-Ar , rumi>r/t2),V(3>>> / AVV 
IF C EP .IT, T s 5" T0 1!*9 

1:2 continue 

*PP6 ('6# 13 ) 

139 CALL ClTPUT 
UU ?c 100 
fc.O 


ms 


-nT) 
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BVtJtrr cm — 

CUN^C* f P3 / V^«(20)»T«;AT(?0)i l JVPR|T'.iV(kn) # HV(20)|NwV, 

5 UM20>,DM2Q>#‘*D'V:CC.FS<1.' ) . CT ( 1 0 ) , N VV , VS ( in ) » VT ( 10 ) 

T »Cp 1 120 ) ,'TcP (23 ) , PC p 1 6) •K r P p » MCPT iCJ'H 20 ) » CR2 f 20 ) 

C ••• yapcp r r<fcSsijPE ( CR - psU ) 

DATA K\*pr: / ! 4 / , ( TRat(I), I s 1 * 1 ** ) 

? / 25 ,» 27 . 5 ,?C.i 32 . 5 , 3 r >.. 37 . , 5 , 4 C,»*Z, 5 , 4 ’ 3 ,. 47 l < 5 l ' 5 Cl M 52 , 5 , 55 ,,!> 7 , 5 /, 

* ( VPH(l)# l # 1,14 > 

is /1,1»2,1,4, ,*,7, ll.» 17, . ?5, > ^4,, 46, b, 62,, 82,, 104, ,130 ,,160,/ 

C *** Hi AT Cp VAPCRI2*T10 m ( 0 M • I* TU/LB ) 

IjATA \HV / 1* /, ( TMV(i:, I * 1,1? ) 

' t <5t i 2C , , 20«* 22,, 37 % , 4t>,, 45,,- S> u t * 53,, 56 , , 56 , , 59 , , 55’ , 4/ , 

* t PV(I>, I s i,is ) 

5 /IV? , , j 95 ., iv6 ,, J 95 ,,191 , ,11*7 ,,175, .153,, 134,, 105 ,»70,, 40., 0, / 

C *** UCUt I' k ^ S T T V ( S'* - l.p/PT 3 ) 


DATA 

ITN 

/ 2 : 

/. ( 

TPN(I> 

. I 

8 1,20 ) 





W ?5. , 

31. , 

3b, , 

. 4 o, . 

44, 

. 45, • 

46 , , 

47, , 

*3. , 

49, , 

' u 

•• •* 

0. , 

M. , 

5k, , 

i 53, * 

84, 

. 55, , 

56, , 

^7, , 

36. » 

59, /, 

1 

( 

FM I ) 

1 I 8 

1 . 20 > 







S> / 4 

.**c. 

4.62, 

4,44 ( 

. a, 29. 

4, ■ 

c, 4,C2, 

3.97, 

3,91, 

3,35, 

3,73, 

5 3 

. 71, 

3,64, 

v$ f b 6 j 

i 3 , 4 fc , 

’.3 

9, 3,29, 

5,17, 

3,n2, 

<*,85. 

?,5 / 

5 AS 

CCfxfUCTlVlTY ( 

OR - RTL/rT-hfi-OK 

) 




DATA 

(CT( 

T ) • I 9 

1.3) / 

' 35., 1 

4"., 

240, /, 

KCC / 

3 /. 



* (C 

6(1). 

I al , 3 

W .* 

;Ch2. .0 

3? 7 . 

,1644 / 






c *** (i AS VISCOSITY ( :R • lP/FT-HK ) 

DATA <VT< 1), 1*1, i: ) / 3T.,4n,,A5,,52,,fti.,i71,,o?,,10l,,l3G,,2no, / 
* * <VS(U, 1 = 1,13) / , 00216 , ,0i270, ,i:0«Cfc» ,00360. tC04x*» .0C4B5. 

5 ,C057f, ,0Cfc47, ,;C701. .OH" / , NW / 10 / 

C ***■ IjAS S F !r C 1 1 I C hE Af ( nR.pSji . HTL/Uh-gR ) 

JATA ly / 2, 55, 50,2, 49,6*2, 48, 7, 49,2,51,7,55.2,60,2,77.2,97 
1 , 3, 21, C.. 2*2. 7, 67, 2, 6 r ,?. '6, 2, 54,; ,^3,2,5?, ?, ^1,3*2. 50, 2, Si, 

7 2, 53,2 / 6, 2, 61, 2, 77,?, 9 *>,3, 21,0, »?*"', 96, 2, 7«s, 7, 66 »?, 6?, 2, r '8, 2, 5^ 

3 , 2,55 , 2,5 4, 2, 53, 7*7,52, 2, 55, r ,*6,2, 67, 2 ,78, 2, 99, 3 ,22#Q,» «*3, ?B, 
a 2.9«. 2. F2, 2. 73, 7,63,2,64,2, a, 2. 5«, r l ‘57,2 l 56 f ?,57,2, 59, 2,64, 2,7V 

5 ,3 , 00, 3, ?3, C., 4 *4, 20, 3. 49, ?, 10, 2,9?, ; ,P2. 2, 7S, 7,69,2,65,2,63,2,^0 

6 ,2, 60,2.62,2.66, ?,!?!, 3, cl, -*,^4.0 ,,5* * , 30 , 3 . 5b, 3 , 16, 2 ,95, : , R4 , 2 . 77 

7 ,2,?*, 2, 68, 2.0*, 7, 63, 2.65,7,63,2,82, 3,02,3,05,0, /. TCP / 25, ,30. 
f #35,»4o.,45,.50 l ,5 t: ,,6r'.»65,i7r’,.75,iCO,»90,,100,»110,il?0,»l 4 O l # 
9 ItbiilbR.iC, /, TCP / 1,47. 7,35, 1-1, 7, 29,4,44,1. 56. 6 / 

DATA NCFP. ixf.pT / 6 . IV / 
fcl 2 
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"SUBP OUTTNF "EXTH1 

C***8 CALCULATE: EXTRANbCUS HEAT TRAILER THROUGH TANK TOP 

COMMON / Ol / Tb#TH2,T l Tl.DC»Cl.iCTUe6iP,TSE l QR # QI:XTiEMDiTHEX|Y 

" * iTT»7TTT7ITT7*K I T . AKTT , PTT , HT»T7ffUC7n »TC |TCr, DU PH6 X , CH I f a C » RE T A 
11 * IE 
TS * TH 2 

T 2 ■ T3 « T4‘«”ITI7T5y*7I " 

JO lcl I *l>i 

CALL P RfcCTT ( All CTl*T2)*,9i AHS(Tl*T2)i P ) 

A3 s -akttv 

CALL f&ECTT ( A4, <T4*T5)*,5, AHS(T4-T5)» PTT ) 

If ( TT 1 f 6T, 0. > GO TO 10? 

T? * {A3*»4*T9**r*ri3*14y*m t {Ai*A3*Al*A4*A3*A4) 

T 3 » 12 
lit' 1U 103 

1?2 CPKTJMT ' 

A2 * Ak ! r ✓ TT I * 12, 

T2 8 (Ae*A3*A48T‘j*Al8(A28A3*AL 1 *A4*A38A4)8Tl) / 

* ( l Al*A2)*U2»A3»S7*A4§A3*8tn#A2*»28{A38A4) J 

T3 « ( A3<*A4*T«5*A2*(A3*A4)*T?> / ( A2* A3* A2*A4* A3*A4 ) 

13 CPWIMJF. 

T 4 8 CA4*T5*A3*T31 '7"TA3*A41 ' 

'»F XI * a3 * ( T3 - T4 ) * ATT 
Ki COMiNut- 

HFTLRK ’ ‘ 

CND 


S UB t* OCT 1 NT- -TRFtUN- r~fTVr,UTiXiRKU ,BFTiCP|nCN,VtS ) 

C***» CALCLLATfc F RfeF CONVECTION COEFFICIENT AT TANK HALL 

CrMNUl' / Cl / Tb»TH2|T,TliD0iFl,DTUFE,P,TSE,QR|QexT»EMP|THEX»V 
T *TTf'',TTI,ATYi»K!TVAI<'TT l PTT,PTK',ACC,XI,XC l XF l O!*PHEX,C^iFAC»PETA 
U R s hhC**2 * Bbl * ACC*1 1 246F *07 * DT * X**3 / V!S**2 
P»> * Lp • VIS / COH 
Hh s C,r5*CCN/X * ( GR • FR ) *• ,25 
ht-TLWN 
bl.D 
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st'BHjiTirm rprcrr ( mrrTfV'CT, Ft* r 

c***» CALCULATE FRF-.F CC^VKTION COEFFICIENT AT TANK TOP 

COPf'CN / 01 f Tt.THS.T.T^cr.FL.C’UFE.P.T JE»a«.1tXT.EMDiTNEX,Y 
r- »TTf'|TTI t ATf,AKITtA , <TY,PTT#FTK»ACC»Vl|Xn f XF#DT»PHEX»CWlFACiRETA 
CCMbLN / C3 / VFK(20 ).TSA t (?0),‘:VFH»T-V<20)»MV(?0)»NMV i 
v T t r. C ? , D r < 2 0 > • '!)?•. , UCC . ;”S < K ) , CT ( 1C > i » VV , VS ( 10 > # VT ( 1C ) 
i ,CP(i20),TCP<20),PCP(4),N'!Pr l NCf ; T,Cf , l<20),CF2(20) 

UATA Ik,!ViIC / i*\ / 
ijATA U L R / D.iHoJ / 

KfHT s i,CF • PG / TF 
L* t T (s * i, / rr 

call Gi? Tab t tf , 3 , i*, t.-^, cp>’, ncpt » 

CALL GFfTaF < Tf , CPFG, 3. TC, CT, C5, fJCC ) 

CALL t^TAfc < T* • VIKii. 3« IV. Vl, VS, f.VV ) 

PH s , Or •CCNf:*<Kn : **»2»hET r **rP r * ACC *1 ,?9ftc*07/m SG*CQNG ) <*CT ) * *,333 

K F T l K A 

cUD 


Si:0F CLT 1 NF. GF\TAK ( T, P, L* I# *• V, N ) 

C**** uLNtML TABLL FbA fLJMA" Ti-TfcSFrL.'T I C-. OH STfcP FUNCTION) 
i'lKtASiCN X ( 1 ) , Y ( 1 ) 

IF” ( T — V ( 1 ) > ICC, 104 , 

lor IF < I ,Ee, 1 ) G.' TO f l r A. 1C4, 1C 6 1 104, 106 ), L 
l 8 ! - 1 

IF ( T »• y c I > > l'O, 104, 1 '’3 
1:1 ir t 1 ,ec, f ) '*c < io4, im, icr*, 105 , iqa ), l 

1*1*1 

if ( t pp y < i ) ) 1 : 2 , ic4, in 
1:2 l * ! • 1 

v:i pf ic < ic 4 , ion lot,, icn i c >, l 

IT: A K * VC P 
Kr TLNN 

1* ✓ I 8 I • 1 

irt « b y c i ) ♦ n-yun / r<t un-Min • (V(i*i).v(i)) 

<“l tlnn 
cNC 


I 
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MM-H* ( »-f . T C , T! ) 

— M-A I - f -w /h& f 1 ^ 1 *’ — f-gi fe- Hf — — t-n- H fc 7- ~ H."1 <HH 1} €FHL^~) — 

, - wl t Tl , t.r,’,T,T! ,or. ti t .LTu, B^ .TSE f'J'?, ikXT.fcMu.TWEX.v 

• .1 TH t 'f| , 4 T ' , A" ! T, A | tT,i t T , . TK.AC'-, VI ,xf # X» »nt ,PHtX,CHlFAC,-’f^4 

/ Co / V^Hlt«)»TSATC9C)»*iVPf«.T <V t £'•) ,HV( ?0) 

" TLr,(2:).C^(S8)« l ?f ,F.-CC,CS<1' » .CTCICJ VSCin) » tfT(lfl) 

* »CP(1JC) .r’Cf*(6) .NCPr, !CCT,Ct l(2'i)»CP2(20> 

10 - v - i -h— / 

HTL a c.Cl 
1)>T a ir/?4» 

C **♦ i lGHfcH T6PH6hATUH!r 

TlF a (TM ♦ T>- 2 ) * .h 

CAUL ufcMAb ( TLF, V-LF. 3. !C, TEN, U '» ) 

0+4= — • — j~r H ■■*- 0 " r fr - 3 — * — ( — +trF — •> — 94H — 5 

w ISL 8 i,GC‘=-2 - J . 1 * ( T^F’ - <0. ) 

ClH * 6,?7fc-L ♦ /.:<9E*4 « ( TLF • 20. ) 

CALL F+^CCk ( HI, A3S<T*+«TH2>,l'i<T,rLF » -tTL , C^L . CONL , V ! SL ) 

C *** w nUtji Tc»PcPATl*Kt PPGION 
TLF s (TH.7 + TC) * .i> 

CAUL CHk - TA - D ( ■ Wf- -afcfr-arH Hh TE * , ■ ■ &-, ■ N??N - 4 

dL « 1.2 ♦ C.L7 * t TLF - ?0, ) 

Cl KL 8 (,37b-? + '/.^9E*4 * ( TLF • 30, ) 

''ISL a 3, rot-? - 1,11F%3 * ( TLF • 40. > 

CALL FmeC CM ( l ; 2,A{c > (T^2-Tc),rXT,£LF, k a ,CPL , CONL , V I SL ) 

C ** AV&kAtb h;,aT TF«fM>FFR COEFF ! C H.>.'T 

+4 ■ — { — >< 1 * + 2 — ) — * — j-5 

i» Tlf-K 

ti-r 




j^LlITY OP TSS 
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sub w o on w g wou w -t *irTftrTtt r~rr-7 

C*M* HEAT TRANSFER COEFFICIENT CN outer WAIL (GAS ONLY) 

COMMON / Ot / TE,TH2.TiTl*CO.EL#DTUee»P»TSE*OR#GEXT#FKU#THEX#Y 
$ .TTMiTTI# ATT#JKfT» AKTT»PTTiHTKi'ACC»Xr#yC*^r#DI #PR'EX»CHir*C»RFTA 
COMMON / 03 / VPR(20 )|TSaT( 20).NVFK»T) V<20)*HV<20).NMVi 
% 70N<20)»Dn(20)#nON»NCC.CS(10)iCT<10),NVV.VS(10)|VT(10) 

I .CPU20),TCP<20),PCP(6>iNCPP»NCrT,CFK20>»CP2{20) 

Data WOR / 0,1860 / 

DATA iH, ID, IC* IV, IK / 5 * 1 / 

KF « 0 ’ ' " ' 

TW * < TM ♦ TC ) • ,5 
IF ( TW ,G6, TSE > GO TO 101 
C ••• CONDENSATION MEAT TRANSFER 
KF 9 1 

XLF « DO/12. 

TLF « { TW ♦ TSE ) * ,S 

CALL GENTAB ( TLF, DLF# 3# ID* TDK, PL, NDN ) 

CALL GENTAB ( TLF, H» 3, JR, THV, HV, NHV ) 

CONL * 0,1?7E«4 ♦ 0,2056-6 * < TlF •• 20, > 

VISL ■ 8.34E-6 •» 3,075B«7 * < TLF - 40, ) 

HE « 3600. *( ACC*DLF**2*H*C0NL**3 / C4 , *XLP*VISL*(TSEnTW ) ) )**0,2? 
RETURN 

C FRFE CONVECTION MEAT TRANSFER 

101 CONTINUE 

IF ( (TM-TE) ,LT. 0, ) GO TO 102 
C *♦ FOR WALL REGION WHICH IS ABOVE EXTERNAL TEMPERATURE 
DXT • Do/24, 

TF • (TM+TE)*,5 

CALL GENTAB ( TF, CPG, 3» IK, TCP, CP}, NCPT ) 

CALL GENTAB C TF. CONQi 3# IC# CT, CS, NCC > 

CALL GENTA8 C TF, VJSG# 3» tV» VT, VS, NVV > 

BETG « 1,/TF 
RHOG « WOR * P / TF 

CALL FRECON ( Hj, (T m*TB>, DXT, RHOG, BETG, CPG, CPF8.VISG > 

C •• FOR WALL REGION WHICH IS BELOW EXTERNAL TEMPFRATURF 
TF * ( TE*TC ) • , 5 
DT 8 ABS<TE-TC) 

FF n 1, 

GO TO 103 

102 CONTINUE - 

C ** ENTIRE WALL COLDER THaN EXTERNAL TfcMPLRATURE 
DXT * MTK 

• 0. 
r’F 8 2, 

TF 8 ( T£ ♦ TW ) • ,5 
OT * TE • TW 

103 CONTINUE 

CALL GENTAB < TF, CPG, 3, IK. TCP, CP*, NCPT ) 

CALL GENTAB < TF, CONQ, 3» IC* CT, CS, NCC ) 

CALL GENTAB C TF, VISG, 3# IV» VT, VS, NVV ) 

BETG • 1,/TF 
fiHOG ■ WOR '* P / tP 

CALL FRECON ( H2.DT, DXT, RHOG, BETG. PPG, CONG, VISG ) 

C •• AVERAGE HEAT TRANSFER COEFFICIENT 
ME s ( HI ♦ FF*H2 > * ,5 
return 

end _ 


i 


101 



r 


subrout i n»r Htust i K - - 

e**«* heat transfer coefficient inside tube 

COMMON t 01 / TE.TH2.f#TI#eO.EL.DTUBE»P#TSE.OR,QFXT#IMD»THEX,Y 
S . TTM till, ATT »A* I r.AKTTTFTT, HTKTfflXYXl HOTJ XT, FT . FFEX VCKTFA C» RET A 
COMMON / 02 / AKliiAKTliHEl«Ul#HFl»HTi,bNl»TAl,TCS.|TMl»EU»QTl#TSi 
$ ,AKl2,AKTa.H62,u2#Hr2,HT2,fcN2#TAe#TC2,TM2,EL2.0T2,TSE 

COMMON f 03 / VPRI20),TSAT<20I.NVPR,TRV(2in,HV(20),NHV‘, 

$ TOm;20),DN(2O)#NDN#NCC,C5(10)#CT(3LO),NVV»VBUO).VT(1O) 

S .CPU20)»TCP(20)iPCP<6)fNCPP»NCFT r CH<20)»CP2l20) 

DATA IPtIW,miH,lC»IW « • 1 /'• FOF 7 071P80 V, PI t 3,14199 t 
CALL GENTA0 ( THEX. HTVT, 3, IH, TUV, HV. NHV ) 

EMD b ABS(OTl) / ( (XC«X( )*PTvT) 

call gentab ( thex. nso", T. tvtvt, vs.' nvv } 

IF ( K ,EQ, 2 ) GO TO 101 
C m» REGION 1 • ANNULAR FLOW 
C •• FORCED CONVECTION COMPONENT 

CPL * 1,2 ♦ 0,07 • ( TWEX • 20. ) 

CONL • 6,37E*2 * 7,3*6M • ( TMEX • 30, ) 

VISL » 3.00E-2 • 1 ,11E«3 * ( THEX * 40, ) 

PR a CPL • VISL f CON|. 

BETA 8 48, • END / ( PI • VISL * DI * Y ) 

HC s 0,276 • CONL / DI ♦ RETA**0,fl * PRM0.4 
C *• BOILING COMPONENT 

CALL GENTAB ( THEX. PHEX, 3, IP, TSAT, VPR, NVPR ) 

RHOV ■ WOR • PHEX / THEX 

CALL GENTAB ( TC , PSW, 3. !W» TSAT, VpH, NVPR ) 

DP a ABS( PSW • PHEX ) • 144. 

CALL GENTAB ( THEX* RHOL# 3# IR. TPN. DN, NDN ) 

DT s ABS( TC • THEX ) 

GE • ACC * 1.296E«07 

SIG a 1,6E*04 • O.D6E.04 • ( THEX - 32, i 

PRCOEF « < RH0L/«W0R*PH6X/THEX) )*»0,P * ( VIBC7V ISL ) , 1 

F a C l, ♦ 2,29*CHlFAC*PRCCEF >**0.§ 

Sal, / I 1. * ,305®«09 * ( F**I.W * BETA )••&.£? > 

HB » S* ,00122bCONL*b» 7 9 *CPL** «45*PH0L** ,49/V!bL**,29*(GE*0P**3/ 

% SIGb*2)**,25*<DT/HTvT/RHCV>b*,24 

C *• COMBINED HEaT TRANSFER COEFFICIENT FOR REGION 1 
Hr a HP * F * HC 
GO TO 102 

101 CONTINUE 

C *aa REGION 2 • DROPLET FLOW 

CPG a 2,0 * 0,0556 • ( THEX - 20. > 

CALL GENTAB ( THEX# CONG. 3. IC, CT, CS. NCC I 

PR a CPG • VISG / CONG 

RE * 48, a EMD / < PI • VlSG • DI * V ) 

HT a 0,276 a CONG 7 01 a RE**0,8 a FRaaQ.4 

102 CONTINUE 

HT a AMAX1 ( HT, 0.1 ) 

RETURN 

ENP 
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c«*«* input# initialize and print data 

COMMON / 01 / TE#TH2»TiTl»i:0,6L.riTUBE»P»TSE»QR.QEXT#fMD.THEX,Y 
S TTM , TT I , ATT”. AKI TpSKTTVFTT, HTX^ Arr#*nyC7XT7DT , PPEI .THTFACi RETA 
COMMON / 02 / AKIl,AKTl,HEl,lll,Hri,HTl,ENl,TAl,TCl#TMi,EU#QYl,TSl 
S » AK12* AKT2 1 HE2 »U2* Hr 2 »hTg»kN2» TA2»TC2 i TM 2*bl 2*QT2,T5g 

COMMON / 03 / VPR(20),TSAT(20)»NVrP|T^VI?in,HVr2CI)#NPV. 

$ tdn( 2 o)»dn( 20 ),ndn»ncc,cs(iO)#ct(io) ,Nvv»veuo> • V T C 10 ) 

$ •CP(120)»TCP(20)»PCP<6)iNCPP,NCPT.Cfl(2Q)*CP2(20) 

DIMENSION Vf23),-0r25) 

data ast, bla / ih*. ih /. v / 2 *>*o /, is / l / 
loi format < 13,ei2,; > 

113 FORMAT ( 43H0* * • * T«E A20V6 CASl DID NCT PUN *•**/) 
CHIFUN(A) * < A / ( 1 < • A ) >**0,9 

100 CONTINUE 

ksf s i 

i chi 8 o 

EMD s 0, 

DO 110 I * 1.10 
110 Oil) ■ 0. 

C READ DATA 
10« CONTINUE 

READ I 5*101 ) II. VV 

IF ( II ,60, 0 .AND, KSF ,E0. 1 ) STO| 

IP I II ,EQ, 0 ) 60 TO 103 
IF ( II ,EQ. 20 .OR, II .EC. 21 ) ICHI * 1 
VIII) « VV 
0(11) » 1, 

KSF 8 0 
00 TO 104 

114 CONTINUE 
WRITE I 6.H5 ) 

GO TO 100 

103 CONTINUE” " 

C •** ASSIGN VARIABLE VALUES 
DO a V ( 1 > 

EL i V|Z) 

Y a V ( 3 > 

DTUBF * V ( 4 ) 

TWTR ■ V(5) ~ 

DI a DTUBE - 2.*TWTH 
T s V ( 6 ) 

TI a V(T) 

ATT s V ( 8 ) 

TTM a V ( 9 ) 

rri'T vutir ” ' 

HTK 9 V(U> 

AKTT » VH2) 

AKT 1 a AKTT 
AKT2 4 AKT? 

AKIT a V ( 13 ) 

“ -”"AK!l • AKJT .... 

AKI2 ■ AKJT 
TE a V C 14 > 

TSE a V(l5) 

P a V f 16 ) 

TH2 a V ( 17 ) 
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PTT a V ( 18 ) 

XI a VC19I 

XC a V(20 ) 

XF a V(21) 

ACC a V(?2> * 32, i*' ~ 

c • * calculate or assign default values 
IF ( Y ,LE. 0, ) Y a 1. 

IF < TSE ,GT . 0, ) GO TO 102 
IF ( 0(16) ,GT i 0, ) 0(15) a 1. 

CALL PENTAB ( p, TSE. 3, IS# VpR, TSAT, NVPF ) 

102 CONTINUE 

IF ( TE ,LE, 0, ) TE 8 TSE 
IF ( 0(15) ,GT. 0, ) 0(14) a 1. 

IF ( P ,GT, 0, ) 60 TO 10? 

IF ( 0(15) , GT i 0, ) 0(16> > 1. 

call gentab ( tse, p, 3. is# tsat# vrp, nvpr ) 

107 CONTINUE 

IF ( TW2 , GT , 0, ) GO TO 106 


IF ( 

0(18) ,GT # o, 

) 0(17) 

a 1. 




call 

GENTAB ( PTT. 

TMg» 3# 

IS* 

VFF# 

tsat, 

NVPF ) 

106 continue 






IF ( 

PTT .GT. 0, ) 

GO TO 1 





IF ( 

0(17) ,GT, 0, 

) 0(16) 

« 1 . 




call 

GENTAB ( TH2, 

PTT. 3# 

IS. 

TSAT, 

VPR. 

NVPF ) 


100 CONTINUE 

IF ( XI ,cE. 0, ) XI a 0,05 
IF ( XF ,LE, 0. ) XF a 1. 

IF ( ACC ,LE, 0, ) ACC a 32, l6 
C *• mark new input values 
DO 111 I a i,22 
IF (0(D) 112#H2,113 

112 0(1) a BLA 
GO TO 111 

113 0(1) a AST 
111 CONTINUE 

C **• PRINT INPUT DATA 


Call 

w 

< 

6M DO 

$ 

DO 

. 6HIN 

9 

C(l) 


CALL 

w 

( 

6H EL 

9 

EL 

. 6MrT 

9 

0(2) 


call 

w 

( 

6H Y 

• 

Y 

t 6H 

9 

C<3) 


CALL 

Ml 

( 

6H DT'JBE, 

DTUBE 

# 6H l N 

9 

C ( 4 ) 


Call 

w 

( 

6M TWTH 


iTWTW 

# 6M IN 

9 

0(5) 


call 

w 

( 

ftH T 

f 

T 

• 6HIN 

9 

0(6) 


call 

w 

( 

6H T I 

1 

T I 

# 6h In 

9 

0(7) 


call 

m/ 

( 

6W ATT 

• 

ATT 

. 6HS0 FT 

9 

0(6) 


call 

Mi 

( 

6H TTM 

1 

TTM 

> 6MjN 

9 

0(V) 


call 

W 

( 

6H TTI 

1 

TTI 

# 6h In 

9 

0(10) 

) 

call 

Mi 

( 

6H HTK 

1 

HTK 

# 6WFT 

9 

cell) 

) 

call 

w 

{ 

6M AKT 

I 

akti 

# 6H9/FDF 

9 

0(12) 

) 

call 

Ml 

( 

6H AKI 

• 

AKI1 

# 6H8/FDH 

9 

0(13) 

) 

Call 

w 

< 

6HTE 

9 

TE 

# 6F0EG P 

9 

C(14) 

) 

call 

w 

( 

6HTSE 

9 

TSE 

# 6HDEG P 

9 

0(15) 

) 

Call 

«l 

( 

6HP 

9 

P 

# 6PPSIA 

9 

0(16) 

) 

call 

W 

( 

6HTH2 

9 

TH2 

. 6H0FG F 

9 

0(17) 

) 

call 

w 

( 

6HPTT 

9 

ptt 

# 6FPSIA 

9 

0(18) 

) 

call 

Mi 

( 

6WXI 

9 

XI 

, 6H 

9 

Q(19) 

) 

call 

Mi 

< 

6HXC 

• 

XC 

» 6W 

9 

0(20) 

) 

call 

Ml 

( 

6HXF 

1 

XF 

# 6h 

9 

0(21) 

) 

call 

Ul 

t 

3HACC# 

ACC/32,16 

# 6PG 

9 

0(22) 

) 
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IF < TE .16. TW2 ) QO TO 114 

tft xr ; Ler^r T tm itJ -n « ~ * ~~ 

C ASSIGN INITIAL VALUES TO START ITERATION 

TMl « 1,01 * TH2 __ 

tm 2 • m " " ' ~ ~~ 

TCI ■ 0,99 * TH2 
TC2 * TCI 
THEX 9 TCI • 1,0 
T§1 s 0, 

TS2 » 0, 

IF ( T1 ,LE. 0, ) GO TO 109 
TSl s ( TSE ♦ Tg ) • 0.5 
TS2 « TSl 
109 CONTINUE 

Eli 8 amaxk xc**3, o,i > 

EL2 « EL * ELI 

C CALCULATE OTHER QUANTITIES USED IN SOLUTION 
CALL PTABLE ( P, CPI ) 

CALL PTABLE ( PTT, CP2 ) 

CALL EXTHT 
CALL RADHT 

IF ( ICHI ,EQ. 0 ) RETURN 
H 8 ( VC » XI ) / 100, 

X b XI 

YSUM * 0,5 • CHIFUN(X) 

DO 105 I a 1,99 
X s X * H 

105 YSUM 8 YSUM * CMIFUN(X> 

IF < XC .LT. 1.) YSUM s ySLM * CHIFUN(XC)*0,5 

chip ac * ysum / ioo, 

RETURN 

END 
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SURHLTINF OUTPUT 
C***» H I f' T SCLLT I O' UATA 

wU^^ON / Oi / TF t T|H?» T» T 1 » DTi FLi C^URfc |P» TSfci QR, OfeXTifcHQi TMtX » Y 
T #TTF,TT!,ATT#AKIT ( AKTT,PTT,HTK|ACC»yi,*R,Xr # Dl»PHEX.PMirAn.»ETA 
/ 02 / AKU. ( 4KTl l MFl l Cl l HFl,FTl,r-M,TAj l TCl#THl,F.U,CTl,TSl 
" .AM?# AKT2«Hi:2iLi » <F8 »hT>ifcN2,TA? l TC2»TM2 # lrl.2,cTrf l TS'* 

?e r CRT at { 1H ) 

‘U m/rfat ( ••#/< *h— .-•)/ > 

( 6 . ) 


1 

* 

CALI 

k ( 

6FEh 

. FrP . 

AHI.F/HA , 



I 

CALL 

A ( 

AFTl-cX 

. TCFX . 

AHlfcO fi , 

1H 



CALL 

k ( 

F h F k> F A 

* RFEX . 

AhfblA , 

iH 


1 

CALL 

k ( 

Af.CtxT 

» CCXT i 

AHTTif/kR, 

1H 


1 

» 

call 

k ( 


. 3k , 

AhFT'l/hR, 

1> 


i 

CALI 

MS 

( 3KI » 

CL1. i ; 

L* * 6hM 

) 


f 

CALL 

*1? 

( ?hCU 

TT1, 0T2. 

AkiPTU/FP 

) 



CALL 

Mr 

( 2hlS» 

ISl, Tl>2» 

?hj hr. R ) 




CALL 

M? 

( 

3hC 

• 

^ * < 

1 

TC»* 

6HLTG H 

) 

CALL 

k i 2 

( 

3UTK 

9 

TM1 1 

T U ? » 

6MCFG K 

) 

.HU 

( 

9. 

Q vi ) 






U^LL 

M2 

( 

3 H h L 

9 

lJ r 1 • 

"fc?» 

6Hb/F?Du 

) 

CALL 

"1? 

( 

3M-f 

1 

MF1. 

nr?» 

6:iE/P?Dr 

) 

CALL 

•'ll 

( 

3f.LT 

# 

HT1, 

ut?» 

6UB/F?Uh 

) 

CALL 

M2 

< 

3k 

» 

Ul. 

u?» 

fcHd/F2tf.< 

) 

CALL 

^1? 

( 

•5LTA 

1 

TA1, 

Ta" * 

(iWCFG R 

) 

CALL 

"I? 

( 

•Ilc:. 

1 

EM. 

t K?f 


) 


CALL FI’MIP 
-MTE ( f, 90 ) 
f* L T L h ^ 
r.i-D 


) 



,^ 1 . 1 -iTV OFT?!" 

OWG IS 
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SUPROUTfNfc PDROP 

C**»* CALCULATES PRESSURE DROP in the TUBE (FUR OUTPUT ONLY) 

common / 01 / T6,TH2»T.TI#CO.eL.DTUB6#P»T6E*ORiOEXT»FMU»THEX#Y 
$ ‘.TTM.TTI V T TT; A K ITyJRTT , FT TTWTK7A"CC75CI » VC # XF , D I , PBEXyCHTFffC »RETA 
COMMON / 02 / AK!l.AKUiHEl,l4#HFl,hTl,fcNl.TAi,TCliTMi,iHiQU,TSl 
$ ,AKl2,AKT2,HE2,i!2,Hr2,pT2,fcN2.TA8,TC2,TM2,EL2»0T2»TSE 

COMMON / 03 / VPR(20>tTSAT(20),NVrP»TrV(20>»HV<20r,NHV# 

S TDN(20>#DN(20),NDN#NCC,CS(10).CT(10),NVV,VS(10),VT(10) 

* «CP( 120 ) i TCP < 20 ) #PCR(6 ) • NCPP* NCFT » Cpl( 20 ) »CP2(20 ) 

DATA ID / i / “ 

RHOV * 0 1 1860 * PHEX / THEX 

CALL GENTAB ( THEXi RhOL* 3. ID* TpN, DN, NDH ) 

PL b 64, / RETA 

IF ( RETA ,GT, 3000, ) Fl » 0.316 / R{.TA**0.25 
E ■ (4,*EM0/Y/(3 . l4l59* (Dl/1?«)*»2) )**2 / < ACC*1 , 296E*07 ) 

A > E * FL / ( 2, * RHOL • DI/12, ) 

B b A * ( RMOL/RHOV * 1, ) 

C ■ *6 / ( RHOV * PHEX*144, ) 

D b ALOG ( 1. ♦ C • (XC«Xl) / <l,*C*XI) ) 

DELPF * EL/(C*(XC»XI )) • <<A*B*Xn*DH/C*<C*(XC*'X!)-<l.*C*Xl>*D)> 

delpp « PHEXMRHOV/RHOL-1.) • D 

CALI W12 < 6«D P F D, DELPF/144., DFL p D* 4FPSIA ) 

RETURN 

End 


SUB* CUT INF PTASLh ( P, Zf X ) 

C***» t N fc N * T P A TaPLE CF SPECIf \r *£*7 AS A FUNCTION OF T p M P f H A T Ij R f AT 
r i f >Pl T PRESSURE p A T"0 «VaR JARLE (TfMRLRATjPc,pHtSSLMF) TApLb 

CCrFC* / C? / VPR(2.1).TSA T (?0) l , iVFR,TMV(2f’),HV(?0),M,V, 

* TCM20.DM2C)*’ DM|NCC,CS(lt )*CT(19)|NVV,VS(1Q)* VT<1Q) 

.Lp(le’O) , TCF(20),rcP(6) iNCPl , ‘ICF T , C»’l( 20 ) , CP2 (20 ) 

DIMENSION UX(29»»),CPX(l) 

MUl vALENCb ( CP, OX 1 
P1F » i, 

DC 1C* I b 2, \CPP 
II = I 

IF ( F C F ( I ) ,LT, P ) CiQ T( 10!' 
p if a t p - prru-i) ) / ( pcr(n - ffp(i-d ) 

IF ( PCF<3 ) ,GE. P ) P 1 F s (*, 
uC. TC xP6 
1C5 CPNT JMJF 
1C6 CUNTIK'L 

l»C K<? I * l.NCPT 

107 CPXtp S rx<I,Il* 1 .) + PIF * ( CX ( 1 1 f I } - CX ( I * I T "1 ) ) 

HhTLRN 

Li'D 
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Sl;6*CCTINr HAI*WT 

CALCULATE R AC 1 AT I '-N >*MT TRANr.fPN 

/ n / Tt-#TH?,T f Tl*Lr!#tl.,CTUF-fc # P # 1Sfc # 0P # 0fcXT.t;V.D#THEX.Y 
a •TT^.TTI,ATTi*K!T l AKrT,rTT,»-TK#*CC.y|»XC,XF iDI'PHEX.C^irACiRETA 
,f> » c, 

*1 TUN 
it D 


SirwctTirjr Ir ( t, c, ii, r , 

(.'•••* *» K 1 1 tl !- 1 N (* L t CLH= T \ i ! > i ~i t* K 

i,r KRMt ( 4X,Ai,Afc,3H s V,-J,;-X,A6 ) 
i i FTfi^AT ( 4\ t Al i ft • 3>: * • i u • ? » i: X , A^ } 

A » Ah', ( C ) 

n i a , l t, 3,ji , n w, a jcto, ) -.v, 13? 

'- ! nTb t 0,133 ) 0, T, !, 

*ETUN 
l‘lr CrM JNU 

'^ITE ( 5,1C1 ) 0, T, 2, t: 

Tlnf> 

u n 


r ^O^ c ®ff^ r %v,> ( 
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1 

ft 


I 


1 


1 


I 


subr o utin e- wirt -r r Pir- p ftn r r 

Cmm WR|Tg TWO OUTPUT NUMBERS ItSUAUV H6G10N 1 AND ?) 

101 FORMAT C 5X,A6»3H » ,F9,4,2X,F9.4,2X#A9 > 

102 FORMAT ( 3X, Hit » E^ i 2J7% i F9Y4775T# A? 1 

103 FORMAT ( 5X,A6,3M ■ .F9,4,2X,E9,e.BX,A6 ) 

104 FORMAT ( 5X,A6,3M ■ . E9 , 2 , 2 X , 69 , 2 , EX ,A6 ) 

N ■ 1 

A • ABS ( Dl ) 

IF ( A ,UT» 0.01 .OR . A ,Qi. 100P. ) N 8 N ♦ 1 

- g u ^aBS t &2 J ^ ~ ^ ^ 

IF ( 8 ,'wT, 0,01 ,OR. B ,GE, 1000. > N 8 N ♦ 2 

00 TO ( 105, 106, 107, 108 ), N 

105 WRITE ( 6,101 ) T,n#02» U 
RETURN 

106 WRITE ( 6,102) T,Di,D2, U 
REtUPN " 

107 WR I V* < 6,103) T,D1,"2» U 
RETURN 

108 WRITE t 6,104) T, 01,02, U 
RE’ URN 

ENp 
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Append'* G 

I DU/TVS PERFORMANCE COMPUTER CODE TEST CASES 

A listing of the Input deck for the test cases is given In Table G-l . 

Definitions are given In Table E-l . A complete set of Input data Is given in 
the first case; each subsequent case Includes only the data to be changed from 
the previous case. Each rase Is terminated by a blank card; an additional 
blank card ends the run. The end-of-file note on the listing is for clarification 
only and is not a part of the deck. 

In these cases, the external temperature is varied over the values 40, 45, 50, 

60, 80, 100 °R for two intern: 1 pressures, 20 and 25 psla. For the conditions 
of the last case of this series (100 °R external and 25 psia internal), the 
transition quality in the tube flow is given values of 1.0 (no mist flow 
regime), 0.05 (no annular flow regime), and 0.0501 (negligible annular flow 
regime) . 

The tank configuration consists of two wall-mounted tubes, each 35 feet in 
length and spaced 2.5 inches apart on the tank. The tube outside diameter is 
0.25 inch and wall thickness is 0.032 inch. The tank height is 2.8 feet and 
the area of the top is 4.4 square feet. The wall thickness is 0.5 inch and 
the top thickness is 1.0 inch. Insulation thicknesses are 0.25 and 0.5 inch on 
both wall and top. Conductivities of the wall and insulation are 90 and 
0.01 Btu/ft-°R-hr, respectively. 

For all cases, the external pressure is 15 psia; the external saturation tem- 
perature is taken from the vapor pressure table. The internal Lh^ temperature 
is also taken from the vapor pressure table, at the internal pressure. The 
transition quality in the two-phase tube flow is 0.85 for the first 12 cases. 
Default values were used for the initial and final quality in the tube and the 
acceleration level. 

The output from the above cases is given in Table 2; cases I through 13 are for 
the 0.25 inch insulation, and cases 14 through 28 are for the 0.5 inch insulation. 
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Table G-1 

TEST CASE INPUT DATA LISTING 


1 

3.5 

14 

40.0 

2 

35.0 

14 

45.0 

3 

2.0 

14 

50.0 

4 

0.25 

14 

60.0 

5 

0.032 

14 

80.0 

6 

0.5 

14 

100.0 

7 

0.25 

18 

25.0 

8 

4.4 

14 

40.0 

9 

1.0 

14 

45.0 

10 

0.25 

14 

50.0 

11 

2.8 

14 

60.0 

12 

90.0 

14 

80.0 

13 

0.01 

14 

100.0 

16 

15.0 

20 

1.0 

20 

0.85 

20 

0.05 

18 

20.0 

20 

0.0501 


END OF FILE 
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; t . Table 0-2 (continued) 

IDO/TVS PERFORMANCE COMWTER CODE RESULTS 
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Table G-2 (continued) 

I DU/TVS PERFORMANCE COMPUTER CODE RESULTS 

CASE 13 
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Table 0*8 (continued) 
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Appendix H 

I DU BREAKDOWN LEVEL AND FLOW LOSS PROGRAM 
H.1 Breakdown Level 

Sample oscillograph data obtained on liquid hydrogen breakdown is given in 
Figures H-l to H-10; the figures are annotated to illustrate liquid levels, 
the breakdown points, flowrate, valve actuation, pressure difference between 
the main tank and IDU, and temperature T4. Additional temperature data 
for test J-39 are shown in Figure H-9. It is interesting to note that 
breakdown in Test J-39 is closely correlated with the temperature increase 
in the ullage due to the warm, incoming GHo, and occurs approximately five 
seconds after the top of the LHAD is exposed to the warm ullage gas. 

However, a more precise correlation between '•.he ullage gas temperature 
increases and the breakdown point could not be obtained since the Dymec 
system sampled 33 data channels in approximately 35 seconds. Therefore, 
a given data point would be sampled every 35 seconds; within this time 
gap, more accurate curves of i^lage temperature could not be obtained. 
Further, the actual maximum temperatures may not have been recorded. 

Breakdown in all cases was determined by examining the continuous level 
sensor voltage vs. time, and then determining the height vs. time by use 
of Figure 4-1. Flowrates were determined by calculating the corresponding 
flow area in the IDU, using the equations given in Appendix I. It should 
be reiterated that due to the strong possibility of leakage out of the IDU 
through PV2, the actual flowrate through the LHAD and out through PV1 may 
have been less than measured Inus, the discrepancy in agreement between 
actual and predicted breakdown heights vs. flowrates shown by Tests J-ll, 
J-19, J-25A, J-27B, and N-ll, may have been partially due to failure of 
PV1 to open properly and PV2 to close properly. 

H.2 Flow Loss Program 

The flow loss program used to design the IDU LHAD, select the screen mesh, 
and correlate breakdovn data is given in Tables H-l and H-2. 
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Figure H-6 
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Table H-l 

I DU BREAKDOWN LEVEL AND FLOW LOSS PROGRAM 

NOMENCLATURE 

D3 * Outer Screen Diameter (Inches) 

T * Annulus Width (Inches) 

LI s Screen Height (Inches) 

F s Friction Factor For Flow Within Annulus 
Q = Flowrate (Lbm/Sec) 

D2 = Liquid Density (Lbm/Ft 3 ) 

D * Screen Pore Diameter (Ft) 

B * Screen Thickness (Ft) 

A » Surface Area to Unit Volume Ratio of Screen Wire (Ft" 1 ) 

E = Screen Volume Void Fraction 

N1 = Multiplying Factor For Pleating 

N2 = Number of Screen Layers 

U - Viscosity (Lb/Ft- Sec) 

PROGRAM 

10 INPUT D3s Ts L 1 
20 F * .025 
30 INPUT Qj D2 
40 INPUT Cs Bs A* E 
50 INPUT N1/N2 
60 Q*Q/D2 
70 U». 000008 

80 A4 » 3. 14 * < D3 * T - T * T > / 144 

90 R1 « A4 / ( 6. 28 * D3 - 6. 28 * T ) * 12 

100 D1 * 4 * PI 

110 F0R G * .3 T0 .9 5STEP .05 

120 X » LI * < 1-G ) 

130 A3 a 6. 28 * X * ( D3 - T ) / 144 * N2 

140 V » Q / A3 

150 N = V * D2 /U/A/A/D 

160 Fl 3 8.6/N + .52 

170 PI a N1 * 1.3 * FI * D2 *B * V * V / E /E/ D / 32.2 

18 C VI a Q / A4 

190 P3 » D2 * VI * D1 / U 

200 P2 * F * ( LI -X / 2) / 1 2 / Cl * C2 * VI * 1 / 2 / 32.2 

210 P3 a D2 * VI * VI / 2 / 32.2 

220 P4 « D2 * ( LI • X ) / 12 

230 P5 « P1 + P2+P3+P4 

240 PRINT X; FI j P21 P3; P4i P51 P3 

250 NEXT fi 

260 END 
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Table H-2 

PROGRAM OUTPUT 


CONDITIONS 

250 X 1370 Screen 

D3 = 24"; T * 1.5"; LI = 18.37"; Q = 2 Lb/Sec; D2 = 4.37 Lb/Ft 3 ’ 

D = .000026 Ft: B = .0004 Ft; A = 22773 Ft* 1 ; E * .2409; N1 » 1; N2 « 1 


SUBMERGED SCREEN DEPTH (INCHES) 


PI -FLOW LOSS THROUGH SCREEN (LB/FT^) 


P2- ANNUL US FLOW LOSS (LB/FT ; 


P3- DYNAMIC LOSS (LB/FTp 


P4-HYDR0STATIC HEAD (LB/ FT 2 ) 


P1+P2+P3+P4 
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Appendix I 

INTERFACE DEMONSTRATION UNIT AND MDAC 260-GALLON 
lh 2 TANK GEOMETRIES AND INSTALLATION CONFIGURATIONS 

Tne rate of outflow from, or inflow to, the IDU is derived from the liquid level 
rise rate as measured by the capacitance liquid level sensor. The IDU internal 
hardware and condition of the annular screen device (full or empty) determines 
the effective free surface area of the liquid, which must be known in order to 
determine volume flowrates of liquid hydrogen. The net free surface area as 
a function of height measured from the inside bottom of the IDU is given in 
Table 1-1. A.| corresponds to outflow with the liquid retained in the LHAD 

annular region. A 2 and A 3 correspond to specific regions for which the area is 
determined by the presence of hardware alone. A^ corresponds to refill or out- 
flow when the LHAD does not retain liquid above the nominal liquid level. 

Figure 1-1 depicts the layout geometry and associated volumes of the IDU and 
the MDAC 260-gallon tank. 
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Table 1-1 

FREE INTERFACE AREA OF IDU 


1. 0 £ Z £ 23 In. (prior to breakdown) 
A 1 " T t°T 2 “ (D S0 2 ' °SI 2) ‘ D CC 2 


A, = 498 in. 2 = 3.458 ft 2 


2. 23 £ Z £ 24 in. 


a 2 8 J *- d t 2 ' D cc 2 ' d pl 2 " d psl ‘ d tcc^ ' A C 
A 2 = 584.8 in. 2 = 3.811 ft 2 

3. Z > 24 in. 

A 3 = T ^°T 2 ’ D CC 2 ’ °PL 2 ’ D PSL 2] ’ A CP 
A, * 611.8 in. 2 = 4.25 ft 2 


°PL 2 ‘ °PSL 2 ' O ‘ A CP 


4. 0 £ Z < 23 in. (refill with gas venting from top of LHAD) 

A 4 » A 1 + T ' D S0 2 ‘ °SI 2) = 604.4 in. 2 = 4.197ft 2 

Nomenclature: 

Dj = Internal diameter of IDU tank = 28 in. 

Djq = External diameter of screen device = 24 in. 

D^j * Internal diameter of screen device - 21 in. 

Dqq = Diameter of central column = 3 in. 

Dp|^ = Diameter of pressurization line (with Teflon insulation) = 1.5 in. 

Dpsi = Diameter of pressure sensing line - 1.25 in. 

Dy = Effective diameter of instrumentation wiring bundles = 0.75 in. 

2 

A rD s Cross-sectional area of capacitance probe = 1.25 in. 
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Appendix 0 

PARKER 5640027 SUBMERGED LIQUID HYDROGEN SHUTOFF VALVE 
0PERA1 IONAL ANALYSIS 

J.l INTRODUCTION 

The Parker Submerged LH£ Shutoff Valve, shown in Figure U-l, is a 2-inch 
diameter, pressure actuated, ball type shutoff valve. The valve was designed 
for use inside a liquid hydrogen aircraft propellant tank. The following 
description is based on Parker Company Report SPD-132. 

Its features include: 

1. The design configuration which allows valve operation while completely 
submerged in liquid hydrogen, which lessens the LH£ boil off due tc 
external lines. 

2. The valve may easily be converted from a normally open to a normally 
closed configuration. 

3. Internal sealing is provided in either flow direction by a specially 
designed double ball seal. 

4. The valve can be pressure actuated for both opening and closing. 

The unit meets the following requirements (P/N 5640027) 

Fluids Liquid or gaseous hydrogen 

Inlet Pressure 0-50 psig 

Flow Capacity 2.8 psi pressure drop at 30 lb/sec LH 2 

Temperature -423 to +165°t 

Actuation Fluid Helium 

Actuation Fluid Pressure 500 to 750 psig 

Leakage 10 sc>m typical at -420°F 

J.2 OPERATION 

The valve consists of a main valve mechanism, an actuator, and a position 
switch assembly. Figure J-l shows two of these three elements for the valve 
in its normally closed position. The position switches, not shwon, monitor 
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Figure J-l . VALt/E SCHEMATIC 


valve positions by means of a cam located on the ball shaft. To open the 
main valve, GHe pressure is applied to the control port. The resulting pressure 
difference across the piston opens the valve against the force supplied by the 
return spring. The force supplied by tne actuator is transmitted to the main 
valve by means of a linkage which converts the linear piston motion to rota- 
tional motion of the ball. As the control port is pressurized, the small bleed 
orifice through the piston allows the pr?ssure in the cavity to the outside 
diameter of the bellows to increase. To close the valve again, control port 
one is vented and the spring, along with the residual pressure acting on the 
bottom of the piston, closes the valve. 

The following analysis was performed to determine the operating factors of 
safety of the Parker Submerged Liquid Hydrogen Ball valves. This analysis is 
based on a preliminary informal analysis obtained from Parker Aircraft Co.; 
operating pressures of 500 and 750 psig are considered. The operating factor 
of safety, FS, is defined as the ratio of applied actuation pressure to 
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minimum actuation pressure required to (1) start acceleration, (2) complete 
actuation, (3) jtart deactuation, and (4) complete deactuation. 

A summary of the analysis is presented below. The valve schematic and appro- 
priate dimensions are shown in Figure J-2. 

J.2.1 Pressure to Start Actuation 

A. With pressure drop across piston 


f + F k + F c + P 9 A. 

b-j s<| 2 b 


. ,76 p S , g 


FS 

_ 500 „ 

- m 

2.8 

FS 

_ 750 _ 

" rrs ' 

4.3 


B. No AP across piston 


f + F. + F c + P 9 A. 

b-j Zb 


• . 373 P sig 


FS 

500 

s 

1.3 

FS 

_ 750 
' 173 

2.0 
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Figure j-2 
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* Piston Diameter = 2.5C In 

U, * Effective Bellows Diameter * 1.72 In 
b 

Ap = Piston Area = 4.92 In^ 

= Effective Bellov/s Ar a • 2.32 In^ 

f = Seal Friction Force - 67C Lbs 

F^ = Dellov/s Force » 50 to 75 Lbs 

F s = Sprinq Force = 130 to 305 Lbs 

Pj * Pneumatic Pressure = 500/750 Psiq 

p£ * Line Pressure = 50 Psiq 

PARKER VALVE SCHEMATIC - ACTUATION CYLINDER 
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0.2.2 Pressure to Complete Actuation 
A. With AP across piston 


f + V F s 2 + Vb 


670 + 75 + 305 + 50 (2.32) 

02 


FS ■ |§7 = 2,1 


FS = 237 = 3,2 


B. No AP across piston 


f + F. + F, + P 0 h, 
b 2 s 2 2b 


670 + 75 + 305 + 50 (2.32 

02 


pc - 500 _ -I ft 

FS ■ W " 1,0 


pc — 750 _ ic 
FS " 3B? ' 1,5 


% 


232 


237 psig 


= 502 psig 
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J.2.3 Pressure to Start Deactuation 


A. Line pressure 50 psig; P 1 = 0 


I 




> 


670 - 75 - 305 - 50 (2.32) 

(4.92 - 2.3'2'j 


67 psig 


B. Line pressure - psig; P-j = 0 
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670 - 75 - 305 - 0 (2.32) 
4.92 - 2.32 


111.5 psig 


Note that calculated pressure P, is ( required to start valve to deactuate. 
Hence, valve will not move unless P, is greater thanK by the values noted 
above. Also valve will not fully deactuate unless P-J is greater than: 

A. Line pressure 50 psig; P-| = 0 


_ 670 + 50 - 130 - 50 (2.32) 
*1 ' 4.92 - 2.32 


182 psig 


B. Line pressure 0 psig; P 1 s 0 


_ 670 +50-130-0 (2.32) 
1 " 4.92 - 2.32 


227 psig 


J.3 CONCLUSIONS 

A. Successful deactuation of valve is dependent upon venting actuation 
pressure (P^) quickly enough such that deactuation pressure (P|) 
exceeds P.j by 227 psid. 

B. The following conditions may prevent this AP from developing: 

1. Restriction in pneumatic vent, i.e., obstruction in orifice or 

screen located in actuation pneumatic port, slow acting or sluggish 
solenoid. 
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2. Ineffective seal on actuator i.e., damaged seal, damaged cylinder 
surface, contamination, tapered cylinder. 
k C. The following conditions may require a higher AP to accomplish 

complete deactuation: 

% 1. Binding, i.e., contamination, incomplete temperature stabilization, 

excessive seal drag on ball. 
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J.4 RECOMMENDATIONS 

A. Provide a separate pneumatic "close" system for each valve. Reasons: 

1. Full pneumatic pressure will be available for deactuation, i.e., 

pc' = S.Q0 = y 2 

2. Leakage past damaged piston seal is not critical. 

3. Operation of actuator is not time dependent. 

4. Reliability of actuation increased greatly. 

B. Add second seal to actuator piston to reduce piston leakage when actua- 
tion (open) pressure is applied. 

C. Plug orifice in actuator piston to reduce piston leakage. 

D. Purging of the volume external to the valve ball should be accomplished 
either by partially opening of the valve during the vacuum/purge cycle, 
or by installing a separate purge system attached to the ball shaft boss. 
Partial opening of valves may be accomplished by applying a reduced 
pneumatic pressure to the open port. 

E. All pneumatic helium vents should be protected from breathing air by 
trapping and holding exhausted helium at vent opening. An inverted 
chamber (opening on the bottom) over a vent opening may be suitable, 
and/or a check valve should be used. 
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